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Random layer lattice structures are considered which consist of layers arranged parallel 
and equidistant, but random in translation parallel to the layer, and rotation about the normal. 
We call a and 6 the axes in the layer, and c the axis normal to the layer. In this notation there 
will be crystalline reflections of type (00/), two-dimensional lattice reflections of type (hk), 
and no general reflections (hk/). Equations are developed for the intensity distribution in a 
two-dimensional powder reflection, and for the integrated intensity. Equations are also de- 
developed for the particle size in terms of the peak breadth, and for the displacement of the 
peak. The powder pattern of a heat treated carbon black is presented as an illustration of 


two-dimensional lattice reflections. 


I. INTRODUCTION 


HE x-ray diffraction patterns of certain 

layer type materials indicate that they 
may exist in a form which is intermediate be- 
tween the amorphous and the crystalline states. 
For example, the x-ray patterns of certain heat 
treated carbon blacks indicate that they are 
built up from individual graphite layers arranged 
parallel to one another at about the normal 
graphite spacing but random in translation 
parallel to the layer, and rotation about the 
normal. The powder pattern of such a material 
consists of two kinds of reflections—crystalline 
type reflections, and diffuse two-dimensional 
lattice reflections. The qualitative form of a two- 
dimensional lattice reflection has been given by 
von Laue! for the limiting case of large two- 


‘ dimensional gratings. It is the purpose of this 


paper to develop the general equations for x-ray 
powder patterns of random layer lattice materials. 


1M. v. Laue, Zeits. f. Krist. 82, 127 (1932). 


II. DIFFRACTION BY RANDOM LAYER 
LATTICE BLOCK 


The scheme of repetition in each layer is 
defined by the translation vectors a;d2. The 
separation between layers is given by the vector 
a3 perpendicular to a; and ay. The layers are 
identical and equidistant, but random in transla- 
tion parallel to the layer and rotation about the 
normal. It is enough to consider only the random 
translation since the only reflections allowed on 
the powder pattern are unaffected by an addi- 
tional random rotation. The position of atom 
in cell mymg, in layer ms is given by the vector 


m" = M10, + M22 (1) 


where 7, is the basis vector for atoms of type m, 
and 6 and ¢« have completely random values. 
Let us represent the directions of the primary 
and diffracted beams by unit vectors so and s. 
The intensity of diffracted radiation in electron 
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units is then given by 
T=>  f, Rm" ] fn expl — /A)(s —50)+ Rm” (2) 
Inserting (1) into (2) and summing, we have 
—So)+ Nia, (Ss —S0)- Node 


ms m3 


where Nia; and No2a2 are the dimensions of the layer, assumed to be a parallelogram, and F is the 


structure factor 
F=> f, exp[(2mi/A)(s — So) + tn]. (4) 


The two sine quotients will be small unless the two Laue equations are approximately satisfied. 
(5) 


If h and k are both zero Eq. (3) is independent of the randomness, and assumes the usual form for a 

crystalline reflection. 

(s—so)+ Nia, Node —So)+ 6) 


001 = 


If h and k are not both zero, the phase factors in the summation over ms; are completely random, the — 
scattering from the individual layers will be incoherent, and the intensity from a single layer will be 


Nia, 
(7) 


The diffracted radiation will be of two kinds; crystalline reflections of type (001), and two-dimen- 
sional lattice reflections of type (hk). There will be no reflections of type (hk/). We are interested in 
the powder pattern of a sample containing a large number of blocks with random orientation. The 
crystalline peaks of type (00/) are given directly by the usual powdered crystal theory. The next step 
is to develop the theory of the diffraction pattern of a two-dimensional lattice which takes with equal 
probability all orientations in space. 


III. PowpDER PATTERN OF Two- plane of the layer and 6; will be perpendicular to 
DIMENSIONAL LATTICE the layer. The solution of the two Laue equations 
Since the vector a3 is perpendicular to a; and = 8'V€S 
az, the reciprocal vectors and will be in the S—So=AH+(s—S 0) +33 (8) 
x where H=hb,+kbo. 


The significance of Eq. (8) is best seen in the 
reciprocal lattice construction shown in Fig. 1. 
AH The terminal point of the vector s—so can lie 
anywhere on the line KHK passing through the 
point H and perpendicular to the layer. For a 


two-dimensional layer the reciprocal lattice be- 
comes a series of parallel lines perpendicular to 
the layer. Since |s—so| =2 sin@, for each index 
Q pair hk there isa minimum value of sin@ given by 
Fic. 1. Reciprocal lattice construction for a , 
two-dimensional lattice. 2 sin@y=A/7. (9) 
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There will be a continuous distribution of 
diffracted intensity for all values of sin@ greater 
than sin@o. In the vicinity of 6) there will be a 
concentration in the diffracted intensity since 
|s—so| changes slowly for small values of ¢. 
Each (hk) band will have the form shown in 
Fig. 4. 

In the general case the layers will not be 
indefinitely large, and there will be an appreciable 
intensity contribution for the terminal point of 
s—So somewhat off the line KHK. To get the 
exact intensity distribution in a (hk) reflection, 
will therefore require the usual integrations in- 
volved in the derivation of an integrated inten- 
sity. Let P be the power, or total energy per 
second of diffracted radiation for a particular 
index pair (hk). 


p=f (10) 


P is obtained by integrating the intensity from 
each layer over the area of the receiving surface, 
and over dM the number of crystals in each 
range of orientation with respect to the primary 
beam. To be able to discuss the shape of a two- 
dimensional lattice reflection, it is necessary to 
know the distribution of power with respect to 
angle. Let us call this quantity Po», then P may 
be expressed by 


P= (11) 


If we represent any value of s—So in terms of a 
value (s—so)x which exactly satisfies the two 
Laue equations we have 


(12) 

From (7), the intensity per layer becomes 
Nya, Node 


By restricting our attention to one particular 
reflection (hk), the usual approximation can be 


made 
sin?Nx N? 
—N? esol 


sin’x 


Tne = exp{ — (#/A*)[(As- Nias)? 
(14) 


Fic. 2. Displacements of s and 5s» involved in the 
integration for diffracted power. 


We can express As in terms of the three variables 
Pipeps as follows 


1 2 


In Fig. 2 let the general vector (s—so) be the 
polar axis. The direction of the primary beam so 
varies by small azimuth and latitude angles 
doda and the direction of diffracted radiation by 
the orthogonal displacements d8dy. The solid 
angle covered by the primary beam is 


dQ=cos 6dade (17) 


and the area of receiving surface covered by the 
diffracted radiation is 


dA = R°dpdy. (18) 


If M is the number of two-dimensional layers in 
the sample, the number dM for which the 
primary beam lies within the solid angle dQ 
having the same equivalent orientation as regards 
a reflection (hk), is given by the proportion 


dM/mM =dQ/4r, 


where m is the two-dimensional multiplicity. 
Substituting in (10), we obtain 


ff ff (19) 


By holding ¢ constant, but allowing small 
changes da, dB, dy, the vector As will receive an 
increment d(As) given by 


d(As) = (da—dy) sindi+d8j 
+(da+dy) cosdk, (20) 


z | 
Y 
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where 7 j & are unit orthogonal vectors. If we 
replace the variables a, y by a new set 4, ¢€ at 
45° to the former axes, we have 


5=(a+y)/v2, e=(a—y)/v2, (21) 
d(As) = v2 sin@dei+dB8j+v2 cosédék. 


For increments ded8dé, the terminal point of As 
sweeps out a volume 


dV =sin26dedBd5. (22) 


From Eq. (15), for increments dp;dpedp3; we 
find the terminal point of As sweeps out a volume 


dV= (23) 
Equating the volumes, we have 
dedBdi = (24) 
Substitution in Eq. (19) gives 
f f f (25) 
Nov, 


The integration with respect to do can be carried 
out independently. Expressing J in terms of 
(16), we obtain 


_ ff f= 


—(/d*)(pi?+ po") Jdpidpedps. (26) 


By replacing the variables pips by a new set of 
orthogonal variables v,v2 such that v; determines 
the component of As parallel to the vector H we 
obtain 


+02") |dvidvedp3. (27) 
Let G be a unit vector parallel to H, and set 
v,=L(As-G). (28) 


The quantity Z introduced in (28) has the 
significance of a particle dimension given by 


1 
LINN 


|. (29) 


To get the distribution of power P2», the integra- 
tion is restricted to values of |s—so| between r 
and r+dr, where r=2sin@. In terms of the 
cylindrical coordinates used in Fig. 3, 


dv,dp3=Lrdgdr/ 
v,=L(As-G)=L(r 


Since the exponential factor involving v2 de- 
creases rapidly we can approximate by replacing 
the spherical shell by the tangent cylindrical 
shell and integrating v2 from — to +2. 


_ NiNsRim Md 


sind 
Xexp[ — \*)L?(r cose —AH)* gdr. (30) 


For the intense part of a reflection, only small 
values of yg are of importance, and cos¢g can be 
expanded. By letting 


aL 
(r—AH), (31) 


F(a) = f exp[ — (x?—a)? ]dx, (32) 


Eq. (30) becomes 


) 
sind 


F(a) coséd(26). (33) 


The power P29 is spread over a circle of circum- 
ference 27R sin26. Calling P29’ the power per 
unit length of circle, and introducing the in- 
tensity of scattering per electron, so that P29’ is 
expressed in absolute units, we have 


Px» 1+cos?20 
( ). (34) 


Py! =—— 
2rR sin26 m*c*R? 2 
TABLE I. Function F(a). 

a F(a) a F(a) a F(a) a F(a) 
—1.2 0.12 | —04 0.61) 0.1 0.96 | 0.6 1.06 
—1.0 0.21 | —0.3 0.69; 0.2 1.01 | 0.8 1.04 
—0.8 0.32 | —0.2 0.77/| 0.3 1.04 | 1.0 0.99 
—0.6 0.46 | —0.1 0.84} 0.4 1.06 | 2.0 0.67 
—0.5 0.53 0.0 0.91 | 0.5 1.07 | 3.0 0.52 


nr 
the 


(30) 
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be 


(31) 
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Fic. 3. Cylindrical coordinates used in the integration. 


Letting 
= (35) 
we obtain 
+cos*26) / L 
Px» =Km ) F(a), (36) 
2(siné)! \ 
where (37) 


The function F(a) as defined by Eq. (32) is 
given in Table I. 
For large values of a (a>3) 


F(a) =(x/4a)! 
and Eq. (36) reduces to 
KmF? (1+ cos*26) 


Poy! = . (38) 
4v2 sin@[sin@(siné —sin@) |! 


The usefulness of Eq. (38) is restricted by the 
expansion of cosg in Eq. (30). For sin@ <1.2 sin 
the error in Eq. (38) is less than 5 percent. A more 
useful equation for large values of a can be 
obtained from (27). From Fig. 3, we see that if r 
is appreciably larger than A//, the contribution 
to the integral for P29’ depends only on the 
length of the reciprocal lattice line KHK inter- 
cepted by the two spheres. The integration over 
and ve can be carried out — to +, 


and 
d =2 


Using the relation sin@)>=cos¢ sin@, and dividing 
by 2xR sin26, and multiplying by the intensity 


per electron, we find Eq. (26) becomes 


Km F?(1+cos*26) 
Po! = (39) 


4 


Equation (39) gives in absolute units, the result 
which has already been derived by von Laue.' 
To plot the distribution of intensity in a two- 
dimensional reflection Eq. (36) should be used 
for small values of sin@é—sin@o, and Eq. (39) for 
the larger values. Over small ranges of angle, 
F* can be treated as a constant, otherwise it 
must be considered as varying moennty with 
sin@, and given by 


2a32 
+" | (40) 


IV. PARTICLE SIZE FROM PEAK BREADTH 


The size of a two-dimensional layer in the 
plane of the layer can be determined from the 
breadth of the two-dimensional reflection at half- 
maximum intensity. The dimension LZ which is 
determined from a reflection (hk) is the effec- 
tive dimension of the layer in the direction 
H=hb,+kb2. Over the small range of angle in- 
cluded between half-maximum intensity values, 
the variation in (36) can be considered as due 
wholly to F(a). From a plot of Table I it is seen 
that the half-maximum intensity breadth corre- 
sponds to a change in a of amount Aa=3.26. 
From Eq. (37) 


27 
3.26= 


L 41) 


cos 6 


It should be pointed out that Eq. (41) for the 
two-dimensional particle size has exactly the 
same form as the equation for the crystalline 
particle size ‘except for a numerical coefficient 
about twice as large. 


V. DISPLACEMENT OF PEAKS 


The maximum in F(a) comes at a=0.55. 
Calling the position of the maximum @m, we 


have 
0.55 = — sind). 


dr. 
— 
| 
de- & 
ing 
ical 
«(| 
| 
j 
| 
(33) 
um- | 
per | 
is 
| .06 
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Fic. 4. Microphotometer record of the diffraction 
pattern of a heat treated carbon black, showing two- 
dimensional lattice reflections. Radiation Cu Ka=1.539A, 
monochromated by reflection from rocksalt. 


The maximum is therefore displaced from 6 
toward larger angles by an amount 


A(sin@) =0.16A/L. 


For small particle dimensions L, this displace- 
ment of the peak toward larger angles can be 
quite important. If the two-dimensional char- 
acter of the reflection were not realized, and the 
peak were treated as a normal crystalline powder 
pattern line, erroneous conclusions about lattice 
contractions could easily be drawn. 


VI. INTEGRATED INTENSITY OF A TWo- 
DIMENSIONAL REFLECTION 


In order to determine from a powder pattern 
the amount of material present in the form of 
independent two-dimensional layers, it is con- 
venient to measure the integrated intensity for 
a two-dimensional reflection. Although each re- 
flection extends to 2@=180° it is better to 
measure only the part extending to a convenient 
and arbitrary value of sin@é. We have 


20 
integrated intensity = f Pd(20). (43) 
0 


From Eq. (26) and Fig. 3, we see that if 20 is 
taken enough larger than 269, such that the 
intensity distribution at 26 is nearly independent 
of particle size, the contribution will depend only 
on the length of the reciprocal lattice line KHK 
within the sphere. Integrating p; and p2 from 
—« to + we find Eq. (26) becomes 


1 
f 
Va 


sin@ 
d(As) = |dp3= 2 sin@y 


By introducing the intensity per electron and 


dividing by 27R sin2@, to give the power per unit 
length of arc, and introducing the abbreviation 
K from Eq. (35), we obtain 


F?(1+cos?26) 
P’=Km f sec*ygdy. (45) 


2 sin@ sin2@ 


If the range 20—24 is kept small, an average 
value for F? and the polarization factor can be 
used. Since sin§y=sin@ cosy Eq. (44) becomes 


im 


2 sind, \ 


cos gdy 


and the integrated intensity is then given by 


0 2 
—sin?6,/sin20]! 
sin?@)/sin?6 ] ). (46) 


Equation (46) is in a convenient form for 
evaluation from microphotometer records, and 
for direct comparison with the integrated in- 
tensity of a crystalline reflection. 


VII. EXAMPLE OF Two-DIMENSIONAL 
LATTICE REFLECTION 


Figure 4 is the microphotometer record of part 
of an x-ray powder pattern of a heat treated 
carbon black. In terms of the hexagonal axes of 
graphite, the reflections appearing are the two- 
dimensional grating reflections (10) and (11), 
with the three-dimensional reflection (004) super- 
imposed on (10). The form of a two-dimensional 
grating reflection is shown very well by the 
(10) curve,? rising very rapidly near 6) and then 
decreasing continuously toward larger angles. 
From the breadth of the (10) peak at half- 
maximum intensity the layer dimension in this 
particular sample is calculated from (41) to be 
L=64A. 


2The two-dimensional lattice character of certain 
carbon black reflections has been “.— by a number 
of investigators. U. Hofmann ~ Wilm, Zeits. f. 
egg ong 42, 504 (1936); Arnfeld, Arkiv f. Mat. 
Astron. Fys. B23, 1 (1932; E. Berl, K. Andress, L. Rein- 
hardt ae W. Herbert, Zeits. f. physik. Chemie A158, 
273 (1932); A166, 81 (1933). 
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Non-Laue Maxima in the Diffraction of X-Rays from Rocksalt—-Equatorial Maxima 


G. E. M. Jauncey anp O. J. BALTzER 
Wayman Crow Hall of Physics, Washington University, St. Louis, Missouri 


(Received February 1, 1941) 


The existence of non-Laue maxima in the diffraction 
of x-rays was first recognized in 1938 by Wadlund for a 
continuous spectrum of x-rays and by Laval for monochro- 
matic x-rays. A brief description of Laval’s exhaustive work 
is given. The main points of the theories of Zachariasen, of 
Raman, Nilakantan and Nath, and of Preston, Bragg and 
Jauncey are also described. In the authors’ experiments 
CuKa, CuK, or MoKa x-rays with a weak background of 
continuous radiation were diffracted from an etched 
cleavage face of a thick rocksalt crystal. With the crystal 
face set to give the 200 reflection of CuKa; a total diver- 
gence of 1.6’ was found in the reflected beam, indicating 
very little warping in the layers near the crystal surface. 
With the crystal face turned by steps of 1° from a glancing 
angle of 45°07’ to 28°07’ the positions of the non-Laue 
maxima associated with the 400 and 620 Bragg reflections 
were found. The non-Laue spots are called associated 
Bragg spots. All photographs showed a shift of the maxi- 


mum intensity of the associated Bragg spot from the 
position of the corresponding Bragg spectrum line which 
is greater than zero. This result differs from that of Raman 
and Nilakantan, who report zero shift with rocksalt. The 
observed shift is always less than Zachariasen’s theoretical 
value for a given A, the angle by which the crystal is set 
away from the position for a Bragg reflection. Possibly for 
A<1° the Zachariasen formula holds. From the value of 4 
at which the associated Bragg spots disappear the shortest 
effective wave-lengths of the elastic waves concerned are 
about 13A to 17A. These are between 2 and 3 times the 
minimum wave-length in Debye’s specific heat theory. On 
the Preston-Bragg-Jauncey view a shift less than the 
Zachariasen value implies that the groups of atoms which 
give rise to the associate Bragg spots are considerably 
larger than an atom plus its 12 nearest neighbors as sup- 


posed by Preston. 


1. INTRODUCTION 


HE first realization that at least some of 

the diffuse spots and streaks which some- 
times appear on a Laue photograph and which 
do not belong to the Laue pattern may not be 
accidental in the sense of depending upon the 
particular crystal sample under test seems to 
have occurred in 1938. In that year Wadlund! 
using the continuous spectrum from a tungsten 
target tube found that a Laue photograph of 
rocksalt taken with a longer exposure than usual 
showed diffuse radial streaks in addition to the 
customary Laue spots. He obtained the same 
pattern of streaks with different samples of 
rocksalt. Similar streaks appeared with different 
samples of sylvine (KCI). In the same year 
Laval? began a systematic experimental study 
of the x-rays diffracted from the front face of a 
crystal when it is turned away from the Bragg 
position for reflection of monochromatic CuKa 
or MoKa x-rays. The results of further experi- 
mental study were reported during the first half 
of 1939 in the Comptes rendus* and a full and 
exhaustive report of the study of the effect with 


vA P. R. Wadlund, Phys. Rev. 53, 843 (1938). 
af: Laval, Comptes rendus 207, 169 (1938). 
Laval. Comptes rendus 208, 1512 (1939). 


monochromatic x-rays was published during the 
first half of 1939 in the Bulletin de la Société 
Frangaise de Minéralogie.* Lonsdale® has called 
attention to this excellent and exhaustive paper 
by Laval. We believe that the effect we are 
discussing was first established by Wadlund for 
continuous x-radiation and by Laval for mono- 
chromatic x-rays. 

Laval’s work* was done with MoKa and 
CuKea x-rays reflected from a bent crystal of 
quartz. The rays were then scattered by crystals 
of sylvine, rocksalt, diamond, calcite, and 
aluminum into an ionization chamber which 
could detect the arrival of one photon every two 
seconds. Most of Laval’s work was done with 
the non-Laue diffraction maxima for sylvine 
associated with the 200, 400, 600, 220 and 440 
reflections of monochromatic MoKa x-rays at 
room temperature (290°K). He also investigated 
the non-Laue maxima for sylvine associated with 
the 200 and 400 reflections at 550° and 665°. 
Laval describes his results in terms of the 
reciprocal lattice. Consider the vectors kp and k, 
each of whose magnitudes is 1/A and whose 
directions are, respectively, those of the incident 


fe Laval, Bull. Soc. Franc. Mineral. 62, 137 (1939). 
Lonsdale, Nature 146, 806 (1940). 
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and scattered beams of x-rays. Laval calls the 
vector k—Kkp (see Zachariasen’s paper®) the 
vecteur de diffusion. Drawing this vector from 
the origin of the reciprocal lattice, Laval calls 
the point at the other end of the vector the 
pole de diffusion. We shall translate this by 
“diffraction pole.’’ The magnitude of k—kp is 
(2/A) sin}¢, where ¢ is the angle of scattering. 
When the diffraction pole falls on the Akl 
reciprocal lattice point, the hki Bragg reflection 
occurs. According to Laval, each reciprocal 
lattice point is surrounded by a region of strong 
scattering. If the diffraction pole falls within 
this region, the intensity of scattering is large 
but of smaller order of magnitude than that of 
Bragg reflection. Between the regions of strong 
scattering surrounding each reciprocal lattice 
point there exists the region of weak scattering. 
If the diffraction pole falls within this region, 
the intensity of scattering is of the order of 
magnitude of Jauncey, Harvey and Woo’s 
diffuse scattering.’ Laval states that each region 
of strong scattering is divided into subregions 
by equi-scattering (isodiffusion) surfaces. If the 
diffraction pole is moved along an equi-scattering 
surface, the intensity of the scattered rays 
remains constant. Laval experimentally deter- 
mines these surfaces to be surfaces of revolution 
about the line joining O, the origin of the 
reciprocal lattice, to M, the Akl point of the 
reciprocal lattice. An equi-scattering surface is 
not spherical. Its dimension is greater in the 
direction perpendicular to OM than its dimension 
along OM. Further, the distance from M to an 
equi-scattering surface in the direction from 
to O is less than that in the opposite direction. 
The maximum intensity of a non-Laue diffraction 
spot will be determined by the position of the 
point of tangency between the sphere of reflection 
and an equi-scattering surface surrounding the 
hkl reciprocal lattice point for a given A, the 
angle by which the /k/ planes® are turned from 
the Bragg angle @ for the hk/ reflection. As a 
result of the shape of Laval’s equi-scattering 

®W. H. Zachariasen, Phys. Rev. 57, 66A (1940); ibid. 
57, 597 (1940). 

7G. E. M. Jauncey and G. G. Harvey, Phys. Rev. 37, 
1203 (1931); Y. H. Woo, Phys. Rev. 38, 1 (1931) and 41, 21 
(1932); G. E. M. Jauncey, Phys. Rev. 42, 453 (1932). 

8For brevity we shall speak of the hk/ planes even 


though we are dealing with, say, the 622 reflection which is 
a second-order reflection from the 311 planes. 
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surfaces, the shift ¢.2—26 is less than that 
given by Zachariasen’s theoretical formula® 


2 sin@ cosj 
tangs = (1) 
1—2 siné@ sinj 


or Zachariasen’s approximate theoretical formula 
das — 20=2A sin*?, (2) 


where ¢cz is the angle of scattering for the rays 
arriving at the position of maximum intensity 
of the non-Laue diffraction maximum associated 
with the hkl Bragg reflection, j is the angle of 
incidence on the hki planes, and A=j—8. 
Equation (1) was derived independently by 
Raman, Nilakantan and Nath.*® Both deriva- 
tions, however, depend on the assumption of 
spherical equi-scattering surfaces surrounding 
the Aki reciprocal lattice point. 

The experimental work of Laval with mono- 
chromatic rays has been confirmed by Preston,'® 
Siegel and Zachariasen,'"' Raman and Nila- 
kantan,"-'" Lonsdale, Knaggs and Smith,” 
Jauncey and Baltzer,'* and Siegel.'7 Raman and 
Nilakantan’s results have been discussed by 
Zachariasen.'*: 


2. DIscUSSION OF THEORY 


Zachariasen* and Raman and Nath’ give the 
vector relation 


(3) 


where the vectors b; (i=1, 2, 3) are reciprocal to 
a; (i=1, 2, 3), the edges of the unit cell of the 
crystal. The magnitude of < is 1/A, where A is 
the wave-length of elastic waves of frequency 


°C. V. Raman and P. Nilakantan, Proc. Ind. Acad. Sci. 
11A, 379 (1940); C. V. Raman and N. Nath, Proc. Ind. 
Acad. Sci. 12A, 83 (1940). 

0G, D. Preston, Proc. Roy. Soc. A172, 116 (1939). 

1S. Siegel and W. H. Zachariasen, Phys. Rev. 57, 66A 
(1940) ; 57, 795 (1940). 

12 C. V. Raman and P. Nilakantan, Proc. Ind. Acad. Sci. 
11A, 389 (1940); Curr. Sci. 9, 165 (1940). ; 

13C, V. Raman and P. Nilakantan, Proc. Ind. Acad. Sci. 
11A, 398 (1940). ; 

4 C, V. Raman and P. Nilakantan, Proc. Ind. Acad. Sci. 
12A, 141 (1940). 
(1940). Lonsdale, I. Knaggs and H. Smith, Nature 146, 332 

16 G. E. M. Jauncey and O. J. Baltzer, Phys. Rev. 58, 
1116 (1940). 

17S. Siegel, Phys. Rev. 59, 371 (1941). 

18 W. H. Zachariasen, Nature 145, 1019 (1940). 

19 W. H. Zachariasen, Phys. Rev. 59, 207 (1941). 
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in the crystal. Raman and Nath’ write (3) in 
the form 
(n/d)+(1/A)=n/dn, (4) 


where 1/d, 1/A, and 1/d,, are treated as vectors 
and d is the grating space for the hkl planes, and 
where m is the order of reflection for the hkl 
planes. The effect of the elastic waves of the 
crystal is to modify (‘“‘modulate”’ according to 
Laval) the grating space d so that it assumes 
the value d,, for a small fraction of the time 
that the x-rays are incident on the crystal. 
According to Zachariasen the elastic waves are 
a consequence of the thermal agitation of the 
atoms of the crystal as in Debye’s theory of 
specific heats of solids.*° According to Raman 
and Nath the elastic waves are excited by the 
incident x-rays by means of a quantum process. 
Zachariasen’s theory gives + parallel or anti- 
parallel to k. Raman and Nath make two 
tentative assumptions concerning +. Their first 
assumption corresponds to the result of Zachar- 
iasen’s theory just mentioned and so they arrive 
at Eqs. (1) and (2). In their second assumption 
« is taken perpendicular to the line joining the 
origin and the /k/ point of the reciprocal lattice 
and in the plane of incidence for the hk/ planes. 
This second assumption leads to 


sin}dap COS(j — = sind (5) 
and the approximate relation 
—20=0, (6) 


even though j7#@. Laval’s experimental results 
give the shift (¢as—26) as being between the 
value given by (1) or (2) and that given by (5) 
or (6), with perhaps the experimental shift 
closer to (6) than to (2). 

Discussing Wadlund’s results with rocksalt 
and sylvine, Zachariasen*! suggested that the 
diffuse streaks might be due to a two-dimen- 
sional lattice appearing on the surfaces of the 
mosaic blocks which constitute these erystals. 
Preston’® concludes that his results with alumi- 
num do not support this idea. Instead, Preston 
suggests that the crystal is broken up into 
groups of atoms by the thermal vibrations of the 
lattice. The groups form three- but not two- 


20P. Debye, Ann. d. Physik 39, 789 (1912). 
21W.H. Zachariasen, Phys. Rev. 53, 844 (1938). 


dimensional lattices. Bragg’ has examined 
Preston’s idea for the case of scattering by a 
group of eight point-atoms arranged at the 
corners of a cube of edge a. Jauncey™ has further 
examined the scattering of a simple cubic array 
of point-atoms with N atoms on an edge and 
with a distance a between adjacent atoms. 
Jauncey has found that the theory gives rise to 
(2) for small values of A=(j— 86) and for values 
of N=2, 3, 4, 5. This is an astonishing result. 

Equation (2) fairly well fits the results ob- 
tained for rocksalt by Siegel and Zachariasen," 
for diamond by Raman and Nilakantan,” and 
for sylvine by Lonsdale, Knaggs, and Smith.” 
Since Eq. (2) can be obtained by at least two— 
if not three—different theories, experimental 
agreement with (2) does not support any one of 
these theories above any other. It will be 
necessary to study the intensity as well as the 
law of position in order to choose between the 
theories. However, the experimental law of 
position is in some doubt. In addition to Laval’s 
finding of a shift (¢42—20) less than given by 
(2), Raman and Nilakantan report results with 
sodium nitrate and rocksalt'* which fit (6) 
rather than (2). We have therefore undertaken 
an experimental study of the law of position of 
non-Laue maxima for rocksalt. 


3. EXPERIMENTAL METHOD 


Instead of the method of penetration through 
a crystal plate used by Raman and Nilakantan,™ 
we have, following Laval?‘ and Siegel and 
Zachariasen,"' used the method of scattering 
from the 100 cleavage face of a thick rocksalt 
crystal. Raman and Nilakantan using the 
penetration method for rocksalt were forced to 
use MoKa and MoK§ rays in order to obtain 
sufficient intensity. The optimum thickness in 
the penetration method equals 1/, for scattering 
in the forward direction, where yu is the linear 
absorption coefficient. For MoKa rays in rock- 
salt, «= 18.0 cm~ and the optimum thickness is 
0.55 mm, while, for CuKa rays in rocksalt, 
u=165 cm™ and the optimum thickness is 0.06 
mm. It is impractical to use CuKa rays in the 
penetration method because of the great loss of 
2 W. H. Bragg, Nature 146, 509 (1940). 


3G. E. M. Jauncey, Phys. Rev. 59, 456 (1941); Nature 
147, 146 (1941). 
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FILM 


Fic. 1. Diagram 
of apparatus. 


intensity if a reasonably thick plate is used. On 
the other hand, when MoKa rays are used, the 
shift ¢22— 20 becomes small if Eq. (2) represents 
the experimental law. For instance, sin@ for Beso 
(the 620 Bragg reflection) is 0.4 and ¢as—26 
= 1.6° for 7—6@=5°. But if CuKa rays are used, 
siné=0.863 and ¢as—20=7.45°. If it is desired 
to show that ¢.2—2600, this can be shown more 
convincingly with CuKea rays than with MoKa 
rays but in the case of CuKa rays the front 
face method should be used. 

The experimental arrangement is shown in 
Fig. 1. The primary beam of x-rays makes a 
glancing angle of incidence 7 on a 100 cleavage 
face of rocksalt. Scattering occurs in various 
directions ¢. The circle represents the photo- 
graphic film which is bent into a cylinder about 
a vertical axis through C. The radius of the 
cylinder is 9.64 cm. Each slit, is rectangular, 
being 0.275 mm wide and 1.45 mm high. The 
slits are 10 cm apart. From the values of 7 at 
which various Laue spots disappear, Amin Was 
estimated to be 0.58A when the copper target 
tube was used so that the tube was operated at 
21 kv peak. A molybdenum target tube operated 
at about 35 kv peak was used for a few pictures. 

All crystal surfaces have been etched with a 
twenty-percent solution of alcohol in distilled 
water in order to remove any warped surface 
layers which might be present. The 200 Bragg 
reflection of CuKa; gave a line 0.7 mm wide on 
a film placed 150 cm from the axis about which 
the crystal may be turned. The width of 1.6’ 
calculated from this data compares well with 
that corresponding to the half-width of 87” at 
half-maximum found by Kirkpatrick and Ross* 
for rocksalt with a double crystal spectrometer. 
The surfaces of our crystals were thus reasonably 
free of warping. 

It may -be objected that CuKa x-rays are so 


*P. Kirkpatrick and P. A. Ross, Phys. Rev. 43, 596 
(1933). 
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highly absorbed that in the front face method 
the rays only penetrate a few atomic layers and 
that therefore in this method a surface phenom- 
enon is observed. This might explain such 
difference as occurs between any results we may 
obtain and those of Raman and Nilakantan. 
In the arrangement of Fig. 1 the intensity which 
is scattered from a layer of the crystal next to 
the surface is proportional to 


1—exp[ (7) 


where ¢ is the thickness of the layer. The thick- 
ness which gives rise to one-half the intensity 
in an aB (associated Bragg) spot is therefore 


0.693 /u{ (1/sinz)+1/sin(@—1)}. (8) 


An spot occurs at 1=44°40’ and 
=124°06’ with CuKa rays. The _half-value 
thickness is therefore 1.6X10-* cm. A linear 
dimension of the mosaic blocks of which rocksalt 
is believed to be composed is of the order of 
10-* cm. Hence the layer which contributes one- 
half of the intensity of the aBg2o spot is about 
16 mosaic blocks thick. Had the thickness been 
a fraction of that of a mosaic block, we would 
have abandoned this method but, as it is, we 
feel some justification in assuming that the layer 
is thick enough for the crystal surface effect to 
be relatively small. 

In this paper we have limited ourselves to a 
study of the equatorial (4k0) associated Bragg 
spots. In each photograph a Bragg line is placed. 
on the film by rocking the crystal through a 
small angle about 7=7,, where 7, is the glancing 
angle of incidence on the 100 cleavage face 
which gives rise to an hkO Bragg reflection. The 
exposure is for an interval of a few seconds up 
to a minute depending on the particular Bragg 
line. For the higher orders of reflection the Bragg 
line is a doublet due to the presence of CuKa; 
and CuKae. The crystal is then turned away 
from the Bragg position to another value of 4. 
With this new value of 7 the film is exposed for 
some 10 to 20 hours. For equatorial spots, the 
Lixo spot (hkO Laue spot) is displaced from the 
position of the Bixo line (ARO Bragg reflection) 
by 2(¢—ig). All angular measurements are made 
from the original films or from microphotometer 
traces of these films. For equatorial spots 


(t—tz) =(j— 4). 


| | 
| 
| 
| 


NON-LAUE MAXIMA IN 


4. EXPERIMENTAL RESULTS 


A sample photograph is shown in Fig. 2. 
The various Laue spots are indicated on the 
figure. In addition the Byoo, Beoo, and Begzo lines 
for CuKa; and CuKag are shown and these are 
used as fiducial marks for standardizing the 
photograph. We note the three strong associated 
Bragg spots, @Byooa, and aBgyo3, and a 
faint diffuse spot aBg2o3. The subscripts a and B 
refer to the CuKa and CuK@ spectrum lines, 
respectively. The angle of glancing incidence on 
the 100 cleavage face for this photograph was 
i=30° 50’. In Fig. 3 is shown the reciprocal 
lattice construction for 7=30° 50’. Two circles 
a and 6 with centers, respectively, at P, and Pg 
are drawn. These represent the spheres of 
reflection for the CuKa and spectrum 
lines, respectively. It is to be noted that the a 
circle passes close to the 400 point. Therefore 
an @Byo¢ spot should appear on the correspond- 
ing photograph. As seen in Fig. 2 this spot does 
appear. Further, the 6 circle of Fig. 3 passes 
close to the 400, 620 and 640 points. In Fig. 2, 
the corresponding @B4o93 and aBg4o3 spots appear. 
If Fig. 3 is carefully drawn, the point 620, 
though close to the 8 circle, is further from the 
circle than either of the points 400 or 640. It is 
a rule that the further a reciprocal lattice point 
is from the sphere of reflection the fainter, but 
more diffuse, the corresponding aB spot is. We 
should then expect the aBg2o3 spot to be faint 
and as seen from Fig. 2 it is faint (it may not 
show on the reproduction). 

In Fig. 4 are shown strips cut from 18 photo- 
graphs similar to Fig. 2. The first or top strip 
gives the equatorial spots for i=45° 07’. The 7 
for successive strips decreases by 1° until the 
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a 4008 \ B400a 
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Fic, 2. Sample of photograph showing the associated 
Bragg reflections of CuKa and CuK x-rays for a glancing 
angle of incidence of i=30° 50’ on the 100 cleavage face of 
rocksalt. The figure is reduced to about one-third of the 
size of the original photograph. 
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640 620 


Fic. 3. Reciprocal 
lattice construction for 
the experimental con- 400 
ditions of Fig. 2. Geo- 
metric construction 
shows that the diffuse 
spots aB400a, aB4008, 
aB6208, and aB6408 
are likely to occur for 
¢= 30° 50’. 


oho 


18th or bottom strip is reached when i= 28° 07’. 
The Bragg reflection lines for CuKa rays fall in 
the vertical lines headed by Byooa, Beooas Bezoa in 
Fig. 4. These lines are used as fiducial marks 
for the respective strips. The 400 and 620 Laue 
spots fall into the slant lines headed by Lago and 
Lez in Fig. 4. It is seen that as the Laoo spot 
approaches the Byoog line, an associated Bragg 
spot (@Byooq) first appears in the 6th strip (this 
is so on the original negative but may not be 
evident in the reproduction). The aB spot is 
always between the B line and the L spot. 
Experimentally the aB,,.. spot is associated with 
the B,,; line. If an Lyn spot is made to ap- 
proach a B,,; line, no aB spot appears unless 
h’k'l’ =hkl. When the L4oo spot falls on the Byoo 
line, photographic reversal and halation appears 
as shown in the 13th strip. Following the progress 
of the Les spot down from the first strip, we 
note that the aBgoq spot disappears in the 11th 
strip (on the original negative). A slight indica- 
tion of photographic reversal is shown for Le2o 
in the 5th strip. The center of the Laue spot is 
not quite in the position of the Bgzoa line but the 
Laue spot has a width due to divergence of the 
incident beam as determined by the slit system. 
The divergence is sufficient to cause some Bragg 
reflection. 

The experimental values of the shift (¢an— 26) 
are given in Table I. The value of igso0a= 94000 
is 33° 07’ and that for igg20e is 41° 15’. 

From Table I we conclude that the shift 
(¢ap— 26) is less than that given by Zachariasen’s 
Eq. (2). The shift, however, is not zero as 
reported by Raman and Nilakantan." Possibly 
at small values of A (1° or less), the shift is 
given by (2). Our result is in general agreement 
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B600a B620a 1.620 


Fic. 4. Equatorial reflections of copper x-rays from 
cleavage face of rocksalt. Glancing angle of incidence is 
4#=45° 07’ for top strip with 7 decreasing by 1° for each 
successive strip. The lines B400a, B600a, and B620a are 
placed as fiducial marks on the strips. An associated Bragg 
spot @B620a is visible on the upper strips. The associated 
Bragg spot aB400qa is first apparent (on the original 
negative) for the 6th strip, becomes more intense and less 
diffuse approaching the 13th strip, and then gradually 
decreases in intensity. A diffuse reflection aB4008 appears 
on the last four strips. 


with Laval’s results. The equi-scattering surfaces 
surrounding the 400 and 620 reciprocal lattice 
points are not spherical but are shaped in a way 
similar to that found by Laval. 

We have obtained a few films with a molybde- 
num target tube. The microphotometer trace 
for the aBgooa spot obtained with A=3° 00’ 
shows a shift of 0° 22’+0° 11’. This is greater 
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than zero but is less than 0° 51’ as given by (2). 
From (8) and ».=18.0 cm~ the half-value layer 
for the aBgooa spot is 6.5X10-* cm or about 65 
mosaic blocks thick when MoKa rays are used. 
Since our result with MoKa rays qualitatively 
supports our results with CuKa rays, we feel 
that our conclusion is not a consequence of the 
effective layer of rocksalt being too thin. 


5. DIscussION OF RESULTS 


From Fig. 1 of Zachariasen’s paper® it can be 
shown that the magnitude of ¢ which gives rise 
to strong scattering in the scattering direction 
¢ is 
=1/A=(2/d) 

—2sin@ sin}¢ cos(}¢—j)}*. (9) 


This reduces to the approximate relation 
A=)/{8+4A(A—8) sin20}}, (10) 


where 6=@— 20. If we wish to find the A which 
gives rise to the maximum intensity of an 
associated Bragg spot, we replace ¢ by ¢az in the 
formula for 6. The values of 6=(¢ag—26) and 
A for the least visible aBiooa spot are 2.4° and 
7.0°, respectively, so that from (10) A=13.3A. 
The respective values for the least visible aBg2oa 
spot are —5.4° and —5.13°, so that A=17.5A. 
In the Debye” theory of specific heats, the 
minimum value of A for a simple cubic crystal 
is twice the grating space, so that, treating rock- 
salt as a simple cubic crystal at least for the 
equatorial associated Bragg spots, we find 
Amin =5.628A. As Fine*> points out the solid is 
treated as a continuum in the Debye theory until 
Amin is approached. At this place in the theory it 
is suddenly remembered that the solid consists of 
atoms separated by pretty definite distances. In 
the case of tungsten crystals Fine finds that the 
vibration spectrum, instead of having a maxi- 
mum at a frequency equal to Debye’s ymax, has 
two sharp maxima at frequencies 0.82 and 0.55 
times the Debye ymax. Though rocksalt is a 
polar and tungsten a nonpolar crystal, it may 
well be that the A spectrum exhibits a maximum 
in the neighborhood of A=12A. With such a 
maximum one could not expect a 400 reciprocal 
lattice point of rocksalt or sylvine to be sur- 


% P. C. Fine, Phys. Rev. 56, 355 (1939). 
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rounded by a spherical equi-scattering surface 
which would be tangent to the sphere of reflection 
when A=7°. Zachariasen’s theory as at present 
stated must break down for rocksalt at A=2 or 
3 times 5.628A. Furthermore, Brillouin®* points 
out that the instantaneous distribution of density 
in an elastic wave can no longer be treated as 
sinusoidal when Debye’s Amin is approached. 
Reflection of x-rays from elastic waves is limited 
to first-order reflection only when the waves 
may be treated as sinusoidal. If A approaches 
Amin, the reflection takes place according to 


nd = 2A siné. (11) 


One may expect that in (4) the term 1/A should 
be multiplied by an integer m’ when A approaches 
Amin- 

From Table I it may be concluded that the 
Zachariasen formula for the shift holds for the 
aByooa spot when A=1°. The value of A for this 
case is 96A, which is 17 times Amin. In the region 
0°=A=1° the corresponding elastic waves may 
be treated as sinusoidal and in this region 
Zachariasen’s theory may hold. However, there 
is the question of the alignment of the mosaic 
blocks. Compton and Allison*’ in discussing the 
work of others with the double crystal spec- 
trometer quote half-widths at half-maximum of 
7 seconds, 87 seconds, and 30 minutes with 
different samples of rocksalt. On account of 
these facts there is great uncertainty in the value 
of A when A=1°. For this reason Laval com- 
pletely avoided the region of A=1°. We believe 
that the experimental field separates into three 
regions to be studied—the region of A=1°, the 
region of strong scattering for A>1°, and the 
region of weak scattering. It would be very 
interesting to study the effect of temperature in 

Ann. de 17, 88 (1922). 


27 A. H. Compton and S. K Allison, X-Rays in Theory 
and Experiment (Van Nostrand, 1935), pp. 401, 729. 


the first region. It is well established that the 
effect of a rise in temperature is to increase the 
intensity of the scattering in the second and 
third regions. 

So far in this section the discussion has been 
in terms of Zachariasen’s theory. However, there 
is the Preston-Bragg-Jauncey theoretical view- 
point. In Jauncey’s paper™ curves are shown for 
the scattering from a simple cubic lattice with 
N atoms on an edge. For A=5° and the 620 
associated Bragg reflection of CuKa x-rays from 
rocksalt (treated as a simple cubic crystal), the 
curves show a maximum at (¢ag—26) given by 
Zachariasen’s Eq. (2) for N=2, 3, 4. There is a 


definite breakdown at N=6 and the shift 


TABLE I. Shift (¢c8—20): CuKa x-rays. 


aByoo aBe20 
—20 ¢aB —20 
A Exp. Eg. (2) A Exp. Eg. (2) 
7.0° 2.4° 4.17° 3.87° 4.5° 5.76° 
6.0 2.4 3.58 2.87 3.2 4.27 
5.0 2.1 2.98 1.87 1.8 2.78 
4.0 1.6 2.38 0.87 0.9 1.30 
3.0 1.3 1.79 —0.13 0.0 —0.19 
2.0 0.8 1.19 —1.13 | —1.2 — 1.68 
1.0 0.6 0.60 —2.13 | —2.3 —3.17 
0.0 0.0 0.00 —3.13 | —3.6 —4.66 
-—1.0 —0.4 —0.60 —4.13 | —4.5 —6.15 
—2.0 —0.6 —1.19 —5.13 | —5.4 —7.64 
—3.0 —1.2 —1.79 
—4.0 —1.8 —2.38 
—5.0 —1.8 —2.98 


becomes less than Zachariasen’s shift. From 
our experimental results we would in terms of 
this theory conclude that N>6. It is hoped to 
discuss this matter further in a later paper. 
Note added in proof.—In references 16 and 23 
the term “modified Bragg reflection” was used. 
In the present article this term has been replaced 
by the term “associated Bragg reflection’ and 
the symbol mB by the symbol aB. The reflection 
indices in the references should be doubied. 
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Evidence for the Radioactivity of Slow Mesotrons 
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A beam of mesotrons is selected by means of a fourfold coincidence system of counters, and 
is allowed to pass through a block of iron 10 cm thick. The absorption of mesotrons by the 
iron is recorded by means of a battery of anticoincidence counters. It is found that a certain 
fraction of the stopped mesotrons is associated with the emission of an ionizing particle from 
the absorber. A special coincidence recording system, whose resolving time is of the order of 
one microsecond, enables one to establish that the emission of the particles is delayed with 
respect to the passage of the mesotron by a few microseconds. The delayed particles are 
interpreted as the electrons resulting from the 8-decay of the mesotron. The present experiment 
establishes only the order of magnitude of the mean life of the mesotron at rest, but more 


quantitative measurements are in progress. 


ECENT experiments on the absorption of 

mesotrons! by the earth’s atmosphere and 

by condensed materials have brought convincing 

evidence for the hypothesis of the decay of the 
mesotron. 

Much less satisfactory is the present situation 
concerning the detection of the disintegration 
electrons. Williams and Roberts? obtained two 
cloud-chamber photographs which seem to de- 
monstrate the process of 6-decay of a mesotron ; 
on the other hand, other assumed mesotron 
tracks have been seen to end in the chamber 
without producing any ionizing particle. The 
evidence from the equilibrium between the elec- 
tron and mesotron components in the atmosphere 
is, at best, inconclusive.* 

The detection of delayed coincidences from the 
disintegration electrons produced by stopping 
mesotrons in a block of lead was attempted by 
Montgomery, Ramsey, Cowie and Montgomery, 
with a negative result.‘ 

I wish to report some results which seem to 
bring positive evidence for the delayed emission 
of the disintegration electrons and to afford a 
rough measurement of the mean life of the 
mesotron at rest. 

The counter set-up is illustrated in Fig. 1 
(circles represent the effective cross sections of 


1 See B. Rossi and D. B. Hall, Phys. Rev. 59, 223 (1941), 


also for earlier literature. 
oui J. Williams and G. E. Roberts, Nature 145, 102 


(1940). 

3G. Bernardini, B. N. Cacciapuoti, B. Ferretti, O. 
Piccioni and G. C. Wick, Phys. Rev. 58, 1017 (1940). 

4C. G. Montgomery, W. E. Ramsey, D. H. Cowie and 
D. D. Montgomery, Phys. Rev. 56, 635 (1939). 


the counters; counters designated by the same 
letter are connected in parallel; cross-marked 
counters are used as anticounters). 

Counters A, B, C and D define a beam of 
mesotrons, which impinges upon the battery of 
anticounters F. The number of anticoincidences 
(ABCD-F) is about 1.5 percent of the number 
of coincidences (ABCD), when no absorber is 
present between D and F. A block of iron 
10X2.5 X40 cm’, placed as indicated by Fig. 1, 
brings the number of anticoincidences to about 
5 percent, in agreement with the known absorp- 
tion coefficient of mesotrons at sea level. 

A certain fraction of the anticoincidences 
(ABCD-—F) is associated with a discharge of one 
of the counters E. This may be due partly to 
mesotrons scattered by the iron, partly to dis- 
integration electrons from mesotrons stopped by 
the absorber. Convincing evidence in favor of 
the latter process will be obtained if it can be 
proved that the passage of a particle through E 
is delayed with respect to the discharge of 
counters (ABCD). 

For this purpose, the counters were connected 
to a system of amplifier units, indicated in 


Fig. 1 as AU1, AU2, AU3 and h.r.p. (for high ~ 


resolving power) unit A4. Unit 1 is a fivefold 
coincidence set (resolving time 3.6X10- sec.) 
connected to counters (ABCDE). The output 
pulse is passed onto unit 2, and operates recorder 
R2 only if it is associated with an anticoincidence 
from counters F and G. The output pulse from 
unit 2 is passed onto unit 3, which is a double 
coincidence set, and operates recorder R3 only 
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when a simultaneous (within 10~* sec.) pulse is 
fed to it from the h.r.p. unit 4. The latter is a 
double coincidence set, connected to counters 
(DE), and has a resolving time of 1.2 X10~° sec. 

If the discharge in the E counters is not de- 
layed, recorder R3 will record the same number 
of pulses as R2, since every process (A BCDE-FG) 
is also a double coincidence (DE). But, if the 
discharge in E is delayed by more than 1.2K 10-° 
sec. and less than 3.6X10- sec. with respect to 
the fourfold coincidence (ABCD), then it will be 
recorded by R2 and not by R3. Thus, the 
difference in the number of counts R2—R3 will 
give the number of particles emitted from the 
absorber within the above time limits after the 
stopping of a mesotron. 

Two obvious causes of error must be avoided. 
First, one must be certain that practically all 
true systematic coincidences (DE) are recorded 
by the h.r.p. set. More than twenty separate 
checks, distributed through the series of experi- 
ments, showed that the efficiency of the h.r.p. 
set for true coincidences was always better than 
99.5 percent. 

Another source of error might reside in the 
random coincidences between systematic four- 
fold coincidences (A BCD) and pulses in E, which 
fall within 3.6 X 10~ sec. but not within 1.2 x 10-* 
sec. Their number ought to be only about one in 
fifteen days. However, it was considered safer to 
take readings alternately with and without ab- 
sorber. In this way, all random coincidences and 
other possible spurious effects would show in the 
blank runs. 

An easily measurable number of delayed 
pulses in E was observed with the iron absorber 
in place, as contrasted with an almost negligible 
number in the blank runs. Table I summarizes 
the results. 

Since the number of mesotrons absorbed by 
the iron per hour is approximately known, the 
observed number of delayed coincidences is 
found to correspond to about seven percent of 
the number of absorbed mesotrons. Considering 
geometrical factors, absorption of the electrons 
by the iron, and the fact that probably about 40 
percent of the mesotrons decay within 1.2 micro- 
seconds, the observed percentage of electrons 
does not appear to be inconsistent with the 
assumption of one emitted electron per mesotron. 


The results would afford a measurement of the 
mean life + if one knew that all additional co- 
incidences produced by the absorber, falling 
within 1.2 microseconds, are due to disintegration 
electrons. Actually, some of them will be due to 
scattered mesotrons, and the value of 7 will 
represent a lower limit. From the data of Table I, 
one finds: 


exp (—1.2X10-*/r) >42/111; 7>1.2X10~* sec. 


Iron was used instead of lead because of the 
more favorable stopping/scattering ratio. Some 
measurements were also taken with an Al ab- 
sorber and gave similar results. Besides the 
strong scattering of slow mesotrons, the absorp- 


O © RI 
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Fic. 1. Arrangement of counters and amplifier units. 
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TABLE I. Experimental results. 


NUMBER OF PULSES PULSES 


ABSORBER Hours R2 R3 R2—R3 
None 231.6 96 93 342 
10 cm Fe 231.6 207 162 45147 


tion of disintegration electrons due to radiative 
losses is another reason against the use of a heavy 
element, and may possibly explain the failure of 


Montgomery ef al. to observe the disintegration 
electrons. 

The experiments are being coritinued with an 
improved arrangement that is expected to give a 
reliable value of the mean life.® 


5 Note added in proof.—Preliminary results with the 
new arrangement have been announced in a letter to 
The Physical Review (Phys. Rev. 59, 613 (1941)). Experi- 
ments in progress, however, seem to indicate a mean life 
shorter than already reported, of the order of 2 micro- 
seconds or possibly less. 
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Without specifying otherwise the stress-energy distribution inside of a material particle, the 
assumption is made that the scalar quantity T, known in electromagnetism as the ‘‘scalar of 
Laue,” shall vanish, in harmony with the Maxwellian stress tensor, and confirmed also by other 
considerations. It is shown that under this condition a theoretical explanation can be given of 
the long established empirical fact that the mass of any particle is necessarily positive. All 
metrical contributions to the total mass come out as necessarily positive, second-order quan- 
tities, with the only exception of an eventual spin, which has a decreasing instead of increasing 
effect on the mass. The strict proportionality of gravitational and inertial mass can still be 
established, at least for the spherically symmetric case. The factor of proportionality, however, 
is numerically different from the value ordinarily assumed. 


1. INTRODUCTION 


NE of the peculiar mysteries of nature is 

the fact that all masses are positive. 
Electric charges can take positive and negative 
values. The masses of elementary particles show 
a great diversity in values, but the sign is always 
positive. The Newtonian force of universal 
gravity is always a force of attraction and not a 
force of repulsion. Why is it that we do not 
encounter negative masses in nature? The 
present paper attempts to give the theoretical 
explanation of that long established empirical 
fact, on the basis of general relativity. 


2. Basic ASSUMPTIONS 


In the Newtonian theory of gravity the mass 
m of a particle appears as a constant of inte- 
gration. The potential equation is 


Ag=0. (2.1) 


The spherically symmetric solution of that 
equation is 


o=m/r. (2.2) 


The constant m is in physical interpretation the 
mass of the attracting particle. 

In general relativity the potential function ¢ 
is replaced by the much more comprehensive 
scheme of ten giz quantities, defined by the line 
element 


ds? = (2.3) 


The field equation (2.1) is replaced by the 
Einsteinian field equations: 


Riu =0, (2.4) 


where Rx is the contracted Riemann-Christoffel 
tensor. The spherically symmetric solution of the 
field equations (2.4) is the line element of 
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Schwarzschild : 


ds?= 


1—m/r 


m 
+(1 (2.5) 
r 


This is the solution that corresponds to the 
solution (2.2) of the Newtonian theory. Here 
too the constant m is merely a constant of 
integration which a priori may assume any 
positive or negative values. Any specification of 
the constant m is just as impossible on the basis 
of the homogeneous field equations (2.4) as it is 
on the basis of the homogeneous field equation 
(2.1). An explanation of the necessarily positive 
character of m requires a generalization of Eqs. 
(2.4), together with the requirement that only 
regular solutions of the new field equations shall 
be permitted. 
Even so the tensor 


Tix=Riu— (2.6) 


will remain the dominating quantity of the field. 
It represents in physical interpretation the 
matter tensor. It is feasible to assume that the 
field equations (2.4) are still good outside of a 
certain action radius p of the particle. The 
spherically symmetric solution (2.5) of the field 
equations (2.4) will still have its significance, 
although restricted to distances r which are 
larger than p. However, m is now no longer an 
arbitrary constant of integration. The solution 
(2.5) has to be adjusted by the law of continuity 
to the inner field. Thus the value of m is de- 
termined by the structure of the inner field. 

The problem of the mass seems thus hopelessly 
removed to a region into which our present 
mathematical capacities cannot penetrate. Owing 
to some unknown properties of the generalized 
field equations, the regular solutions of these 
equations should yield under all circumstances 
a positive m. 

The present paper intends to show that the 
necessarily positive character of the mass may 
be explained, even if we have no detailed 
account concerning the inner structure of the 
field of the particle. Only one scalar restriction 
has to be made on the inner field. This restriction, 


which leads to fundamental consequences par- 
ticularly for the dynamical behavior of a particle, 
may be obtained from the following consideration. 

Let us assume that the final field equations of 
nature are deducible from an invariant action 
principle, which shall contain the curvature 
quantities Rx quadratically. The only invariants 
which contribute to such an action principle, 
are R? and RxR*. The other 3 algebraically 
possible invariants do not add anything new.! 
Now the field equations deducible from such an 
action principle lead to an exact equation for 
the scalar curvature R alone. We get 2 


AR=0, (2.7) 


where A is the invariant four-dimensional 
Laplace operator. This equation (potential 
equation) has no regular solution outside of 


R=const., (2.8) 


if a static field is required. Moreover, considering 
the fact that the metrical field far from matter 
must be practically Euclidean, the constant of 
the right side must be of cosmological smallness. 
For atomistic purposes we can put: 


R=0. (2.9) 


This relation, which has to hold for the entire 
field, implies for the matter tensor the vanishing 
of the scalar 7, i.e., the diagonal sum 7,*, 
sometimes called the ‘‘scalar of Laue’’: 


T=0. (2.10) 


The Maxwellian stress tensor is an example for 
such a matter distribution. 

On the other hand, the so-called ‘kinetic 
stress tensor” 


T= (2.11) 


does not satisfy the condition (2.10). This, 
however, cannot be construed as an argument 
against the plausibility of that condition. The 
form (2.11) of the matter tensor cannot be taken 


1C. Lanczos, Ann. of Math. 39, 842 (1938). At this 
occasion the author wants to acknowledge that he was not 
the first toestablish this theorem. As Professor F. Juettner of 
Breslau has informed him in a letter, the theorem has been 
moe before by F. Juettner, Math. Ann. 87, 270 (1922), 
d on previous results of R. Bach, Math. Zeits. 9, 110 
(oaa} and some Academy papers of R. Weitzenboeck 
2 C. Lanczos, Phys. Rev. 39, 716 (1932); Eq. (5.3). 
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literally as a microscopic distribution of matter 
inside of a particle. It has only statistical 
significance for the particle as a whole, as it is 
clearly indicated by the velocity vector «', which 
has no microscopic significance. Now for an 
averaging process over the entire particle the 
condition (2.10) loses automatically its signifi- 
cance because a nonlinear differential operator 
does not preserve its validity if a statistical 
process takes place. Just as the average value 
of the square of the velocities of a gas is not 
equal to the square of the average velocity, so 
the relation (2.10) breaks down if the micro- 
scopical matter distribution of the particle is 
replaced by an average distribution. One should 
not think that in this process the quantities of 
second order are negligible. As this paper will 
show, the mass of the particle, which is just the 
decisive quantity of the kinetic tensor (2.11), 
drops down to a quantity of second order, if 
the condition (2.10) is to be fulfilled inside of 
the particle. 

Apart from the quadratic action principle, 
the condition (2.10) seems to be an interesting 
hypothesis from another view point. The com- 
plete Riemannian curvature tensor R iim», Vanishes 
only for Euclidean manifolds. Einstein has 
applied a contraction on the Riemann tensor, 
arriving at the tensor Ri. The vanishing of this 
tensor requires much less than the vanishing of 
the complete curvature tensor. But even Rj, 
vanishes only outside of matter, while inside of 
a material particle even this tensor is different 
from zero. It seems interesting to apply a second 
contraction, now arriving at the scalar R, and 
assume that at least this scalar curvature will 
vanish quite generally, not only outside but also 
inside of matter. 

The present paper draws the consequences of 
this hypothesis, without restricting in any other 
way the microscopic distribution of the matter 
tensor. 

We enumerate here all the basic assumptions 
of the following deductions. 

(a) The hypothesis (2.10) is the only specific 
hypothesis which goes beyond the customary 
postulates of relativity. All the other assumptions 
are dictated by strong physical plausibility. 

(b) It is assumed that all field quantities are 
continuous and differentiable. Points of singu- 


larities in which the gj, or their derivatives get 
infinite or even discontinuous, are principally 
excluded. The relation between the curvature 
tensor Rx and the matter tensor of physics is 
assumed according to the fundamental principle 


of relativity: 
Tix = (2.12) 


tx is here the density of the physical matter 
tensor, while x is a gauge factor whose value 
depends on the units employed. In c.g. units 
k=1.87X10-*' cm/g. 

It was shown in an older investigation of the 
author’ that the field equations (2.12) can be 
integrated by a successive approximation process 
which expands the field quantities gi, and also 
Tix in powers of the gravitational constant x, in 
which «x is considered a variable parameter: 


(2.13) 


The only condition of this integration process is 
that the conservation law of the matter tensor 
has to be satisfied, which is a mathematical 
identity of Riemannian geometry : 


V a 


This relation yields a vectorial condition for 
each one of the successive matter tensors 74, 
tix, «++. That means 4 equations between 10 
quantities, which can always be satisfied without 
the danger of over-determination. 

(c) It is assumed that the value of « is small 
enough for sufficient convergence of the series 
(2.13). In view of the smallness of « this assump- 
tion seems quite justified, even if we do not 
know what the actual distribution of matter 
inside of the particle is. We assume that the 
quantities of second order are negligible in 
comparison with the quantities of first order and 
the quantities of third order negligible in com- 
parison with the quantities of second order. 

(d) Outside of a certain practically very small 
action radius p of the material particle, the 
matter tensor shall vanish. 

(e) The field shall be static, i.e., the giz shall 
not depend on the fourth variable x. 


3C. Lanczos, Zeits. f. Physik 13, 7 (1923). 


(2.14) 
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From these postulates the following two 
results will be deduced: (A) the total mass of a 
particle comes out as a quantity of second order; 
(B) that quantity comes out as_ necessarily 
positive. 

In the following formal developments an 
operation will return rather frequently, for which 
a suitable abbreviation will be useful. We want 
to introduce the following notation: 


{X}. (2.15) 


(For infinitesimal fields g=1 and the factor \/g 
may be omitted.) The integration on the left 
side shall be extended over the entire volume of 
the material particle, i.e., over the entire region 
in which the 7 are different from zero. On the 
boundary of the region, and outside of it, 7, =0. 


3. INFINITESIMAL FIELDS 
_ For infinitesimal fields of first order we put 


and normalize the reference system in the 
customary manner: 


—-—=0. (3.2) 
OX, 2 Ox; 


We then have, according to Einstein: 
— Six, (3.3) 
but in consequence of (2.10) we may write: 
Ayin=2Tix. (3.4) 
A is here the ordinary Laplace operator: 
A=0?/dx_" (a=1, 2, 3), (3.5) 


the index value a=4 being omitted, owing to 
the static character of the field. 

For sufficiently large distances r from the 
origin the potential functions yx can be expanded 
in negative powers of r. The lowest term is 


We now want to show that all {7%} have to 
vanish. 


The conservation law (2.11) of the matter 
tensor simplifies at present to: 


OT (3.7) 


Let us multiply by x, and write Eq. (3.7) in 
the following form: 


O(x.T ia) /OXa= (3.8) 


The index & is here restricted to the values 1, 2, 3, 
just as the sum index a, while the index i may 
take all the 4 values 1 to 4. 

We multiply Eq. (3.8) by the volume element 
dr=dx,dxedx; and integrate over the entire 
particle, applying on the left side the Gaussian 
integral transformation. This yields at once 


{Tix} =0, (3.9) 


considering the vanishing of 7, on the boundary. 
The only exception is | 744} which is ordinarily 
not zero and responsible for the mass of the 
particle, 74, being in physical interpretation the 
mass density. 

However, owing to the condition (2.10) we 


have 
=0. (3.10) 


Applying operation (2.4) to this equation we 
obtain 
=0. (3.11) 


We notice that the hypothesis (2.10) necessitates 
positive and negative mass densities which balance 
each other and produce the total mass zero. 

Although the solution (3.6) vanishes com- 
pletely, an outside field is still possible in the 
form of poles of higher order than first. Practi- 
cally, only the dipoles are of interest. Poles of 
higher order decrease too rapidly with the 
distance to become noticeable. Now for dipoles 
the solution (3.6) has to be replaced by the 
following solution: 


1 
Xal ik —{- (3.12) 


In order to obtain information about the value 
of the new bracket expressions, we multiply the 
divergence equation (3.7) by x.x%m and obtain by 
the same method as above the relation: 


mi} + =O, (3.13) 


= 


— = 
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The index 7 can again assume all the 4 values 1 
to 4 but at first let us abandon the value 4 and 
restrict all the 3 indices i, k, m to the choice 
1, 2, 3. We write down the relation (3.13) for a 
cyclic permutation of i, k, m and repeat the 
permutation once more. Then Eq. (1) plus Eq. 
(2) minus Eq. (3) yields: 


sx} =0, (3.14) 
in view of the symmetry of the matter tensor: 
T T xs. (3.15) 


The relation (3.14) establishes the fact that 
the yx. field vanishes even to poles of second order, 
as long as the indices i, k are restricted to the 
values 1, 2, 3. Moreover, by the same method 
as above, making use of the condition (3.10) 
we can prove 


44} =0 (3.16) 


and thus include even y4, among the vanishing 
quantities. Hence, the expansion of the yx starts 
with quadrupoles. 

However, it is different with the 3 quantities 
vis (t=1, 2, 3). Here only the relation 


(xiT ra} + ia} =0 (3.17) 
can be established, while the 3 quantities 
xa} — is} (3.18) 


which form the components of an anti-symmetric 
tensor, need not be zero. In physical interpreta- 
tion the o, may be considered as the 3 compo- 
nents of the angular momentum of a spinning 
symmetric top which is in stationary rotation 
about the axis of symmetry. Although the total 
momentum is zero (the axis of rotation passes 
through the center of mass), the total angular 
momentum is different from zero. If we identify 
the axis of symmetry with the Z axis, we may 
write our solution in the following form: 


Oy 
a= =(- -), 
4ri Ox 


v31=0, (3.19) 


where the constant o represents in mechanical 
interpretation the “‘spin’”’ of the particle. The 
factor 7 which stands for the imaginary unit 
4/—1, enters the formulae (3.18) and (3.19) due 
to the fact that the fourth variable x, of the 
line element is an imaginary quantity, being 
equal to ict. 

It seems remarkable that in sufficiently large 
distances from the center, where the poles of 
higher order become negligible in comparison 
with the dipoles, the spin quantities yis are the 
only existing quantities of the infinitesimal field 
of first order. 


4. THE FIELD OF SECOND ORDER 


Spherical symmetry. Now the field of second 
order shall be included in our considerations. 
Since the scalar condition (2.10) is the salient 
point of our conclusions, we try to integrate 
that differential relation. We do not attempt to 
do that exactly. This is probably impossible. 
But for weak fields the integration is possible in 
successive steps, and only two steps are necessary 
for the conclusions in which we are here inter- 
ested. At first the special but important case of 
spherical symmetry will be treated. The next 
chapter deals with the problem of an arbitrary 
matter distribution. The results are more con- 
clusive in the special case, as one may expect. 

The line element of a spherically symmetric 
metrical field of four dimensions may be written 
as follows: 


ds? (4.1) 


p and q are here some arbitrary functions of r. 
For this line element the nonvanishing compo- 
nents of Ry, come out as follows : 


2 


R22 = = (1—e??+rqg—rp), 
r 
2 . 
ip). (4.2) 
r 


* See, e.g., A. S. Eddington, The Mathematical Theory of 


Relativity (Cambridge, 1 1924), p. 82; W. Pauli, ‘‘Relativi- 
taets-theorie,” Encyk. der math. Wiss. Vol. 19, p. 728. 


MASS 


For sufficiently weak fields e?? may be replaced by 
1+2p+2p’, (4.3) 


if we neglect quantities of third and higher 
order. In this approximation we obtain 


—2p 


[r°g+2r(q—p) —2p 
— 2p?+(q?— gp)r?]. 


The linear terms within the bracket may be 
written thus: 


R= 


2 


r 
(4.4) 


d 
—(r°q—2rp) (4.5) 
dr 

and the condition (2.10), which is equivalent to 
(2.9), results in the following differential equa- 
tion of first order: 


d 
—2rp) =2p?—(u?—up)r’, (4.6) 
r 


where we have put 
(4.7) 


The two arbitrary functions u and p are reduced, 
in consequence of (4.6), to but one arbitrary 
function. It will be convenient to consider u as 
the arbitrary function, and to express p in terms 
of u. 

In first approximation the right side of (4.6), 
being a quantity of second order, may be 
neglected. We thus get: 


ru—2rp=A. (4.8) 


Putting r=0 we notice that the constant A has 
to vanish since singularities are not permitted. 
Hence we obtain: 

p= (4.9) 


For the second approximation it is permissible 
to make use of the relation (4.9) on the right 
side of (4.6). This yields: 


—(r?u — 2rp) = (4.10) 
dr 


Integrating, we have 


r°u—rp=} f uur'dr 
9 


wrdr 411) 
0 
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On the right side u needs to be known only in 
first approximation. But the results of the 
previous chapter have shown that no spherically 
symmetric outside field exists in first approxi- 
mation. Hence u has to vanish outside of the 
particle, although it may take arbitrary values 
inside. Consequently for any r=p the first term 
on the right side of (4.11) vanishes, while in 
the second term the upper limit of integration 
may be replaced by the constant p. We thus 
obtain for any r=p: 


(4.12) 
0 


On the other hand we know, from the theory 
of the Einsteinian vacuum equations for in- 
finitesimal fields, that outside of the particle we 
must have 


e2?=14(m/r), e2¢=1—(m/r). (4.13) 


We know in advance that for our problem m 
must come out as a quantity of second order. 
We can put practically: 


p=m/2r, qg=—m/2r, 


4.14 
u=qg=m/2r* 


and the substitution of these values in the left 
side of (4.12) gives: 


m=} f 
0 


This equation establishes the mass as a necessarily 
positive, second-order quantity. 

One of the fundamental theoretical results of 
general relativity is the strict proportionality 
between the mass of gravity and the mass of inertia. 
This result is not a natural outcome of some 
general principles, and we have to discuss this 
question, in view of the second-order quality 
of the mass. 

An application of the conservation law (2.11) 
to the dynamics of a pariicle shows® that the 
inertial mass of a particle has to be defined by the 
following quantity: 


m=({T}. (4.16) 


5 C. Lanczos, Zeits. f. Physik 59, 514 (1930). This 4 4 
covers satisfactorily only the case of special relativity. But 
the method can be properly —— to general relativity, 
as the author hopes to show in a later publication. 


(4.15) 
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As long as the mass is infinitesimal of first order, 
covariant and contravariant quantities need not 
be distinguished. The equations for infinitesimal 
fields yield: 


1 1 m 
(4.17) 
4 r 


r 


which gives at once the relation: 


m=4em. (4.18) 


Gravitational mass and inertial mass are identi- 
cal, except the factor 47. 

For the situation here considered the above 
proof does not hold. Both m and m vanish in 
first approximation. If we want to show their 


‘proportionality, we must take into account the 


field with second-order accuracy. We have to 
utilize the exact expressions (4.2), neglecting 
only quantities from the third order on. 

The last expression of the equations (4.2) 


yields 
2 

r 


Applying the operation (2.14) on this equation, 
we obtain with an accuracy of second order: 


d 
TH} =m=4e f 


0 
Pup jr. (4.20) 


Now in the first term under the integral sign 
we can integrate by parts and obtain: 


made f (4.21) 


In the first term we can substitute for w its 
value according to (4.14), while in the second 
term the relation (4.15) may be utilized. We 


thus obtain finally: 
m= 222m /3. (4.22) 


The strict proportionality of inertial mass and 
gravitational mass is once more established, al- 
though the factor of proportionality has changed 
from 4x to (22/3). 


5. THE FIELD OF SECOND ORDER. ARBITRARY 
STATIC DISTRIBUTION OF MATTER 


In this last chapter we follow the same line 
of thought as in the chapter before, but aban- 
doning the restriction of spherical symmetry, 
The metrical field shall again be static. It shall 
also be sufficiently weak to permit first-order and 
second-order approximation, neglecting quanti- 
ties of third and higher order. But no further 
a priori restrictions shall be made. The condition 
(2.10) shall again be integrated, with first- and 
second-order accuracy, in a similar manner as 
before. 

The reference system shall be normalized by 
a condition which is very useful in all investiga- 
tions requiring the exact values of the Ry.® 
This condition is: 


=0, (5.1) 
which may also be written in the form 
(5.2) 


For infinitesimal fields the condition (5.2) 
reduces to the condition (3.2), used for the 
normalization of infinitesimal fields. 

In the reference system (5.1) the following 
expression holds for the scalar curvature R: 


og’ 
(5.3) 


R=A logy/g+2l w* 
Xa 
A symbolizes the invariant Laplace operator of 
4 dimensions: 


(5.4) 


which reduces, however, to 3 dimensions owing 
to the static character of the field (u, v=1, 2, 3). 

The condition (2.10) requires the integration 
of the following differential equation: 


A logy/g= — 30 w*dg""/ Axa. (5.5) 


Again we do not attempt an exact integration 
but are satisfied with successive approximations, 
stopping after the second step. 

In first approximation the right side of (5.5) 
may be omitted and we get the ordinary po- 
tential equation for log \/g. The only possible 


®6C. Lanczos, Physik. Zeits. 23, 537 (1922). Equation 
(5.3) above quotes Eq. (11) of that paper. 
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solution of that equation free of singularities is 
a constant. Hence 
g=const. (5.6) 


and the boundary conditions in infinity specify 
that constant to 1: 
g=l. (5.7) 


For the second step let us put 


since y must be a quantity of second order. For 
this quantity the following differential equation 
has to be solved: 


Ay = (5.9) 


Neglecting quantities of third order we may 
write this relation as follows: 


(+ OYva 
Ay 


OY wa OY w OY wv 
Ox, OXe OX 


Ox, OXa 


On the right side the yi, need to be known with 
first-order accuracy only. In this approximation 
the normalization of the reference system reduces 
to (3.2) and simplifies even further to: 


OY ia/IXa=0, (5.11) 


owing to the vanishing of y. Hence the first 
term of the right side of (5.10) may be written 
in the following form: 


OY 
Yua ) (5. 1 2) 
Ox, OXa 


Let us now multiply Eq. (5.10) by the volume 
element dr and integrate over the infinite space. 
Then the term (5.12), being a complete diver- 
gence, and in view of the boundary conditions 
in infinity, will not contribute anything to the 
integral. The left side too permits the application 
of the Gaussian integral transformation. If we 
denote by N the inward normal of a sphere 
with an infinite radius, by do the surface element 
of that sphere, we obtain 


OY uw» 
OXe 


(5.13) 


Now we know that for very large distances 
the dipoles and poles of still higher order become 
negligible, and the line element must reduce to 
the following form: 


+(1—m/r)dx.?. (5.14) 


This yields: 
y=2m/r (5.15) 


and the relation (5.13) gives: 


1 
( ) dr. (5.16) 
4 OXe 


Here is the expression that replaces the simpler 
relation (4.15) found before for the case of 
spherical symmetry. It is of interest to assume 
at first a line element which is positive definite 
in all 4 dimensions. Then all the y% are real and 
the relation (5.16) shows that any deviation from 
Euclidean geometry, wherever it occurs in the field, 
will contribute something positive toward the total 
mass of the particle. 

The actual line element of nature has 3 
positive and 1 negative dimensions. Here the 
three yi4 quantities (i=1, 2, 3) are imaginary and 
their square is negative. All the yu with i,k 
=1, 2, 3 and also ya, contribute something 
necessarily positive to the total mass. But the 
“spin quantities” yi, have a decreasing instead 
of increasing influence on the mass. It is feasible 
to assume that the spin of the particle is a kind 
of correction which cannot be so powerful as to 
change the sign of the mass from positive to 
negative. The spherically symmetric particle 
considered before has, of course, no spin. There 
we did not encounter any terms which had a 
decreasing influence on the mass. 

Let us investigate again the question of 
proportionality of the mass of gravity and the 
mass of inertia. The situation in the general 
case is not as favorable as it was before for the 
spherically symmetric particle. By the same 
method we have applied before, the following 
expression can be deduced for the mass of 
inertia, omitting the details of calculation: 


“) dr. (5.17) 


716 C. LANCZOS 


For spherical symmetry the vis are zero, except 
as. Hence the second term reduces to 


or using the notations of Section 4: 
8x f rutdr. (5.19) 


This integral is reducible to m owing to the 
relation (4.15) and we arrive once more at (4.22) 
which expresses the mass of inertia in terms of 
the mass of gravity. But for the general case we 
do not succeed in reducing the additional integral 
of (5.17) to m alone. The spherically symmetric 
case has only two degrees of freedom and the 
additional restriction T =0 leaves but one degree. 
The general case, after normalizing the reference 
system, has 6 degrees of freedom, reduced to 5 
by the vanishing of the scalar 7. Here the 
proportionality of inertial and gravitational mass 
cannot be a priori established, without additional 
restrictions relative to the distribution of matter 
inside of the particle. 


6. CONCLUSIONS 


The electrical action radius of the electron is 
of the order of magnitude 10- cm, while the 
mass of the electron has the amazingly small 
value of about 10-*5 cm. In view of the almost 
negligible part that gravitational forces play in 
comparison with the electric forces, the con- 
jecture was frequently made that gravity must 
be a secondary effect, a kind of surplus phe- 
nomenon of electrical action. The present paper 
gives a theoretical reason why the mass of a 
particle may have to be considered as a quantity 
of second order. If the distribution of matter is 


subjected to the restriction that the scalar T of 
the matter tensor shall vanish, we get at once 
that the total mass of a particle is a quantity 
of second order. More than that, we are also 
able to account for the fact that the masses of 
nature are always positive. It turns out that all 
the y« which characterize the deviation of the 
metrical field from the Euclidean field, con- 
tribute something necessarily positive to the total 
mass, with the only exception of the yi4 quantities 
(¢=1, 2, 3) which are only present if the particle 
has a stationary spin. The spin part of the 
metrics is the only part which has a decreasing 
instead of increasing influence on the mass. 

The really fascinating consequences of the 
possibility that the total mass of a particle may 
be a quantity of second order, belong to another 
field, which has not been touched in the present 
investigation. Customarily it is assumed that the 
dynamical law associated with the Einsteinian 
field equations Rx%=0 is necessarily the law of 
the geodesic line, no matter what happens with 
the matter tensor inside of the particle. This is 
actually true as long as the mass is a quantity 
of first order. But an entirely new situation 
arises if the value of the total mass drops down 
to a quantity of second order. The dynamical 
law is then no longer of second but of third 
order, and is much more comprehensive than 
the simple motion along the geodesics might 
suggest. The author was able to show that even 
the electrostatic interaction of particles is within 
the scope of the generalized dynamical law (leaving 
apart at first the general electromagnetic inter- 
action), in spite of the fact that the Riemannian 
basis of geometry is retained, the 4-dimensional 
concept is retained, and even the Einsteinian 
field equations for the empty space are retained. 
The author hopes to be able to publish these 
investigations on some other occasion. 
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The levels 3d* *GDH (and ?P?) and 3d? 4p(?F) “GFD (@P) ‘DPS ('G) *HG have been 
identified in Co VII and Ni VIII. The graphs of relative energy against Z for the even levels 
are, on the whole, smooth; however *G 2H 2D show marked and very similar irregularities, which 
remain unexplained. The doublet separations behave in a regular fashion. The sharp decline in 
certain of the odd levels, especially ‘+FDDP, between Co VII and Ni VIII can be explained as 
due to the approach of (undiscovered) levels of the 3p°4d‘ configuration, which are expected to 


cross 3d*4p near this stage of ionization. 


I. INTRODUCTION 


EVERAL of the complex spectra of the 

elements of the first long period with high 
excitation have been partly classified. Wherever 
intensive work has been done on any sequence, 
it has been possible to classify some of the levels 
in the configurations which give rise to transi- 
tions 3d*—3d*—4p in stages of ionization up to 
about VI. Beyond that little has been done. 
Bowen! has pointed out that at about this stage 
of ionization the 3p°3d** levels are expected to 
cross and hence disturb. the 3d*~'4p levels. This 
enhances the difficulty but also the interest in 
the study of the 3d*~'4p levels in this neighbor- 
hood. The present paper extends, through Co VIT 
and Ni VIII, the work of Bowen and others? on 
the three-electron scandium-like spectra. 


II. EXPERIMENTAL 


The plates for this investigation were taken 
on the University of Wisconsin’s 21-foot, 30,000- 
line-per-inch, grazing incidence vacuum spectro- 
graph with a condensed spark source. The 
apparatus and the method of reducing the 
readings to wave numbers are described in detail 
elsewhere.? Independent readings were taken 
upon two plates for each element. The excellent 
agreement between the two sets can be at- 


11. S. Bowen, Phys. Rev. 53, 889 (1938). 

?Sc I: H. N. Russell and W. F. Meggers, Sci. Pap. Bur. 
Stand. 22, 329 (1927); Ti II: H. N. Russell, Astrophys. J. 
66, 1 (1927); V III: H. E. White, Phys. Rev. 33, 674 
(1929); Cr IV: H. E. White, Phys. Rev. 33, 674 (1929); 
I. S. Bowen, Phys. Rev. 52, 1153 (1937); S. Pasternack, 
Astrophys. J. 92, 129 (1940); Mn V, Fe VI: I. S. Bowen, 
Phys. Rev. 47, 924 (1935). 

3J. E. Mack and E. E. Anderson, in preparation for 
Rev. Sci. Inst. 


tributed to the large dispersion (4500 cm~! per 
mm at 1.6X10® cm to 181 cm™ per mm at 
210° cm-'), the use of a rigid plateholder 
accurately ground to shape, and the use of 
polished plate glass, 5 feet in length, for the 
Schumann plates. The wave numbers of the 
lines were determined relative to the copper 
lines of Betz,* which in turn had been determined 
relative to Edlén’s oxygen lines on the same 
spectrograph. Each plate had on the order of 3000 
lines in the wave number region extending 
approximately from 219,000 cm to 926,500 
cm~!. The intensities were estimated visually on 
a crude logarithmic scale of 0 to 10. 


III. RESULTS AND DISCUSSION 


The results of this investigation are given 
below in Tables I and II for Co VII, and in 
Tables III and IV for Ni VIII. For the con- 
venience of any reader who wishes to make a 
critical study of the levels, we present the 
results in the form of rectangular arrays. The 
row and column headings in Tables I and III 
comprise tables of levels. The information given 
in the row-column intersections consists of three 
items: intensity of the line (bold face type); 
discrepancy, vobs MINUS Peale (italics); and a 
number (upper left-hand corner) which corre- 
sponds to the position of the line in Table II or 
IV. The symbols mean: b—broad, d—diffuse, 
c—classified elsewhere in this spectrum. 

Tables II and IV include only lines classified 
in this investigation.® 

‘J. C. Betz, Thesis (1935), University of Wisconsin 


ibrary. 
5 We shall be glad to make our complete line list available 
to anyone who may wish to use it. 
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TaBLe I. Co VII. 


3d3 2G 2D 2H 
24 3} 33 43 2 4} 
v 0 698 1610 2723 21096 «21304 22187 24151 25063 31348 231555 33251 33873 
Av 698 912 1113 18373 208 883 1964 912 6285 207 1696 622 
75 122 
1205 70 7: 133 
348 72 76 79 125 
13} 447,314 7-2 0-1 0-4; 
737 65 69 73 116 
|23 448,051 9c 0 10 0 §+5 We — 
1254 68 
{4} 449,305 80 
135 62 67 71 114 121 
|33 449.440 3+45 10+1 8-4 0O+4 
1425 56 108 117 131 
3d2(3F 2F (25 450,865 10c +33 3-i — 10c + 9+4 — 
124 63 115 119 132 
(45 450,989 0-5 0¢+1 14+14 10ce —12 0-2; 
226 | 50 105 109 129 130 
3d2(3F)4p 451,215 +11 9+6 §-30 — lc +40 1-3 34401 
344 48 104 128 
14 451,559 10+9 2-10 2-34 
251 59 
199 57 107 127 
557 38 44 52 lll 
133 37 43 101 102 
3d2(3F 2D(1} 452,709 1-2 3-17 2+9 1413 — 
1448 39 49 100 1 110 124 126 
(33 454,157 4-16 0-3 14+24 9-7 tc-26 4412 0412 
325 31 36 99 23 
4688 22 24 25 96 97 112 120 
3d2(3F )4p 2G(34 459,170 1411 O-17 O-18 10+3 2-1 5-95 —5 
1417 21 23 OF 95 113 118 
(43 460,587 1-18 2-36 00 10-1 2+1 
4703 17 80 81 98 
3d2(3P)4p 48 14 465,290 — 3bd —36 3+4 4-1 3-4 o-i1 — 
4068 16 64 66 93 
3d2(3P)4p 2469,358 lbd 0 Oc+5 9-3 
1017 13 15 58 60 92 
{14 470,375 § —22 9+ Oc+9 — 1414 
1505 12 14 4 90 91 
{25 471,880 3+2 0-2 0+9 44+40 1-17 
2187 9 10 ll 45 61 87 8y 
474,067 0-4 2+10 3c —8 0-1 0+9 
508 6 34 85 
3d2(8P)4p 4P(} 474,575 | 0+4 4+3 10-3 2c +17 
581 7 29 30 5 82 
475,156 3c —17 0415 4-1 2-16 2-1 
898 4 46 86 
3A('G)4p 476,054 2+4 x-3 — 8+1 
234 5 42 51 84 
(43 476,288 0-21 — 14+41 5c+1 2c 0 10 0 
304 3 27 40 8 
(24 476,592 0+17 2+4 6+4 8+15 0-8 10 —23 
9024 | 2 19 41 47 
3d2('G)4p 2H(44 485,616 | 0+30 — §5+10 10 4-5 0-4 
1968 1 8 28 32 
(5) 487,584 4+26 1+2 0+10 4-2 


The levels were predicted by extrapolation in 
the Sc I-like sequence.? In the curves, Fig. 1, 
we have plotted »(Z+1)—»(Z) against Z, a 
scheme which allows irregularities to be shown 
on a larger scale than a curve of v against Z. Here 
the differences are arbitrarily referred to that 
for the lowest level in each configuration, i.e., 
3d? *F,, and 3d*4p ‘G2, appear as horizontal lines 
at ordinate zero. The separations 3d°4p(°F) 4G2; 
—3d*‘F,, and their first, second, and third 
differences are as given in Table V. This investi- 
gation contributes only the last two columns of 
points on the curves. Irregularities in the other 


columns are probably in the main due to the 
interconfiguration interactions (3d*, 3d*4s, 3d4s*) 
and (3d*4p, 3d4s4p, 4s°4p), with perhaps one or 
two errors in the literature, as in III d*?P. 
Among the (even) d’ levels, the ‘F curves show 
a smooth general upward trend, and the ‘P 
curves a general downward trend with slight 
irregularities in ‘P, and *P2,; between Co VII and 
Ni VIII. The curves are so smooth that it 
appears possible to predict levels by extrapola- 
tion, to an accuracy of about 100 cm~'. However, 
the */7 and °G curves show a remarkably sharp 
rise of almost 1000 cm-' between Fe VI and 
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TABLE II. Classified lines of Co VII. 
CLASSIFICATION CLASSIFICATION CLASSIFICATION 
INT. 3d8 INT. » 3d3 3d24p INT. v 3d3 
4 484,887 ‘Fy — 2 431,622 ‘Py, — 
0 484,036 ay | 10 450,959 — @F)Dy 1 431,418 
0 476,609 — GPP» | 9 450,523 | 4 430,388 4 — 
2 475,360 — (G)PGs, 0 449,702 — | 2 430,245 
0 474,657 0 449,374 ‘Fy, — (§F)*Fa 5 430,089 
0 474,579 — — Fey 9 429,999 — 
94 10 449,288 = 34 9 429,829 ‘Py Fs, 
3 474,441 
2 472,467 3} 34 8 449,087 \4Py 1 429,068 34 
10 471,355 4Fy — | 0 448,994 — | 1 428,403 — 
3 | 3 448,747 — sy 5 — 
5 4 0,353 14 a = 2 42 
10 470,268 |? 48067 { ‘Py | 0 427,265 — 
4 
9 469,684 2h 9 448,051 {api 1 426,852 rT] 
ibd 469,358 *Fiy ‘Py — Fy 
3bd 464,556 — | 10 447,831 — (F)*Fs | 10 426,724 — Fy 
1 462,523 — (G)?Hs | 8 447,695 33 — — 
5 461,475 37 43 | 10 447,353 — Fa, 10 425.914 { — Fay 
1 460,553 a4) | 10 446,971 4 — 
1 458,959 — 8 446,713 — @F) Fs; 0 425,293 — Fy 
1 458,483 3} 7 446,614 3 — 5 423,569 ‘Py — 
2 457,828 a7 4} 5 446,446 3) — 3 422,927 ?2Dy 
0 457,543 34 3 446,263 1} 4 422,821 — 
0 456,429 a7 3} 7 445,761 0 422,259 °Ga — 
2 455,292 4P,, — (@P)*P | 0 445,703 — 0 420,918 — 
6 454,409 23 24 0 445,236 — | 0 420,650 — (@F) Fs 
0 454,343 — ('G)?Hs, |10 445,014 2Dy — | 2 419,970 ?Dy— (@F)D, 
0 454,075 4P, — | 0 444,587 ‘Fy — 1 419,864 — 
4 453,851 ais 444,198 — @P)Sy 3 419,691 
4 453,786 4Fy, —?(F)?Dy | 4 443,985 — @P)Sy 9 419,521 — 
4 453,709 — (G)*Hs, | 2 443,600 — @P)Py | 0 417,762 "Hy — Fay 
4 453,482 4P, — (8P)*P, 3 443,099 0 415,377 2Dy— 
2 452.953 | 443037) — Py 
1d 452,871 ‘Fy, —?@F)?D | 8 442,804 — (G)Gy 
1 452,707 — GFPDy | 0 442,728 — (@P)*Dsy 
2 452,571 — |10 442,415 — 
4 452,531 3} 1 440,853 "Hy — @P)Ds 
8 452,456 — @P)*Px | 4d 440,573 — 24 
+ 452,360 — (GPHy 1 440,308 @P)*Dy 
1 452,178 — 1 438,834 
3 451,994 — 2 437,780 4 
— | 0 436,436 °G3, — 
3 451,872 — (@P)*D3, | 10 435,523 
— 
10 435,022 °Gy, — 
0 451,739 — (G)Hy | 2 434,106 “a7 3} 
10 451,568 — @F)D, 0 433,941 — @P)Sy 
10 451,431 a7 34 2 432,295 
5 451,226 { agit = (6)Ge 


Co VII, and a fall between Co VII and Ni VIII 
such that the difference between Co VII and 
Ni VIII is what would be predicted from Cr IV, 
Mn V, and Fe VI. In other words, considering 


vy as a function of Z, each of these four levels 
shows a shift in Co VII with respect to what 
would have been expected by extrapolation from 
lower stages of ionization, and stays shifted in 
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TABLE III. Ni VIII. 
J 2 3 4} 4% 4 5 
0. tale 3721 23261 23710 24669 26977 28068 34689 35120 36754 
Av 1018 «1178S «1637 19540 449 959 2308 1091 6621 431 1634 781 
64 71 - 119 
3d2(3F)4p 4G(2} 565,124 6-5 — 0414 — 0+6 
264 63 67 113 132 
3d2(3F)4p ‘F(1} 565,388 8+9 6b-—9 0-16 — — 4413 
1443 59 62 66 108 sd 13 
12} 566,831 20 5-2 90 7-16 0+7 — +32 
133 61 110116 1 
13} 566,964 34+1 — O-1 1423 — 0+12 
1782 53 60 65 104112 115 129 
568,746 5+13 4+8 2-2 2-6 O-11 2+13 0+40 
818 56 114 
103 41 97 98 126 127 
3d2(3F)4p 2D(1h 569,667 2%+20 — 2-8 O41 2-11 0-7 
172 40 48 95 101 125 
569,839 14+2 0+4 3+4 2-22 
514 38 91 94 
13 570,353 6bdc +17 1-7 2-12 _ 
193 51 93 107 124 
3d2(3F 2F (24 570,546 8-1 %+7 — 10 — 0+2% 
414 49 105 28 
(44 570,960 100+4 9-3 045 — 
557 37 46 89 92 
[23 571.517 2-2 6+7 2-2 — 
287 43 52 106 121 
(33 571,804 — &+1 60 — — 3+3 
41 32 42 90 120 ©6122 
(24 571,845 3-10 — %+26 — 0+12 0+8 
1124 102 123 
(54 572,969 6 0 4+24 — 
358 87 100 118 
(33 573,327 — 1-15 — 1413 3444 
8010 19 20 85 103 
3d2(3F)4p °2G(3} 581,337 1+5 2-19 — 0+5 
1904 18 1 83 9% 99 
(43 583,241 2-3 2+6 5-4 2-5 0-4 
4064 15 68 70 86 
4S 14 587,305 1-s — 3+8 3-7 247 — 246 
3459 13 54 58 82 
3d2(9P)4p 4 590,764 20 Oc —4 1-2 
1411 10 12 50 55 
{1} 592,175 2-13 6-7 ~ 2+22 Oc —7 
2 9 11 39 45 9 80 
{2} 594,068 0-3 O+2 1-5 6bde +12 1414 — 6-12 §+2 
2702 26 28 
3d2(8P)4p 4 596,770 6 0 00 
135 27 29 74 
3 30 75 7 
{3} 596,908 40 x0 — — 3-8 
1662 2 3 23 2234 
(23 598,570 3-8 1-11 — 2+8 8+10 2-10 
68 4 6 4 33 36 69 76 
3d2(1G)4p 2G(3} 598,638 — 0-14 1+33 O+17 8+17 24+1 2+10 1-5 
441 5 31 35 733 «O77 
(44 599,079 2-2 4-6 6+2 
14338 1 16 17 22 
3d2(1G)4p 2H(44 613,417 e-s — 5+2 8+8 
(54 615,725 — 3-8 — 1+8 


Ni VIII. The *D levels behave somewhat less 
regularly throughout the sequence, exhibiting 
the same general behavior as *G and ?H in 
Co VII and Ni VIII except that the shift is 
opposite in sign. This shift is remarkable since, 
so far as we know, it is the first instance of such 
an irregularity that cannot be attributed to a 
mutual interaction. Other configurations than 
d* can hardly account for the disturbance since 
the nearest configuration of the same parity, 
3d*4s, is already distant by some 300,000 cm-', 
a distance that increases smoothly with Z. 
Within the configuration d*, on the other hand, 


the only levels yet to be found are 14, 
°F 24,33, the J values of which prohibit them from 
influencing *G,,; and *H. Since the spectra in 
this region show a crude qualitative agreement 
with the Slater theory of complex spectra, an 
attempt has been made to correlate the shift 
with a possible change in the Slater® integrals 
F,(3d3d) and F,(3d3d). However it appears 
impossible self-consistently to make such a 
correlation. The LS-approximate values calcu- 
lated for the known levels with respect to the 


6 J. C. Slater, Phys. Rev. 34, 1293 (1929). 
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TABLE IV. Classified lines of Ni VIII. 
CLASSIFICATION CLASSIFICATION CLASSIFICATION 
0 611,230 ‘Fy — | 8 568,361 ‘Fy — @F)Fy | 2 545,148 ‘Py — 
1 3 —-$98,562 4 — | 6 568,083 | 4 544,925 — 
1 597,547 4 | 5 567,747 4 — | 0 $44,588 — 
0 596,440 — | 567,499 {*Ps — | 2 544,071 ‘Px, — 
2 595,356 “a 4 | 9 543,980 “ 
1 594,950 ‘Fy 3 567,380 ‘Fy, 3d 543,731 GF)*Fy 
0 594,065 — | 5 567,240 a — | 1 543,569 
593,187 “an 4 1 567,052 ‘Py — | 9 543,121 ‘Py — Fy 
0 593,058 2 $66,831 ‘Fy — | 0 $42,897 Gy 
2 $92,162 566,570 3 | 0 $42,294 — 
1 $91,879 ‘Fx, — (@P)'Dy | 3 565,953 ‘Fx — @F)Gy | 7 $42,178 — (3F)*Fy 
6 591,156 — 5 565,817 — | 541,758 °Gay — 4 
$90,764 8 565,397 4 | 0 541,662 ‘Py — 
587,649 — | 5 $65,119 4 @FYGy | 2 541,471 — 
1 $87,297 — | 2 565,023 — | 2 540,691 — Fy 
5 586,442 0 564,649 — CF | 1 540,010 — 
0 585,342 “an 4 | 6b 4 | 0 539.861 4 — OF) Fy 
2 «581,054 ‘Fy — | 3 564,052 ‘Py — 3 538,652 *Day — 
1 580,330 | 2 563,959 (G)Gy 0 $38,152 — 
~—-$79,134 | 3 563588 ‘Py— | 0 537,168 — 
1 $78,258 — (G)*Hy| 0 562,954 ‘Fy — | 537,118 * Day — 
8 576,671 a | 2 562,643 ‘Px, — | 0 536,733 — 
2 $74,868 — | 0 $62,333 *Hy — | 5) 536,207 *Hay — OF 
0 573,986 9 567 218 | 0 535.450 %Dy — @FPFy 
8 573,911 a — 4 — (@P)*Dy | 0 535,160 4a — 
$61,879 — (G)Gy | 2 $34,967 — 
6 573,509 ‘Py — | 1 561,599 4 | O 534,540 
5 573,192 14 | 3 560,152 a — @P)Dy 1 534,197 "Hy — (@F)*Fy 
0 573,060 | 6 559,367 2D — 2 | 0 534,097 — 
| 2 556,270 — | 1 $31,743 Dy — Fy 
+ 572,096 — ('G)*Gy | 7 555,663 — | 4 530,281 
3 571,835 — @F)?Dy | 5 555,169 — (@F)*Gy 
8 571,678 — | 5 554,369 
2 571,583 3 — | 2 553,269 
6 571,013 a — | 
$52,190 — | 
77 2 570,571 — | 1 548,643 ‘P,, — | 
1-5 2 570,503 — | 6 548,124 — | 
1 6bd 570,370 \ = 1 547,174 | 
1+8 1 569,841 — @F)*Dy | 
2 $69,687 | 3 546,843 
4,24 lay | 2 546,631 y — (F)D, 
1 569,413 ‘Px — 
an 10 568,831 4 — | 0 545,762 — | 
hift 10 568,780 3 (@F)*Fy | 1 545,272 — | 
2 568,487 — @P)'D, 
rals 
ars 
la ‘F are: *Diower = 20F2+ 75F, 
lcu- 4F=0 (193 F2— 1650F2F,+8325F 2)! 
the ‘P=15(F.—5F,) and the correlation would require changes in the 
*G=105F,+4(F:—5F,) first differences of F; and F, to be chosen so as 
*H=105F,+-9( F.—5F,) to vield the above-described behavior. 
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Fic. 1. First differences in the scandium-like isoelectronic sequence. v(Z +1)—»(Z), where each v is measured with I 
respect to the lowest level of its own parity. ( 
i 
TABLE V. their Of the odd levels, the ‘G curves run quite \ 
smoothly, showing only a slight flattening out t 
| At Ae As between Mn V and Fe VI. However, ‘F, ‘D, and 
| | Se I —4741 4P show a sharp decline between Co VII and 3 
1 33,377 Ni VIII. The sextet and quartet FDP levels of 3 
Till 28,636 56,887 23,510 ~8,573 the disturbing configuration 3p°3d‘ undoubtedly 
a) Vill 85,523 14,937 are the lowest and hence, to the left of the 
71,824 — 2,202 
1 Cr IV 157,347 12,735 crossing, will affect most strongly the quartet 
} 84,559 -944 FDP levels of 3d*4p. Hence the sharp decline in 
MnV 241,906 11,791 
96,350 _636 the levels of Ni VIII may be taken as a qualita- 
| | Fe VI 338,256 “ase 11,155 me tive verification of the prediction of Bowen! that 
' Co VII 445,761 , 11,858 these levels cross near this stage of ionization. A 
119,363 preliminary survey of our copper line lists fails 


Ni VIII 565,124 : 
to show any of the expected strong lines of 
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TABLE VI. List of possible levels and their transitions for the 3d**P levels. 


Co VII Ni VIII 
3d3 2Pyy’ 2Py” 2p,” 2p,” 
? 28,876 ? 29,324 ? 29,244 ? 29,712 ? 31,515 ? 32; * ? 32,261 ? 32,939 
4p 0+17 +s 30 
2 
Oc +24 1-4 0 +34 0+20 0+7 
2+29 3-3 10 
24 1c +8 
30 30 Oc 0 10 ac +15 
24 10 1-15 
2P 3-21 0-1 0-1 
Oc +24 0+22 4411 Sc +3 
—14 —{2 
0-30 2+32 2+23 Tc +27 
‘Py Oc +3 Oc +15 1-5 2-8 Oc —32 1-19 3-19 
4 3c —1 4+il 2-14 & +7 Oc +5 2 +5 
24 & —21 0c +5 Ic —3 4+8 


TABLE VII. Level and line count, and discrepancies of observation, in the configurations 3d* and 3d*4p. Note.—Sc | and 
Ti II are omitted. V III and Cr IV have, respectively, 10 and 6 levels of other configurations, giving rise to 54 and 21 lines, not 


counted here. 

Vill Criv Mn V Fe VI Co VII Ni VIII 
Number of levels 33 37* 41* 42* 43 43 
Number of lines 66 46* 80* 106* 133 132 
Mean wave number 0.8 -108 cm" 1.6 2.4 3.4 4.5 5.7 
Mean | Vops — Veale! 1 4 cm" 6 2 4 2 4.2 9 7 8 2 
Number of inter-combination lines 12 13 14 29 48 45 


* Plus 3 levels giving rise to between 9 and 16 lines in each of the spectra Cr IV, Mn V and Fe VI, discovered by Bowen and reported without 


line data by S. Pasternack (see Astrophys. J. 92, 129 (1940)). 


CulIX in the expected positions, probably 
because of this crossing. 

From transitions alone, it was not possible to 
determine the J value of the levels d*?Dz, 
(corresponding to Bowen's level A) since chance 
coincidences give the lines 444,925 and 429,274 
in Co VII and 548,561 in Ni VIII equal to the 
wave number differences of the *D2; levels and 
the *G,, levels. 

A search for the 3d*?P levels resulted in two 
possible sets of values (indicated by single and 
double primes) with the evidence not sufficient 
in any case to establish either set as real. Table 
VI gives a list of the possible levels and their 
transitions. 

It will be noted that in Ni VIII the level 
3€°(°P)4p*P,, (596,905 is equal to 
3d°(°P)4p'D;, (596,908 cm-'). Question marks 


have been placed on several levels where the 
determination depended upon a single transition, 
or upon several doubly classified lines, or where 
the discrepancies | are large. We 
believe all the questioned levels to be real except- 
ing those listed in Table VI. 

From Table VII it will be noted that the 
average relative discrepancy is of the same order 
of magnitude in all the spectra. The fact that 
more lines have been found in Co VII and 
Ni VIII than in earlier members of the sequence 
is due no doubt both to instrumental factors, on 
the one hand, such as the larger spectrograph 
and the excellent sparking conditions, and to 
real atomic properties on the other hand, such 
as the breakdown in LS coupling shown by the 
greater number of quartet-doublet transitions. 
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Paramagnetic Relaxation and the Equilibrium of Lattice Oscillators 


J. H. Van VLECK 
Harvard University, Cambridge, Massachusetts 
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The present paper is rather negative in character. Its aim is to show that the conventional 
models of paramagnetic relaxation cannot be used at helium temperatures, because the lattice 
vibrations cannot possibly serve as a thermostat in the way ordinarily supposed. Namely, 
if the quanta exchanged between the spin system and the lattice oscillators are of the order 
0.2 cm characteristic of the usual Stark splittings in chrome or iron alum, the spins will be 
in thermal contact only with oscillators at the extreme low frequency end of the lattice spectrum 
which are too few in number to conduct away the surplus spin energy. It is found that this 
limited band of oscillators is interrupted much more frequently by interactions with the 
spins than by (a) collisions with the walls or (b) interplay with other oscillators due to an- 
harmonic terms (calculated in the next paper). Consequently the lattice vibrations in thermal 
contact with the spin system will be at nearly the spin temperature, and not at that of the 
helium container as ordinarily supposed. Impurities with abnormally large Stark splittings, 
of the order 2.0 cm=, would make a wider band of lattice oscillators available and so would 
avoid the difficulty of insufficient lattice conductivity, but it is very doubtful whether the 
impurities could be in thermal equilibrium with the great bulk of the paramagnetic ions. 
By thus showing that conventional hypotheses will not work, we aim to pave the way for a 
6 future paper proposing a rather unusual substitute conduction mechanism which is, however, 


if 

f 

| 
| 


compatible with the thermodynamic formulae of Casimir and du Pré. 


Introduction.—In a previous paper! the writer 
computed paramagnetic relaxation times for 
chromium and titanium alum on the basis of the 
thermodynamical model of Casimir and du Pré,? 
wherein the lattice oscillators are treated as a 
thermostat as far as the interchange of energy 
with the spins is concerned. These calculations 


of all oscillators, has such a small heat capacity 
that it can serve as a thermostat only if it is 
surrounded by a bath of liquid helium, and is in 


equilibrium with the latter. However, the situa- | 


tion is much more acute even than this, because 
of the fact that the work required to turn over a 
spin against either the applied magnetic or the 


| were in fairly satisfactory agreement with experi- crystalline field is very small, not more than i 
. ment at liquid-air temperatures, but failed com- 0.2 cm~', for example, in chrome or iron alum. E 
. pletely at helium temperatures, giving here a Hence the conservation of energy allows ex- P 
wrong order of magnitude for titanium alum, change of energy only between the spins and ( 
| and the wrong sign for the dependence of the oscillators of very low frequency w, which are far I 
| relaxation time of chrome alum on the strength from numerous inasmuch as the density of ‘ 
of the applied magnetic field. It is the purpose oscillators is proportional to w*? and which con- 1 
of the present article to point out that it is sequently cannot conduct energy rapidly to the ] 
! impossible for the lattice oscillators to preserve a_ walls. The same difficulty does not arise in con- Z 
if constant temperature in the face of frequent nection with the calculations at liquid-air tem- c 
energy transfers with the spin. Hence it is not peratures, as here interaction between spin and s 
| surprising that the calculations based on the lattice is secured by means of second-order s 
ordinary assumptions should partially fail and processes wherein a spin scatters a lattice quan- 7 
a | so, apart from basic thermodynamic features, tum of comparatively high energy with only a s 
| | the conventional theory must somehow be relatively small modulation of frequency, like the y 
| F modified. It was emphasized by Casimir himself? Raman mechanism in radiation theory. Such a 
/ " that at low temperatures the lattice, inclusive processes, however, have been shown by Waller* : 
“HEH Van Vleck, Phys. Rev. 57, 426, 1052 (1940). and others to be of minor importance at very Q 
(1988) B. G. Casimir and F. K. du Pré, Physica 5, 507 low temperatures. 
h 


3H. B. G. Casimir, Physica 6, 156 (1939). 


*T. Waller, Zeits. f. Physik 79, 370 (1932). 
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Width of the low frequency band.—If the spin 
of a paramagnetic atom had sharply defined 
energy levels, split only by a definite crystalline 
field and the applied magnetic field, then because 
of the conservation of energy the possible energy 
exchanges between spin and lattice would form a 
discrete spectrum, and only an_ infinitesimal 
fraction of the lattice oscillators would be on 
“speaking terms’’ with the spin. Actually, the 
energy-exchange spectrum will be blurred into a 
continuum for two reasons. In the first place, 
the spin-spin interaction will introduce a fluctu- 
ating effective magnetic field acting upon the 
atom, which is superposed upon the applied field, 
and which has a continuous distribution of 
magnitudes. Secondly, there is a diffuseness in 
energy levels due to the uncertainty principle, 
because the harmonic oscillations do not persist 
over a large number of periods. 

For purposes of estimating the order of magni- 
tude of the spreading in energy due to spin-spin 
interaction, it will suffice to assume that the 
latter broadens the spin levels into bands of 
width 2K8/h where 8 is the Bohr magneton 
number, and K is the root mean square spin-spin 
field. The frequency distribution law for the 
oscillators is Cw*dw, and their total number is 
3C(k6/h)® where C is a constant and k6/h is the 
usual Debye cut-off frequency. Hence the frac- 
tion x of the total number of oscillators which 
is in thermal contact with the spins becomes 
x~2BK- h*w?/4tk®@. Here hw, the center of the 
energy band of allowed exchanges, is of the order 
0.2 cm! in chrome alum, the precise value 
naturally depending on the strength of the 
applied magnetic field, while k@/h is approxi- 
mately 230 cm~', and we previously! estimated 
K=290 gauss, so that BK ~10~ cm. Thus only 
about one in 10!° of the total number of lattice 
oscillators is on “‘speaking terms’’ with the spin 
system. The observed relaxation time 7 for a 
spin is about 10-* sec. for chrome or iron alum. 
The number of exchanges or interruptions per 
second for an oscillator accessible to the spin is 
vi=(1/3xry) where y is the total number of 
atoms per paramagnetic atom (=48in KCr(SO«)s 
-12H,O). This interruption frequency »; is thus 
of the order sec.—. 

The estimates of the preceding paragraph, 
however, overlook the fact that according to the 


uncertainty principle, or its equivalent, one is 
justified in talking of the proper frequency w of a 
lattice vibration only if it persists uninterrupted 
over several periods. Our preceding estimate of »; 
is about 10 times as large as w. Consequently 
there will be a diffuseness in the oscillator levels. 
The corresponding blurring of spin states is 
negligible, as the spin-lattice energy exchanges 
are prorated among the spins of all paramagnetic 
atoms, much more numerous than the limited 
number of low frequency oscillators. This diffuse- 
ness makes the spins capable of exchanging 
energy with more oscillators than otherwise. 
The effect is roughly the same as though one 
permitted transitions in which conservation of 
energy based upon ideal monochromatic vibra- 
tions failed by amounts up to A»;. Consequently 
we should at first sight expect the width of the 
zone of oscillators accessible to the spin to be 
of the order 10'° sec.~' or 1 cm™", as this was our 
estimate of »; in the preceding paragraph. How- 
ever, we must allow for the fact that when the 
burden of spin exchanges is distributed over a 
wider band of oscillators, the interruption fre- 
quency per oscillator is diminished, and so we 
ought to use a smaller value of »; as soon as the 
band width becomes greater than that 10-? cm~! 
used to obtain our estimate 10" sec.~'. Con- 
sequently, when allowance is made for this 
compensating effect, the value of »; is changed to 
10° sec.-', and the correspoindng zone-width to 
0.1 cm. This estimate may seem to be some- 
what high because line breadth tends in general 
to be rather less than the interruption frequency, 
but we may regard this tendency as roughly 
counterbalanced by the fact that there are 
really several frequency bands, instead of one as 
we assume. Namely, in the alums there are four 
paramagnetic ions in the unit cell, and as a 
result there are four orientations of the applied 
magnetic field relative to the microscopic trigonal 
field, and hence four patterns for the decomposi- 
tion of the energy levels ; furthermore even one ion 
may have more than one frequency interval, as its 
ground state is usually split into more than two 
components even without spin-spin interaction 
or blurring due to the uncertainty principle. 

The mechanisms (a) and (b).—From the 
above discussion, it appears that the oscillators 
which can exchange energy with the spin system 
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are confined to a zone of width about 0.1 cm7. 
However, if no more oscillators than this are on 
“speaking terms’’ with the spin, they will be 
unable to transfer energy fast enough between 
the spin system and the helium bath to account 
for paramagnetic relaxation times as short as 
those observed experimentally. Fréhlich and 
Heitler> have shown that the heat conductivity 
of the spin system itself is negligible. Hence the 
heat flow between the spin system and the 
boundary must take place via the lattice oscil- 
lators. Two possible mechanisms for this flow are, 
(a) direct coupling between the low frequency 
oscillators and the helium bath and, (b) interaction 
of the low frequency oscillators with those of higher 
frequency, through anharmonic terms. In (b) there 
must be indirect contact between our low fre- 
quency oscillators and the bath via ordinary 
oscillators of higher frequency, for we have 
already noted that the heat capacity of the entire 
lattice system is inadequate for it to serve as a 
thermostat unless it is in equilibrium with the 
bath. We call vibrations of the “‘low frequency” 
type if their quanta are commensurate with the 
energy required to turn over a spin, and if they 
can thus exchange energy with the spin system. 
By the term “ordinary oscillators’’ we mean 
those whose frequencies do not appreciably ex- 
ceed kT /h, as vibrations with w>kT/h are not 
appreciably excited, because of the Boltzmann 
factor, and so need not be considered. A spin 
quantum is about one-tenth of k7, in the helium 
region, and so the low frequency oscillators are 
much less abundant than the ordinary ones. 
Inadequate conductivity of the low frequency 
oscillators.—To examine whether the mechanism 
(a) can be effective, that is to say, whether the 
temperature T, of the low frequency oscillators 
can be equilibrated by direct contact with the 
surrounding helium bath of temperature To, we 
imagine that we have a spherical specimen of 
radius R and suppose that the elastic waves are 
uninterrupted except when they interact with a 
spin or reach the walls, as this is the extreme 
case most favorable for the penetration of the 
bath temperature into the specimen. As already 
emphasized by Casimir in another connection,*® 


5H. Froéhlich and W. Heitler, Proc. Roy. Soc. A155, 
640 (1936). 
*H. B. G. Casimir, Physica 5, 495 (1938). 
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the mathematical analysis appropriate to this 
model corresponds essentially to radiation theory 
rather than to ordinary conduction. We may thus 
consider that we have a spherical hohlraum such 
that the elastic waves impinging on the wall have 
a temperature 7',, whereas those emitted by the 
walls have a temperature 7». If u(T) is the 
energy per unit volume of oscillators at tempera- 
ture 7, then the analogue of radiation intensity 
is I(T)=vu(T)/4r, where v is the velocity of 
sound (2.3105 cm/sec. in the alums). Unit wall 
area absorbs energy { {I(T cos@dQ=jvu(T) 
per unit time, and correspondingly emits jvu(7>). 
If is small, we may replace u(T,) —u(T») 
by.c(T1—To)f(w) Aw, where f(w)Aw is the number 
of low frequency oscillators per unit volume in 
the interval Aw, and c is the mean specific heat 
per oscillator. Hence the rate at which heat is 
transferred from the low frequency oscillators to 
the helium bath is 


dQ/dT =a(T1—T»), with a=rRvcf(w)Aw. (1) 


We can take c to have the classical value k = 1.37 
X10-'* ergs/deg. since our low frequency oscil- 
lators have quanta small compared to kT. We 
may suppose f(w) - Aw to be of the order 3.4 10", 
since the formula for the number of oscillators 
per unit volume in the frequency interval a, 
w+ Aw is 127v~*w*Aw and since our low frequency 
band involves an w and Aw about 6X10° and 
3X 10° sec.—!, respectively. The conductivity a is 
consequently 3.3 X 10-*R? watts/deg. 

Now according to the thermodynamical theory 
of Casimir and du Pré,? the rate at which heat 
flows from the spin system to the low frequency 
oscillators is 


dQ/dt=a(Ts—T 1), with a=(4rR*/3)cs/r. (2) 


Here Ts and ¢s are, respectively, the spin tem- 
perature and the spin specific heat per unit 
volume, and 7 is the paramagnetic relaxation 
time. The numerical value of cs is 48007~° 
ergs/cc, while 7 is the order 10~ in the helium 
region. Consequently at 1.4°K, the conductivity 
constant a@ involved in (2) is about 10'R’® 
watts/deg. Since the heat capacity of the low 
frequency oscillators is negligible compared to 
that of the spin system, the heat flows (1) and (2) 
must be equal, whence (7',— 7») /(T's—T 1) =@/a. 
Clearly, if these oscillators are to serve as a 
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thermostat, at substantially a temperature 7, 
we must have this ratio small compared to unity. 
With our estimates of a and a, however, a@ is 
thirty times as large as a for a specimen 1 cm in 
radius, and is three times as great even if R is 
as small as 0.1 cm. In other words the low fre- 
quency oscillators are interrupted much less 
frequently by collisions with the walls than by 
spin impacts, and so have a temperature much 
closer to that of the spin system than of the bath. 
The calculation of the effectiveness of process 
(b), ie., energy exchange between the low fre- 
quency and ordinary oscillators is more difficult, 
and is given in the following paper. We may 
assume that the ordinary oscillators are at the 
bath temperature 79, as this supposition gives 
the most favorable case, and is, as a matter of 
fact fully warranted, as the ordinary oscillators 
are good conductors because they are so nu- 
merous. The calculations of the following paper, 
especially its Eq. (21), show that then the rate 
at which heat is taken away from the low 
frequency oscillators due to the process (b) is 


dQ/dt=b(T,.—T»), (3) 


with b=3X10-°R® watts/deg. at 1.4°K if, as 
before, we assume that the band of low fre- 
quency oscillators has w=6X10°, Aw=3X10° 
sec.'. For a specimen 1 cm in radius, the process 
(b) is hence much less important than (a), and 
is nothing like rapid enough to make the low 
frequency oscillators serve as a thermostat. 
Since we have seen that neither process (a) nor 
(b) conducts heat fast enough to make the tem- 
perature of the low frequency oscillators simulate 
that of the container, one immediately wonders 
whether possibly the correct model is instead one 
in which the low frequency oscillators are more 
nearly at the spin temperature, and in which the 
“bottle neck’’ for the flow of heat between the 
spin system and the helium bath is located be- 
tween the low frequency oscillators and the 
container, rather than between the spin and 
these oscillators, as ordinarily supposed. Such a 
proposal, however, encounters serious difficulties, 
and is, in our opinion, not tenable. In the first 
place, the computed relaxation times become 
entirely too long. The values calculated under 
the conventional assumption that the low fre- 
quency oscillators serve as a thermostat were of 


about the right order of magnitude for chrome or 
iron alum, and already too great for titanium 
alum. With the mechanism (b), for instance, the 
formula for + becomes cs(a+6)/ab instead of 
Cs/a, where cs is the spin specific heat, and a, b 
are the conductivity constants for spin-lattice 
coupling and process (b), respectively. The modi- 
fied formula corresponds to the fact that the 
conductivity due to two resistances 1/a, 1/b in 
series is ab /(a+6). With mechanism (a), there is 
no unique formula for the relaxation time, as 
heat will be conducted away more rapidly from 
spins at the edges of the specimen, than those at 
the center. As a result the expressions for the 
dispersion and absorption become complicated 
functions of frequency, corresponding to a distri- 
bution of relaxation times + wherein r+ increases 
with the depth of penetration into the sample. 
Such a behavior is probably not in agreement 
with experiment, as the observed dispersion 
curves agree fairly well with those computed 
under the assumption of a single relaxation time. 
Another difficulty with (a) is that the mean 
relaxation time would increase with the volume of 
the specimen, whereas no dependence on size is 
found experimentally.? These complications due 
to spatial inhomogeneity do not arise in connec- 
tion with (b), as in (b) the ordinary oscillators 
serve as a thermostat throughout the volume. 
However, in (b) there is the further objection 
that the temperature dependence of r is wrong. 
Namely, the observed values of + are approxi- 
mately proportional to 7~*, whereas the theo- 
retical expression ¢s(a+)/ab is proportional to 
T~-* if a>b, inasmuch as* cs~T~, and as Eq. 
(21) of the following paper shows that )~7"“. 
We thus find that neither mechanism (a) nor 
(b) can account for the observed flow of heat 
between the spin system and the container, 
unless one assumes that somehow most of the 
thermal contact between the spin system and the 
lattice is made directly between spins and a 
substantial portion of the “‘ordinary oscillators,” 
whose frequencies are comparable with k7/h. In 
other words, energy exchange between spin and 
lattice takes place primarily in virtue of the 
handing back and forth of quanta which are 


7 W. J. de Haas and F. K. du Pré, Physica 5, 501 (1938). 
8H. B. G. Casimir, W. J. de Haas, and D. de Klerk, 
Physica 6, 241 (1939). 
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considerably larger than those of the order 0.2 
cm~! which represent the prevalent Stark split- 
tings. If the spin system makes thermal contact 
with a sufficiently large fraction of the ordinary 
oscillators, the conductivity of the portion of 
the lattice in communication with the spin will 
be sufficient to dispose of the heat liberated by 
the latter, and the ‘‘bottle-neck”’ will be between 
the spin and the lattice, rather than between the 
lattice and the bath. An idea of the fraction of 
oscillators required can be obtained by referring 
to Eq. (1). We require a conductivity 300, or 
more, times that previously computed on the 
basis of a zone of low frequency oscillators of 
width sec.—!, and mean frequency 6X 10° 
sec.~', as then for a specimen of radius 1 cm, the 
constant a in (1) will become about 10 times that 
a in (2). To obtain this increased conductivity, 
we must use in (1) a value of f(w)Aw about 300 
times larger than formerly. As the oscillator 
density f(w) is proportional to w’, the requisite 
increase will be obtained if we assume, for 
instance, that hw=2.0 cm, and that the band 
width hAw is at least of the order 0.3 cm, or 
alternatively, that hw and hAw are each about 
1.0 

Unfortunately it appears impossible to find any 
reasonable mechanism which will provide for 
the appreciable exchange of such abnormally 
large quanta between the spin system and 
lattice. Temperley'® has suggested that several 


® It may be objected that Eq. (1) based on the radiation 
model, furnishes only an upper limit to the conductivity, 
and so does not tell us the necessary number of oscillators, 
as this model does not take into account the possibility of 
the mean free path of an oscillator being terminated 
except by collisions with the walls. However, one can 
also solve the heat conduction problem for a sphere with 
surface at temperature 7» and containing a source of 
intensity tcs(Ts—To) per unit volume (cf. Eq. (2)). 
The necessary conductivity can then be gauged by the 
requirement that the temperature at the center of the 

here be much closer to 7) than to Ts. If one assumes 
that the conductivity per oscillator in the alums is about 
the same as that found experimentally by de Haas and 
Biermasz (Physica 2, 673 and 5, 619) for other crystals 
at helium temperatures, the estimate of the necessary 
number of oscillators thus obtained, though somewhat 
larger, does not prove to be greatly different from that 
supplied by the simple radiation model. We aim to de- 
termine only rough orders of magnitude, and so throughout 
the paper it is not worth while to distinguish between 
ra and w*Aw, or in the present connection to allow 
or the fact that the specific heat c in (1) does not quite 
ope the equipartition value k when hw is as high as 
0 

10C, N. V. Temperley, Proc. Camb. Phil. Soc. 35, 256 
(1939). 
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spins turn over simultaneously, but in our 
opinion,' the probability of such a process is 
negligible. Another possibility is that we have 
grossly underestimated the blurring due to the 
uncertainty principle. As we intimated in a 
previous paper,’ most of the difficulties would 
disappear if the line breadth due to the uncer- 
tainty principle were very large, of the order 
3X10" sec.—! (=1 cm"), about ten times as large 
as estimated in the first part of the present 
article. However, it is very difficult to see how 
this value could be so seriously in error, especially 
since estimates of line breadth based simply on 
the interruption frequency tend to be too high 
rather than too low. (For instance, the line 
breadth" of a harmonic oscillator due to radia- 
tion does not increase with the quantum number, 
although the amplitude and hence the transition 
probability does.) A blurring as wide as 1 cm=', 
would make the effective oscillator density at low 
frequencies practically independent of frequency, 
instead of being proportional to w as in the 
Debye law, and is doubtless contradicted by 
experimental evidence in other fields, such as 
specific heats. 

Effect of impurities.—A possibility which must 
be more seriously considered as a mechanism for 
the transfer of large quanta is the potential 
presence of impurities. By an impurity, we mean 
any paramagnetic ion whose splitting is very 
much larger than the mean value yielded by 
specific heat data, in other words is of the order 
2.0 cm— instead of 0.2 cm~!. The impurity need 
not necessarily be chemical, and instead can be 
an ordinary atom with an unusual crystalline 
(or perhaps even exchange) potential due to 
lattice imperfections, or in an abnormal valence 
state. The possible importance of impurities has 
beem emphasized to the writer by Professor 
Casimir, who points out that two samples not 
differing ostensibly in chemical and crystalline 
structure sometimes yield relaxation times differ- 
ing by a factor 2 or more. The first question which 
arises in this connection is whether it is possible 
to have impurities which are sufficiently rare as 
to be without appreciable influence on the 
specific heat, etc., but which are abnormally 
potent in spin-lattice coupling. This question can 


1 V. Weisskopf and E. Wigner, Zeits. f. Physik 63, 54; 
65, 18 (1930). 
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in a certain sense be answered in the affirmative. 
Namely, the specific heat is proportional to only 
the square of the Stark splitting v of the spin 
levels, whereas the spin-lattice coupling is under 
certain conditions” proportional to v®. Hence, at 
first sight, it seems that if, say, 1 ion in 10° has 
10 times the normal Stark splitting, the spin- 
lattice conductivity might be appreciably affected 
by such an impurity. However, this is not really 
the case, as the fact has been overlooked that 
such large spin quanta as 2.0 cm™~ cannot 
possibly achieve a common temperature with the 
more prevalent, tenfold smaller, quanta charac- 
teristic of the spins of the great majority of the 
ions. The energy exchanges among spins must 
satisfy the conservation of energy, and so spin 
temperatures can be equilibrated only by proc- 
esses in which the turning over of a large quan- 
tum is offset by the simultaneous turning over of 
about 10 small quanta, due to spin-spin inter- 
action. Such a cooperative or clustering phe- 
nomenon can be regarded as in many respects a 
generalization of the Temperley mechanism, and 
has a negligible probability (of the order e~*° 
or so). One way of escaping this difficulty is to 
assume that the impurities have a large over-all 
Stark splitting v, but that the splitting Av 
between certain sub-states, say a doublet, is 
comparable with the usual Stark intervals. Then 
the relative, but not the total population of the 
two doublet components conforms to the spin 
temperature. Because of this fact, the catalytic 
effect of the impurity is not as great as previously 
estimated, as it can be shown that an extra factor 
(Av/v)? is thereby introduced. Thus if the im- 
purities have about 10 times the usual splitting, 
our previous requirement of 1 in 10° would be 
changed to 1 in 10%. Such an abundance for 
impurities is presumably so high as to be out 
of the question, but quite irrespective of this fact 
it is highly unlikely that they would have the 
right kind of Stark pattern, wherein there are 
small intervals of just the right size super- 


2 Eq. (53) of reference 1 shows that the conductivity 
obviously contains a factor »* because the splittings 
Wy—W, in (53) are proportional to ». In addition 
there may be an extra factor »* because increased splittings 
are usually the result of closer upper states and so are 
usually accompanied by larger perturbing matrix elements 
H®(n';n'’). This can be seen, for instance, from the fact 
that Eq. (49) of reference 1 and Eq. (14) of J. Chem. Phys. 
7, 71 (1939) involve denominators of similar structure. 


posed on the large. The most likely possibility 
of this type would appear to be ions with an 
odd number of electrons having a degenerate 
or nearly degenerate orbital state. Then there 
would be a large total splitting, but small 
doublets due to the Kramers degeneracy which 
are decomposed only by the magnetic field. The 
only example of this situation in the iron group 
is, however, Co** if the crystalline field is of the 
usual nearly cubic type. Ions with an even num- 
ber of electrons might have certain sub-states 
which separate only in virtue of the applied 
magnetic field, or of small deviations of the 
crystalline field from the dominant symmetry, 
but the difficulty is that there would be no 
matrix components of the spin-spin interaction 
between the different sub-states, so that transi- 
tions between the latter would be impossible. 
For example, an ion with S=1 (e.g., V*+**) has 
no matrix elements of spin between the two 
degenerate levels Ms=-+1 in an axial field. 
Thus, all told, any impurity mechanism appears 
to be too highly artificial to be at all reasonable. 

Conclusion.._We seem to be between Scylla 
and Charybdis. We have shown that if the 
quanta exchanged between spin and lattice are 
comparable with the prevalent Stark splittings 
of the spin system, the latter is on “speaking 
terms” with only a very limited number of 
lattice oscillators, which are located at the ex- 
treme low frequency end of the spectrum, and 
which conduct heat to the walls or to the rest of 
the oscillators far too slowly to serve as a 
thermostat as presupposed in the ordinary ther- 
modynamic model. On the other hand, we have 
seen that any mechanism, such as impurities, 
which involves the transfer of larger quanta, 
and which so makes more oscillators accessible, 
is highly improbable. There is, however, one 
more possibility, which we believe is the way 
out, and which we will consider in a future 
paper. It is that the conduction of heat from 
the spin system to the helium bath takes place 
via virtual rather than real lattice states, i.e., via 
a sort of ghost wave. This process, though com- 
patible with Casimir’s thermodynamic postu- 
lates, seems at first sight rather strange, and so 
we have gone to some detail in the present 
article to exhaust the other, more conventional 
possibilities. 
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The rate of transfer of energy between lattice oscillators due to anharmonic perturbations 
is computed when different portions of the frequency spectrum are not in thermal equilibrium. 
This calculation is needed for the preceding paper and is not to be confused with that of 
ordinary thermal conductivity, since the present article is concerned with the flow of energy 
between oscillators at the same point of space, rather than that associated with spatial temper- 
ature gradients. The mean free paths for the two problems are shown to be unrelated, particu- 
larly at low temperatures, inasmuch as spatial thermal resistivity can be created by boundary 
or mosaic reflections, free from energy exchange, but does not arise from energy-transferring 
collisions between oscillators unless they are of the abnormal ‘‘umklapp” type in which mo- 
mentum is imparted to the grating framework. The force constants needed for the calculation 
for the alums are furnished by Bridgman’s compressibility data. The rate of flow of energy 
between a band of low frequency oscillators and the main body of lattice vibrations when 
the two are at slightly different temperatures proves to be proportional to 7‘ at low tempera- 
tures. The corresponding collision frequency between oscillators is about 10° sec.~' even at 
helium temperatures and thus is of such a magnitude that equilibrium between the different 
vibrations can be considered as secured instantaneously from the standpoint of macroscopic 


acoustical experiments, but not at all as far as paramagnetic dispersion is concerned. 


\ 


INTRODUCTION 


T is the purpose of the present paper to 
calculate the transfer of energy between 
lattice oscillators at different temperatures. We 
shall suppose that the oscillators at the low 
frequency end of the spectrum are at a different 
temperature than that of the rest and shall 
determine the rate at which the temperatures 
become equalized. It is to be emphasized that 
we are thus concerned with a different kind of 
conductivity than in the ordinary theory of heat 
conduction, for in the latter one studies the 
passage of heat between two spatial portions of 
the lattice at different temperatures, whereas we 
are interested in the flow of heat between 
different vibrations at the same position in space. 
Most physical problems involve spatial inhomo- 
geneity in the temperature distribution, but 
permit the assumption that all vibrations at the 
same point in space are at the same temperature. 
The reason for the latter supposition is that the 
transfer of energy between vibrations at the 
same point is more rapid than the spatial flow. 
Thus, for instance, in his interesting theory of 
the internal friction of solids, Zener! supposes 
that the rarefactions and condensations are at 


1C. Zener, Phys. Rev. 52, 230 (1937) ; 53, 90 (1938). 


different temperatures, but that the temperature 
distribution over the frequency spectrum at a 
given point in space is uniform, because, as he 
says, ‘‘thermal equilibrium between the different 
normal vibrations is established so rapidly.” It 
is our intent to calculate how great the latter 
rapidity really is, as its numerical value is needed 
in connection with our preceding article on 
paramagnetic relaxation, though it is so great 
that the equilibrium can be considered instan- 
taneous for a macroscopic elastic problem such as 
Zener’s. On the other hand, we are not concerned 
with spatial inhomogeneity. Thus our mathe- 
matical problem is quite different in nature from 
that involved in ordinary conduction theory and 
is in many ways a simpler one, for spatial 
inhomogeneity in temperature is achieved only 
as a sort of ‘‘wave packet”’ or interference effect 
between waves of different frequency. In the 
present calculation, we regard an oscillation of 
given frequency as having the same temperature 
throughout the crystal, and neglect the phase 
effects associated with the superposition of waves 
of different frequency. Thus with our model the 
spatial conductivity would be infinite. Never- 
theless our approximations appear legitimate 
because we are concerned with a much shorter 
relaxation time than in the spatial problem. We 
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can apply our results to problems in which the 
spatial distribution of temperatures is not uni- 
form, because the spatial interference between 
different waves is essentially a coarse-grained 
effect not affecting the total strength of the 
matrix elements connecting oscillators of differ- 
ent frequencies. 

At first sight, it perhaps appears that the 
calculations of the present paper are unnecessary, 
as one can determine a mean free path /x from 
the spatial thermal conductivity K by the usual 
relation K=43Clxv, in which C is the total 
specific heat per unit volume of all oscillators 
and v is the velocity of sound. If one knows the 
mean free path, it is easy enough to make an 
estimate of the rate of energy transfer, for the 
collision frequency is the ratio v/lx and the 
amount of energy transferred between colliding 
oscillators of temperatures 72,7, is, on the 
average, usually supposed to be of the order 
4c(T2—T), where c is some sort of a mean 
specific heat per oscillator. However, this way 
of estimating energy transfer is fallacious, as the 
mean free path entering in the spatial thermal 
conductivity is irrelevant for our purposes. 
Namely, as Casimir? has emphasized, the peculiar 
dependence of thermal conductivity on cross 
section found experimentally* in the helium 
region shows very clearly that the thermal 
resistivity at low temperatures is not due 
primarily to impacts between oscillators of 
different frequencies, but instead to reflection of 
the elastic waves at the walls, or at the surfaces 
of crystalline imperfections due to mosaic struc- 
ture. Such reflections can have no effect on the 
redistribution of temperature between oscillators 
of different frequency, as a reflection merely 
alters the direction of a wave without altering 


its frequency or energy. On the other hand, the 
most common variety of collision between 
oscillators is the type not involving any change 
in the combined momentum of the two inter- 
acting vibrations. Such collisions are of im- 
portance for the reapportionment of energy 
between oscillations of different frequency, but 
do not give rise to any resistivity for spatial heat 
flow. As Peierls‘ shows, only the abnormal 
“umklapp” processes, which impart momentum 
to the crystal grating, can influence the spatial 
conductivity, and at low temperatures such 
phenomena are unimportant compared with the 
wall and mosaic effects. Thus the collisions 
which are important for the purposes of the 
present paper are not revealed at all in measure- 
ments of thermal conductivity, and conversely, 
the latter involves reflections irrelevant for us. 
In other words, the mean free path or collision 
frequency for energy exchange between oscilla- 
tors has in principle no relation to that revealed 
by the thermal conductivity. Hence it is not 
surprising that the mean free path 1.0 cm which 
we calculate at 1.4°K is about ten times as 
large as that /x deduced in the usual way from 
measurements of the thermal (spatial) con- 
ductivity in a similar temperature region. The 
rate of flow of energy to or from a low frequency 
oscillator turns out to be less than one-hundredth 
that which would be obtained by using the 
naive formula }c(7,—T7>2)v/lx, as the latter errs 
not only in the type of mean free path, but also 
in supposing that the energy of two colliding 
oscillators is divided equally after impact. 
Actually a low frequency vibration will suffer 
only a relatively small modulation of its energy at 
collision, for it can only gain or lose a quantum hw 
which is small compared to its mean energy kT. 


DETAILS OF CALCULATION 


Force constants from Bridgman’s compressibility data.—We shall assume that all the atoms of the 
solid are equivalent, and that the interatomic forces are central in character. These approximations 
are necessary to make the calculations tractable, and presumably do not affect orders of magnitudes. 
Since the potential energy V;; between atoms 7 and j is thus supposed a function only of the distance 
ri; between them, the total energy of the solid takes the form 


2H. B. G. Casimir, Physica 5, 495 (1938). 


3 W. J. de Haas and Th. Biermasz, Physica 5, 619 (1938). 
*R. Peierls, Ann. d. Physik 3, 1055 (1929). 


ture 
ata 
s he 
rent 
itter 
| 
on 
reat 
tan- 
‘h as 
‘ned 
the- | 
rom | 
and 
itial 
fect 
| 
1 of 
ture 
lase 
Ves 
the 
ver- 
late 
rter | 
| 


732 J. H. VAN VLECK 


The derivatives of f;; necessary for our calculations can be evaluated from Bridgman’s compressibility 
data,’ which he writes in the form 


AV 1AV B 
Ves 


ie. p=—-——+— (2) 


Vo A Vo A? 


where Vo is the equilibrium volume at zero pressure. For potassium chrome alum, whose compres- 
sibility is no doubt sufficiently typical of the alums used in the Dutch relaxation experiments, 


Bridgman’s values of A and B are 
A=64.9X10-7, B=112.5X10-® (pin Kg/cm’). (3) 


Since p= —0U/0V and since dr;;/dV = }7;;/V we have by development of (1) in the displacements 
Ar;; from equilibrium, 


P= Disk GAV/ Vo) + fis” — V0)? +--+ J, (4) 


where r;;7f;;"’ means [7;;?(0°f/07;;") Jo, etc. In writing (4) we have omitted terms in f;,;’, as their effect 
vanishes in virtue of the fact that p=0 at = Vo, which is equivalent to the relation }-7;;f;,’=0. 
The sums 5>7;;*fi;/’ and ¥°r;;*fi;""" needed for later work can immediately be obtained by comparison 
of (4) with (2) and (3). 

Expansion of the potential to the third order.—-We now develop the energy (1) as a Taylor's series 
in the normal coordinates qi, ---, gsw for the crystalline system. Here N is the total number of 
atoms composing the crystal. It would be very difficult to find the spectrum of the normal coordinates 
exactly, and so we shall make the usual assumption that the frequencies are distributed in accord 
with the standard Debye model. There are, of course, no linear terms in the expansion, and the 
quadratic members, together with the momentum terms expressing the kinetic energy, merely 
determine the normal coordinates themselves. The third-order terms in the qg’s make the motion 
slightly anharmonic, and we regard them as a perturbation. They thus furnish a mechanism for 
the transfer of energy from one normal mode to another. The fourth-order terms we shall omit, 
as their effect in coupling together the different modes of vibration is qualitatively similar to that 
of the third-order members, but much smaller in numerical magnitude. 

We must now compute explicitly the terms of the third order in the expansion of the energy U. 
We write the displacement associated with a normal coordinate ga as Af;=Uaga exp(tk-r;) where 
r; is the equilibrium value of the position vector of a typical atom j, while k, is a vector of magnitude 
k,=22/dq directed along the direction of propagation of the wave associated with the vibration qa, 
and u, is a unit vector along or perpendicular to k, according as the wave is longitudinal or transverse. 
We shall now assume that the wave-length A, is large compared with 7;; (= |r;;|). This approximation 
is warranted at helium and hydrogen temperatures, where only low frequencies are excited, provided, 
as is doubtless true, the important interactions are between adjacent, or nearly adjacent atoms. 


If \./ris<1, we have 
Ar;— Arj= doa (5) 
The corresponding modulation in Ar;;= |r;;+Ar;—Ar;| —7;; is 
~ §P. W. Bridgman, Proc. Am. Acad. 64, 63 (1929). Strictly speaking, we want compressibility data at the absolute 


zero rather than room temperatures, but the variation of the elastic constants with temperature is a second-order effect 
of little importance for the present type of calculation. 
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The terms of the third order in the expansion of the energy (1) in terms of the g’s are thus 


The rather formidable expression inclosed in heavy brackets in (7) is a periodic function of 7 when 
summed? over 7, since all atoms are assumed identical. Hence, by the usual argument,’ the expression 
(7) vanishes unless 


+ks +k, (2x, d)m, (8) 


where m= (71), 2, m3), With m, m2, m3 integers, and where d is the lattice constant. Since the 
matrix elements of gq are of the form An,= +1, nonvanishing elements of (7) can arise only if each 
of the quanta changes by one unit. The energy associated with the vibration a, of frequeucy wa, 
is (n+4)hwa. From standard perturbation theory,’ it is known that an appreciable matrix element 
for energy transfer between different members of the system exists only if there is conservation 
of the total energy, i.e., 


=0, (9) 


where ka= |kKa| =22w,/v and where each of the An's can only be +1. It is possible to have a simul- 
taneous solution of (8) and (9) with m#0 in (8) only if two (or more) of the &’s are large enough 
to be comparable with 22/d, and involve An's of opposite sign in (9). Destruction of a large quantum 
at low temperatures, however, is virtually impossible, for the corresponding initial Boltzmann 
factor is negligible. Hence we can set m=0 in (8), so that (8) reduces to the familiar condition for 
the conservation of momentum. In other words we can neglect the so-called*’ ‘‘umklapp”’ processes 
corresponding to m#0. Clearly Eqs. (8) and (9), with m=0, are compatible only if the three vectors 
k., ks, k, are collinear, and if one of the k’s is the algebraic sum of the other two, and has the opposite 
sign for its An. In other words the destruction (or creation) of a quantum of one frequency is offset 
by the creation (or destruction) of two quanta of lower frequency. The transfer process is thus 
essentially a needle-like one, as all three quanta are directed along the same line. 

In the preceding paragraph, and elsewhere, we have assumed that the velocity of sound is the 
same for all frequencies, so that the wave number & is proportional to the frequency w. This as- 
sumption involves no serious error. If different waves have different velocities, the vectors k,, Ks, k, 
cease to be collinear, but our estimates of the magnitudes of transition probabilities are not much 
affected, as the order of magnitude of the density factor p(w.) in our later Eq. (11) is not greatly 
influenced by the non-collinearity. 

The fact that k., ks, k, are collinear greatly facilitates the summing in (7). We will suppose that 
the r;; can be regarded as directed at random. This assumption introduces no error of consequence, 
as a cubic crystal is very nearly isotropic. It is readily seen that then a nonvanishing result is obtained 
only either if all three of the vibrations a, 8, y are longitudinal, or else if two are transverse and 
one longitudinal. The averaging over the random orientations of r,;; is readily performed. To do 
this, we use the fact that the mean values of expressions of the type forms (r-e)‘, (r-e)®, 
and (r-e)*(r-e’)(r-e’’) are respectively r4/5, r°/7 and r°(e’-e’’)/35 if r is random, and e, e’, e” are 
unit vectors, with e’, e’’ both perpendicular to e. It is thus found that (7) becomes 

+ (6/210) (4 Lars Ley Us). (10) 

* The variable of summation i is not to be confused with the factor i= \/—1 preceding the wave number k in Eq, (7) 
and elsewhere. 


7Cf. N. F. Mott and H. Jones, The Theory of the Properties of Metals and Alloys (Oxford ay tw y Aig | p. 256. 
*Cf., for instance, P. A. M. Dirac, The Principles of Quantum Mechanics (Oxford University Press), edition, p. 160. 
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A superscript / or ¢ attached to a g means that it is a normal coordinate of longitudinal or transverse 
type, respectively. The factor 6 preceding the summations arises from ordering considerations, as in 
(10) the factors are arranged in a definite order, whereas this was unrestricted in (7). The inequality 
a>B> vy, etc., does not necessarily refer to the relative sizes of wa, ws, w, and is merely a device for 
excluding redundant terms. (The probability of repetitive cases such as a=8 is negligible.) 

We will now suppose that we are interested in a particular vibration of frequency we. We wish 
to compute the rate at which this oscillator changes its quantum states due to interaction with the 
other vibrations of the crystal. The latter will have frequencies mostly of the order kT/h, and will 
be supposed in thermal equilibrium. On the other hand, we will suppose that wo is at a different 
temperature from the main body of the lattice spectrum. In other words, we are dealing with a 
particular oscillator which is out of step thermally with the great bulk of vibrations of the crystal, 
and wish to determine how rapidly it acquires the prevalent temperature of the lattice. 

First let us consider the case that wo is a longitudal vibration. Then we, say, associate wo with y 
in (10) i.e. take ky=22wo/v. To determine the total probability A(m,—-n,—1) of a transition of 
the form n,—n,—1, we must integrate over all possible values of wa. If we suppose that n,—n.+1, 
then we may take mg—ng—1, so that we=wg+wo. Namely, if wo is small compared with kT/h (the 
case of greatest interest in the study of paramagnetic relaxation), the great bulk of collisions will 
involve values of wa, ws larger than wo>=w,, so that considerations of energy balance will prohibit 
Na, Mg from both decreasing or both increasing simultaneously. Even if w is comparable with kT/h, 
omission of the case An,g=Ang will not impair the orders of magnitudes of our results. By per- 


turbation theory,® we have 
kO/h 


A (n,n, —1) (| Us(mamgny; | *)w p(wa)dwa. (11) 
0 

Here the subscript Av means that for given we,ws we are to average over the different possible 
values of m,mg distributed in equilibrium fashion appropriate to the temperature 7. The weight 
factor p(wa)dwe is the density of states, relative to energy, for the we, ws system when w, takes on 
all possible values between we and wa+dwa. Eq. (11) involves an integration over wa because transi- 
tions satisfying the conservation of energy are possible for every value of wa. The upper limit of 
integration is the cut-off frequency associated with the Debye temperature @. 

Evaluation of o(@.).—The conservation of momentum is a rigorous consequence (apart from 
umklapp processes) of the periodicity of the lattice, so that the vector relation (8), with m=0, 
is exact. On the other hand, the conservation of energy must not be assumed to apply in computing 
p(w.), as calculation of the density in an energy space requires a knowledge of how states are dis- 
tributed over the energy interval in the vicinity of the point in question, although p(w.) is, of course, 
finally to be evaluated at the point at which the conservation relation (9) applies. The boundary 
conditions at the edge of the crystal require that Raz, Ray, Raz be integral multiples of 27/1, where 
l is a dimension of the crystal, supposed for simplicity a cube, with edges parallel to the coordinate 
axes. There is one longitudinal vibration, for each allowed value of k,, and furthermore, with given 
k,, there is one corresponding value of kg satisfying the conservation of momentum, viz. ks=Ka—Ky,. 
Hence the number of states in the interval kaz, Ray, Raz} Raz +dRaz, Rayt+dRay, Raz t+dRaz is 


dN =[V/(2m)* (12) 


where V is the volume /*. If we introduce elliptic coordinates \=}(ka+ks), u=3(ks—ka), ¢ having 
as foci two extremities of k, which for present purposes is regarded as a fixed quantity, then the 
right side of (12) becomes [2V/(27)*k, ](\°—u?)dAdudy. Clearly (hv/2mr)(ky+2u) is the lack of 
energy balance AW. If we use AW, wa=vk./2m as variables in place of \, wu respectively, we have 


4x? V 
dN=p'dAWdwde with p’=—-— 
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The quantity p(w.) needed in connection with (11) is }- 2p’. Here the factor 27 is contributed by 
the trivial integration over the azimuth ¢, while the factor } owes its origin to the fact that AW 
= (hv/2m)(|Ka—Ks| —ka+kg) is necessarily positive, and so the Dirichlet integrals, or equivalent, 
involved in the familiar derivation’ of (11) extend to only one side of the origin, instead of embracing 
both sides. Thus, since k, = 2rwo/v, etc., we have 


p(wa) = Vwgws/v woh). (a, 8 longitudinal) (14) 


This is the value of p(w.) to be used when the first part of (10) is substituted in (11). For the second 
part of (10), where a, 8 are transverse, we must use a value of p(wa) four times as great as (14), 
as there are then two transverse vibrations for each allowed value of ka or kg. 

Calculation of transition probability.—The matrix elements of the g’s, and the values of the n’s 
appropriate to thermal equilibrium are given by the usual formulas 


q(n; (n)w=1/Lexp(hw/kT)—1], (15) 


where M is the mass of the crystal. We can now substitute (10), (14), (15) in (11), and express the 
results in terms of the A and B constants, instead of the f;;, by comparison of (2) and (4). In squaring 
the second part of (10), one may replace (u,- ug)? by its mean value }. One thus finally obtains 


1) = nD I, (16) 
ith 
and 
Wa’ 2 hwe 
fn f f = 4 (18) 
o Lexp(hwe/kT) —1Lexp(hws/kT) —1] 0 


The expression (17) is independent of crystalline size since M/Vo is the density. The approximate 
value given in the second form of (18) applies to the case that hw» is small compared with kT, where 
wg = We — wo can be replaced by we without serious error. In the preceding paper we needed primarily 
this case, but the approximations used to simplify (18) furnish a crude estimate of orders of magni- 
tudes even when wa and wo are comparable. The replacement of the upper limit of integration by 
infinity is legitimate since we are interested in temperatures small compared to the Debye cut-off. 
Then the exponentials can also be simplified because the presence of the factor w,.’ws? makes most 
of the contributions to the integral come from values of hw/kT somewhat greater than unity. 

Equations (16, 17, 18) apply directly only in case wo is a longitudinal vibration. If it is transverse, 
we must identify wo with 8 rather than y in (10), and take w.=w,+wo. Correspondingly we replace 
a, B by a, y in (14), (18), etc., and in connection with the second part of (10) we must use a density 
twice rather than four times the expression (14). Thus when w» is transverse, the expression for the 
transition probability is an expression identical with (16, 17, 18), except that in (17) the terms 
with factors 54, 45 are deleted, in virtue of the absence of the first part of (10) for transverse waves, 
and the remainder of (17) is halved. 

If we use Bridgman’s values (3) of the constants A, B, and take v=2.3X 10° cm/sec., M/Vo=1.7 
g/cc, then the value of the factor DI for longitudinal oscillators wo is 2.010‘ sec.~'. The corre- 
sponding expression for transverse oscillators is 4.0 10* sec.—'. Thus the longitudinal low frequency 
oscillators are much more effective intermediaries than are the transverse ones in transferring heat 
between the spin system and the crystal proper, and so we will henceforth in numerical estimates 
consider only the longitudinal type. 

Rate of energy flow.—In the preceding paper, we needed the value of the conductivity for the 
heat flow between the low frequency oscillators and the main body of vibrations under the assumption 
that the former and latter are at temperatures T,, T respectively, with |7,—T|<T. If the low 
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frequency oscillators are located in the band wo, wo+Awo, and are Vf(wo)Awo in number this flow is 
dQ/dt =(A — 1) — A (mo — two Vf (wo) Aw, (19) 


inasmuch as an amount of energy Aw is added or removed from the vibration mp» at each collision, 
The two transition probabilities involved in (19) differ solely in their Boltzmann weighting, i.e., 


A (mony — 1)/A (no — no) =expl(hwo/ kT) — (hwo/kT 1) T)hwo/ kT*). (20) 


The mean value of mo needed in connection with averaging (16) is, by (15), very approximately 
kT /hwo since hwo/kT<1. By (16), (18), and (20), the expression (19) becomes 


dQ/dt=0(T1—T) with V f(wy) Aw. (21) 


The band of longitudinal low frequency oscillators in which we were interested in the previous 
paper had f(wo) =42w9"/v', w= 6X 10°, Aw=3 X10" sec.', V =4rR*/3 where R is the radius of the 
specimen. On combining these numerical values with the estimates of DJ obtained in the preceding 


_ paragraph, it is found that dQ/dT =3 X10-*R® watts/deg. at 7=1.4°K. 


Mean life.—Instead of computing the conductivity 0, it is in many ways more illuminating to 
examine the mean free time 7» of the oscillator w» between collisions, which is given by the formula 
to=1/(2A(nmo—no—1))w. Here the factor 2 appears in the denominator to allow for the fact that 
there are approximately equal probabilities of a low frequency oscillator absorbing and emitting a 
quantum. In our example, the numerical value of ro is 0.5X10~® sec. The corresponding mean 
free path 1=rov is 1.0 cm. 

The mean free time 7» is not to be confused with the relaxation time 7, required for the temperature 
difference between the low frequency oscillators and the main body of lattice vibrations to reach 
1/e of its original value, under the assumption that the specific heat of the former (viz., Rf(w)- Aw) 
is small compared to the latter. The formula for 7, is r;=kT'/hwDI and gives r,.=2X10~‘ instead 
of 0.5X10-* second at 7=1.4°K. The distinction between the two kinds of times arises because 
an oscillator loses or gains only a portion of its energy at each collision. 

Even a mean free time 0.5 X 10~* sec. is not nearly as short as required if the conventional theory 
of paramagnetic dispersion is to function properly at helium temperatures. In consequence, this 
theory must be abandoned, or very materially modified, as discussed in the preceding paper. By 
customary standards, however, the mean life ro or the relaxation time +, must be regarded as ex- 
ceedingly brief, at least at ordinary temperatures, and hence the process of temperature equalization 
practically instantaneous. Thus in experiments on acoustical vibrations, etc., all the oscillators can 
be regarded as at the same temperature at a given point of the crystal, though spatial inhomogeneities 
are of course possible. Our estimate ro>=0.5X10-° and r,=2X10~‘ sec. was for wo=6X 10° sec. 
and 7=1.4°K. In acoustical applications one is interested in much lower frequencies and much 
higher temperatures, but this will reduce 79 even further, since the integral J involved in (16) in- 
creases rapidly with 7, and since by (16) the value of r9>=1/(2A(n—>n —1)), is inversely proportional 
to the mean number of quanta mg. The frequencies involved in acoustical problems are so low that 
{no)» will amount to millions rather than being of the order 10 as in our example appropriate to 
paramagnetic relaxation in the helium domain. It will be noted that the interruption frequency due 
to anharmonic perturbations is much larger than the proper frequency when the latter is in the 
acoustic region. The resulting blurring, however, is insignificant, because acoustic vibrations have 
such high amplitudes that the loss or gain of a quantum at a collision represents only a minor 

disturbance. 
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A new method for treating the most important eigenvalue problems in quantum mechanics 
is developed. The solutions can be found immediately once the equations are factorized by 
means of linear operators. These operators acting on a normalized eigenfunction change it 
into a new normalized eigenfunction and all solutions can be found once the basic eigenfunction 
is known. This basic eigenfunction is a solution of a simple differential equation of the first 
order. The underlying theory is explained more fully on a special case (Section 1) and then 
the rules of procedure are formulated explicitly (Section 2). The rest of the paper contains 
applications including the Kepler problem treated according to Dirac’s theory. 


INTRODUCTION 


HIS paper contains a very simple method 

of solving the most important eigenvalue 
problems in quantum mechanics. The method 
applies however only to the discontinuous spec- 
trum. No power series or polynomial develop- 
ment is necessary. The difficulty of normalization 
is avoided; our formulae give the energy eigen- 


values and the normalized eigenfunctions. 


The basic ideas of the method are closely re- 
lated to those recently developed by Schrodinger." 
But there are also essential differences between 
these two methods. In Schrédinger’s language 
the basic difference can be expressed as follows. 
Schrédinger uses a finite number of infinite 
ladders, whereas I use an infinite number of 
finite ladders. All problems on which Schrédinger 
illustrated his method can be treated directly 
by mine. But the opposite is not true. This is 
evident if we compare the two treatments of 
the Kepler problem (in a Euclidean, or spherical 
space). In this paper the Kepler problem fur- 
nishes a direct application of the method. But 
it cannot be treated directly by Schrédinger’s 
method. What Schrédinger does, and what is 
usually done, is to use a transformation involving 
the coordinate and energy and to change the 
Kepler problem into a different, though related 
problem, accessible to his method. 

_ The first section of this paper contains a 
typical eigenvalue problem which I have tried to 
treat slowly and carefully. The discussion of all 
other cases follows the same pattern and I 


'E. Schrédinger, Proc. Roy. Irish Acad. A46, 9-16 
(May 1940). ; ° 


present them in a tempo which accelerates with 
their increasing number. 


1. THe First CASE 


When presented with a problem of the 
Sturm-Liouville type, we shall always change 
it into the form 


u=u(x). (1.1) 
This is possible if in the original form 
v=v(y), (1.2) 


the functions p, p are never negative and if p/p 
exists everywhere in the interval in which the 
solution is considered. The transformation lead- 
ing from (1.2) to (1.1) is? 


u=(pp)'v; v= f ( ) dy. (1.3) 
0 p 


In particular if p=p, the transformation from 
(1.2) to (1.1) does not involve a change in the 
independent variable.* The boundary condition 
will always be the same: the vanishing of u at 
the ends of the interval. In the case in which 
the interval reaches infinity, the condition that 
JS wdx (taken over the whole interval) is finite 
should be added. 


We shall now formulate our first case: 
m(m-+1) 
sin? x 


? Courant-Hilbert, Methoden der Mathematishen Physik 
(Springer, Berlin, 1931), Vol. 1, p. 250. 

®’ The choice of the form (1.1) is only to systematize 
ee The method could be applied just as well 
to (1.2). 


u+(A+1)u=0. (1.4) 
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Indeed the equation (1.4) has the general 
form (1.1); m is assumed to be a positive integer‘ 
(m=0, 1, 2---) and the interval for x is (0, x). 
Putting (A+1) in (1.4) instead of \ seems to be 
unnecessary but it is connected with the history 
of this equation and with its physical meaning 
about which I shall say a word more in Sections 
3 and 7. 

The next step is to factorize the equation (1.4). 
Once this is done, the problem of finding the 
eigenvalues for \ and the eigenfunctions is easy 
and can be solved almost without calculation. 
To factorize (1.4) means to replace it by the 
following two equations: 


{(m+1) cot x+d/dx} {(m+1) cot x—d/dx} 
Xu(A, m)={A—(m+1)?+1}u(a, m), (15a) 
{m cot x—d/dx} {m cot x+d/dx}u(d, m) 
={rA—m*+1}u(d, m). (1.5b) 
After having performed the indicated linear 
operations we find that both (1.5a) and (1.5b) 
lead to (1.4). Therefore (1.4) can be written 
either in the form (1.5a) or (1.5b). In (1.5) the 
dependence of ™ on X, m is expressed explicitly 
and that on x suppressed. The equation (1.5b) 
can be obtained from (1.5a) by interchanging 
the first two coefficients, and changing (m+1) 
into m everywhere with the exception of u(A, m). 
For later use we shall write (1.5) in the more 
general form: 


H™+0+ m) 
m) =(A—L(m))u(a, m), (1.6b) 
where, in our special case 
H™+=m cot x+d/dx; 
(1.7) 
L(m) =m?—1=(m—1)(m+1). 


Theorem I 
f = f (H™+¢)fdx, (1.8) 


if (gf) vanishes at the ends of the interval and 
H+ are given by (1.7). This means: the opera- 


* We omit throughout, the case of m<0 since this does 
not lead to new eigenfunctions. 


tors H- and H* are mutually adjoint. The proof 
is self-evident. 


Theorem II 


Let us assume that for a particular value of 
A=Ao and for m+0, u(Ao, m) is a solution of 
(1.4). (We do not care, for the moment, whether 
u(Ao, m) satisfies our boundary condition.) We 
can prove that 


u(Ao, m), —(1.9a) 
u(Ao, m—1) = H™*u (ro, m) (1.9b) 


are also solutions corresponding to the number 
pairs (Ao, m+1) and (Ao, m—1), respectively. This 
means: having a solution belonging to Ao, m, we 
can find, by applying our linear operators, new 
solutions corresponding to Ao, m+1 and Ao, m—1. 
Thus we can go one step higher and one step 
lower in the m’s. The argument can of course 
be prolonged and thus we obtain a Jadder of 
solution in the m’s, belonging to a definite Apo. 


For the proof we multiply (1.6a) by H@*)- | 


and (1.6b) by H™*. The result is 
(Hom+— 0+) (Qo, m)) 
m)), (1.10a) 
m)) 
m)). (1.10b) 


The comparison of (1.10) with (1.6) proves 
our theorem expressed in the equations (1.9). 
In our argument no use was made of the special 
form of the H*. 

We know that the ladder constructed in this 
way must end when the bottom, that is m=0, 
is reached since m is assumed to be 20. The 
essential point to which we are driving is that 
the ladder in which we are interested has not 
only a bottom but also a top, that is, it is a 
finite ladder. 

The equations (1.6a) and (1.6b) can be 
interpreted with the help of our ladder picture 
in the following way: by going one step up the 
ladder and one step down the ladder (or one 
step down and one step up) we arrive at the 
solution from which we started, but multiplied 
by (Ao—L(m+1)) or by (Ao—L(m)). 
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Theorem III 


If (Ao, m) is our solution, that is a regular 
solution vanishing at the ends of the interval, 
then the solutions obtained through (1.9) are 
also our solutions. The contrary is also true. 
Or: our and not our ladders do not mix. The 
ladder can be either entirely ours or entirely 
not ours. 

If u(Ao, m) is our solution (for m+0), then as 
a glance at (1.4) shows u(Ao, m) must go to zero 
at the end of the interval (0,7), at least as 
rapidly as sin? x. Therefore, H°*»~u(Xo, m) and 
H™*u(Xo, m) must be also our solution and thus 
the argument can be prolonged. This follows 
from the form of H”* expressed by (1.7). If, 
however, “(Ao, m) is not our solution then neither 
can be u(Ao, m+1) =H m), because go- 
ing one step up and one step down leads (except 
for a constant factor) to a solution from which 
we started. 

Though in proving this theorem we made use 
of the special form of our equation (1.4) the 
argument is much the same in most of the cases 
considered later. 


Theorem IV 
A necessary condition that our solution 
u(Ao, m), belonging to Ao, m should exist is that 
ho—L(m+1) 20. (1.11) 


For the proof we multiply (1.6a) by u(Ao, m) 
and integrate from 0 to 7. We obtain, because 
of Theorem I and (1.9a) 


m-+-1)dx 
f m)dx. (1.12) 


But this equation is not possible if \>—L(m+1) 
<0 and u(Ao, m) exists. In the important case 


ho —L(m+1)=0, (1.13) 
we deduce 
u(Ao, m+1) =0 (1.14) 


from (1.12) assuming that (Ao, m) is our solution 
(and, therefore, not identically equal to zero). 
Thus our ladder has its top; it ends for the value 
of m which satisfies \»9—L(m+1)=0, if such 


integer m exists. This theorem is again general 
and holds as long as (1.8) holds. 


Theorem V 


A necessary and sufficient condition for the 
existence of our ladder, belonging to (1.4) is: 


A=A,=L(l+1) =I(1+2); 1=0,1,2,+++. (1.15) 


a. The condition is sufficient—From now on 
let us write instead of u(A,,m), where 
1 (0, 1, 2 ---+) is defined by (1.15). For a chosen 
l there may be various m’s for which u,” is our 
solution. Let us choose 


m=l, (1.16) 


which is always possible, as m is an integer. 
The equation (1.6a) takes then the simple form 


(1.17) 

This equation is fulfilled if 
H+)~y,'=((1+1) cot x—d/dx)u;'=0, (1.18) 
= u;'=asin'*!'x; a@=constant. (1.19) 


This is our solution and as it contains an 
arbitrary constant it represents all the solutions 
belonging to m=/, \,=L(/+1) and satisfying 
our boundary condition. Indeed we see: If there 
were any other solution #,', regular and satisfying 
our boundary condition and for which H“+?~i,' 
+0, then we could go up the ladder to a,'*'+0 
contrary to Theorem IV. 

We cannot go higher towards greater m’s, 
because the higher solutions are zero solutions. 
But we can go down the ladder, obtaining always 
our solutions, always different from zero, as it 
can be seen from (1.12) and from the fact that 
L(l+1)—L(m+1) is always greater than zero 
for m=0, 1---1—1. Thus we have for a given / 
a finite ladder with the top corresponding to 
m=Il and the bottom corresponding to m=0. 
The solutions are: 

u,°, u;', u,', (1.20) 
defined explicitly by 
u,'=a sin 
(1.21) 


= H™*+u," =(m cot x+d/dx)u;,". 


b. The condition is necessary.—We shall prove 
now that these are the only solutions satisfying 
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our boundary condition. Let us assume that this 
is not so. First we see that our solutions do not 
exist for \o<0. Indeed in this case A\y»—L(m+1) 
=Ao—(m+1)?+1 is always negative (Theorem 
IV). The case L(/+1) = (/+1)?—1 for /=0, 1--- 
was discussed before. Let us now assume 
A=L(l' +1) = +2) where I’ is not 
an integer but greater than zero, and investigate 
the corresponding ladder which we assume to 
exist. If / and /+1 are two integers satisfying 


l<I’<Il+1, 


then u,' must exist as it was gained by a legiti- 
mate climbing up from the given solution u;/". 
But u,'*! cannot be our solution because L(l’+1) 
—L(l+2)<0. Thus from our solution we could 
obtain a solution by going up the ladder which 
is not ours, contrary to Theorem III. 

Summarizing we can say: To each \,;=/(/+2) 
(1=0, 1---) and only to such X’s, there belongs 
a finite ladder of solutions the top of which and 
the successive steps down to m=0 are described 
by (1.21). 

Although our argument applied to the special 
form (1.15) of Z(/+1) it will be much the same 
in all other cases. There is, however, one special 
feature of (1.4) which will not appear later. 
In the case of (1.4) we could equally well char- 
acterize the ladder by starting from the bottom. 
Indeed for m=0 the solution is simply 


u,)°=a sin (l+1)x 


and we could go along the same ladder up to the 
top instead of down to the bottom. 

There is one more problem: that of normaliza- 
tion. We shall now introduce the normalized 
functions U,;" and the normalized operators 3(- 
and X*. We define 

= (L(l+1) —L(m))H™+. 
(The new operators have two indices: / and m.) 


Um. (1.24) 


Theorem VI 


If U,™ is our normalized solution then U;"*', 
U;"" are also our normalized solutions. Indeed 
instead of (1.6) we can write: 


= U," 


(1.22) 


(1.23) 


(1.25a) 
(1.25b) 


INFELD 


and multiplying (1.25) by Uj," and integrating 
we have: 


=f (uryae=1. (1.26) 


By using the operators 3¢* instead of H* we 
thus obtain a normalized ladder if the solution 
from which we start, that is the top, is normal- 
ized. Thus the final, normalized solution of our 
problem can be written: 


A=l1(1+2); 
1 /2.4+++(21+2)\3 
Ur 
X(m cot x+d/dx)Uy", (1.28) 
l,m=0,1,2--+; ml. 


sin'*'x, (1.27) 


Although the reasoning leading to this result 
was fairly long, the mechanical calculations are 
practically non-existent, once the equation is 
factorized. Equation (1.28) and A,=L(/+1) 
=l](1+2) are known, once the factorization is 
performed and U;' is a normalized solution of 
an elementary differential equation of the first 
order. 


2. GENERALIZATION 


We shall now summarize in general terms the 
procedure which leads to a solution of all the 
eigenvalue problems to be considered. 

The equation of the type 


ul’ +r(x,m)utr\u=0; m=0,1--- (2.1) 
has to be brought into the form 
(k(x, m+1)+d/dx) (k(x, m+1) —d/dx) 
Xu(a, m)=(A—L(m+1))u(a, m), (2.2a) 
(k(x, m) —d/dx) (k(x, m)+d/dx)u(A, m) 
=(A—L(m))u(A, m).  (2.2b) 


We assume that this is possible and we assume 
further that L(/+1)—L(m+1)>0, for m=0--- 
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Fic. 1. 


Then the eigenvalues for are Xo, A1, Ae, 


where 
1=0,1,---. (2.3) 


To each value of \; there belong solutions 
u,°, u,', ° (2.4) 


The top of this ladder, that is u,' is determined 
by the equation 


(k(x, 1+-1) —d/dx)u,'=0, (2.5) 


from which follows 


u,'=C exp [ (2.6) 


where the constant C is determined by the 


condition 
1. 


The other expressions in (2.4) are given by the 
recurrence formula 


= (L(l+1) —L(m))-! 
x (R(x, m) +d/dx)u," (2.7) 


The eigenfunctions obtained in this way are 
normalized. 

Figure 1 will perhaps make the situation 
clearer. Once having the equations factorized 
we know X, as function of / and can therefore 


regard the solutions of (2.1) as functions of x, 


and two positive integers 1, m. The solutions 
along the line bisecting the (/,m) plane are 
given immediately by a simple quadrature. From 
each of these solutions a ladder leads down to 
the solutions corresponding to the same value 
of X; and obtained by successive application of 


the 35+ operation. No solutions lie ab®ve the 
bisecting line.® 

The questions which have so far been omitted 
and which are important for the application of 
this method are: 1.) When can an equation of 
the type (2.1) be factorized? 2.) How can the 
factors k and L be found if the factorization is 
possible? Here is the answer to these questions :° 
from (2.1) and (2.2) follows: 


m+1)+L(m+1) 
=—r(x,m), (2.8a) 

k*(x, m) —k’(x, m)+L(m) = —r(x, m). (2.8b) 
Replacing m+1 by m, in (2.8a) we have: 
k?(x, m) +k’ (x, m)+L(m) = —r(x,m—1), (2.9a) 
k?(x, m) —k’(x, m)+L(m) = —r(x, m). (2.9b) 
We introduce: 
—2a= —2a(x, m)=r(x, m)+r(x,m—1), (2.10a) 


28 =28(x, m)=r(x,m)—r(x,m—1).  (2.10b) 
Therefore (2.9) can be written 
k?+L=a, (2.11) 
k'=, (2.12) 
or 
2kk'=a'; k'=B (2.13) 
and therefore 
k=a' /28. (2.14) 


Thus the last equation gives us k, if the factor- 
ization is possible. Now the answer to our first 
question follows immediately. We introduce 
(2.14) into (2.11) and we have 


L(m) = a—a’?/4B?. (2.15) 


If L(m) calculated by (2.15) is a function of 
m alone, then the factorization is possible and 
(2.14) with (2.15) give all the data necessary for 
factorization, since a and 8 are known functions. 

Let us take as an example that considered in 
Section 1. We have: 


r(x, m) = —(m(m+1)/sin? x) +1; 
a=(m?/sin? x) —1; 
B=—m/sin? x. 
~ 8 This diagram was suggested to me by Professor Synge. 


*This problem was investigated along different lines 
by my students, Coleman and Lin. 
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Therefore : 
k(x,m)=mcotx; L(m)=m?—-1, 


exactly as before. This method can be used to 
obtain the factorization in all the cases treated 
later. 


3. SPHERICAL HARMONICS 


We shall now solve the differential equation 
for associated spherical harmonics: 


1 d dn m*n 
sin x— +An=0. (3.1) 
sin x dx dx sin? x 


We shall bring this equation in the form (1.1) 
by a transformation 


y = (sin! x)y. (3.2) 


Thus 7 are the ordinary associated spherical 
harmonics and y the density functions belonging 
to them. Equation (3.1) takes the form 


m2— 


(3.3) 


This equation can be factorized : 


{(m+4) cot x+d/dx} 
X {(m+4) cot x—d/dx}y 
={A—m(m+1)}y, (3.4a) 
{(m—}) cot x—d/dx} 
X {(m—}4) cot x+d/dx}y 
={A—(m—1)m}y. (3.4b) 
Therefore we have, in this case, comparing (3.4) 
with (1.6) 
L(m+1)=m(m+1); 


(3.5) 
=(m—}) cot x¥d/dx. 


The solutions for which 7 is regular in the 
interval (0, x) exist only if 


1=0,1,2,---. (3.6) 


The functions y,°, y?, y?-++y.' form a ladder 
the top of which is determined by 


(3.7) 


therefore, because of (3.5) 


y' =a (3 8 


determines our unique solution. By normalizing 
the top solution and introducing the 3’s instead 
of the H’s we obtain a normalized ladder of Y 
functions characterized finally by: 


A,=1(1+1); 
1.3---(2/4+1)\3 
vi'=( ’) sin'+}yx; 
2.2.4-++2/ -(3.9) 


X | (m—4) cot x+d/dx}y,".. 


We can formulate a more general problem 
characterized by the equation 


m?+2my—m 
etrAg=0, (3.10) 


+2¢'y cot x— 
sin? 


where m, as before, is an integer (m2 0) and y is 
an arbitrary positive parameter. Again this 
equation takes the form (1.1) by a substitution 


u=(sin’x)¢ (3.11) 
and we obtain 
(m+~y)(m+y-—1) 


sin? x 


u”’ 


(3.12) 


We see that this equation is a generalization 
of the two cases considered up to now. Indeed 
putting y=1, we have the equation (1.4) and 
vy =}3 leads to (3.3). The functions u are general- 
ized spherical harmonics and the condition that 
¢ are regular everywhere in the interval (0, x) 
leads to the determination of normalized func- 
tions U,,, satisfying (3.12). The result is: 


A=l(/+2y); 
U',4= 
(3.13) 
= 
X {(m+y) cot x¥d/dx}, 
U,, U", ,. J 
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The equations (1.27) and (1.28) follow from 
(3.13) for y=1, and the equations (3.9) for y=}. 


4. THE KEPLER PROBLEM 


The application of this method to the Kepler 
problem is straightforward and does not require 
the usual transformation which changes the 
original eigenvalue problem into a different one. 
The differential equation is: 


m(m-+1) 


2 


x x? 
where, in the obvious notation: 
A=(2h?/uZ*e*)E; (4.2) 


(The usual notation is to use / instead of our m 
in (4.1) and, later » instead of our /+1. But 
in order to make the comparison with Section 
1 simpler, we use here this unconventional 
notation.) 

We perform the transformation 


p=xyV(x) (4.3) 


and obtain in place of (4.1) the equation which 
has the form of (1.1): 


(4.4) 


The usual way to treat (4.4) is to perform 
the transformation 


&=(—A)!x. (4.5) 
Then (4.4) takes the form 


]p—p=0; 
(4.6) 


and ™ are functions of ~ But although we can 
find the solution of( 4.4) from (4.6) these two are 
quite different eigenvalue problems. This can 
be seen if we compare the orthogonality condi- 
tions for these two equations. If p:, p2 belong to 
different values of X and similarly f, je to 
different values of \, we have 


ff and ff (4.7) 
0 0 
Our method gives us, however, the direct 


solution of (4.4) and the use of the transforma- 
tion (4.5) is not necessary. 


The factorization of (4.4) gives: 
m+ 1 
x m+ 


a 1 d 
m+1 


=f A —}p, (4.8% 
(4.8a) 

m dx\ix m 


1 
m? 


which has again the general form of (2.2). 


Putting 
—1/(/+1)?, 


E,= — (4.9) 


we obtain regular solutious vanishing at the end 
of the interval (0, ©). To each d, of the form 


or 


(4.9) we have p:°, where 
pi (4.10) 
is a solution of 
(= 1 d 4.11) 
x dx 


From this top of the ladder we can reach the 
bottom by the recurrence formula: 


m 1 4d 
pi” '=Ht p= (~--+_) pi". (4.12) 
x m dx 


If we try to repeat the uniqueness proof of 
Section 1 we encounter one difficulty: for \>0 
the expression A—L(m+1)=A+1/(m+1)* is 
always positive. The uniqueness theorem can 
be proved only under the assumption \ <0, and 
only then does a discrete spectrum exist. 

The final result for normalized functions R,” 
can be expressed now : 


A= —1/(1+1)?; 
R'=(—) 
1417 
K (4.13) 
(I+ 


m 
x m dx J 


2 
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5. THe GENERALIZED KEPLER PROBLEM AND 
THE OSCILLATING ROTATOR 


We shall now generalize the Kepler problem, 
as before we generalized the equation for 
spherical harmonics. This generalization will 
later be applied to the solution of two problems: 
the oscillating rotator and the Kepler problem 
treated by Dirac’s equations. 

We start from the equation 


2 +y)(m+y—-1 
pry ty 


x 


f+aAf= 0, (5.1) 


assuming y>0O and looking for solutious which 
vanish at the ends of the interval (0, ©) and 
for which fo°f?dx exists. (For y=a=1 Eq. (5.1) 
goes over into (4.1).) 

The last equation can be factorized : 


a d 
x m+y dx 


m+y a d 
x| 
x m+y dx 


(5.2a) 
(m+)? 
a d 
x m+y—1 dx 
m+ 1 a 
x| 
x 


a’ 
= (5.2b) 
| (m++7-—1)? 


If \<0 we have a discrete eigenvalue problem 
for 


= (5.3) 
and the ladder f°, y, is determined by 


exp (—xa/(l+7)); 


-(—- +o) hes m+1 


= If), (5.4) 


INFELD 


Finally, the normalized solutions are 
Fh, y= (2a/y+)) 
Xx't7 exp —[xa/(/+y)]; 
2) (l—m) (5.5) 
m+y a d 


x m+y dx 


which again are a generalization of (4.13). 
The solution of (5.1) enables us so solve the 
problem of an oscillating rotator. Its equation is :’ 
2a = m(m+1)+a 


(5.6) 
x x* 


where a is a given, positive parameter. 
But, for a given m, the equations (5.6) and 
(5.1) are identical if 


y= (5.7) 
Therefore: ¢;” is a normalized solution of (5.6) if 
Fr, 


where y is defined by (5.7) and must be kept 
constant on the way down the ladder from 
F',, to F,,. Introducing (5.7) into (5.3) we 
have for \ the known expression : 


a 
(—A,)'= (5.8) 
n+3+[(m+})*+a]! 
with n=l/—m. This example shows that our 
method can sometimes be used even if a direct 
factorization is not possible. 


6. THE KEPLER PROBLEM TREATED BY 
Drrac’s EQUATIONS 


The Kepler problem treated by Dirac’s 
equations leads, for the radial functions, to the 
following differential equations 


(6.1a) 


dr r 


7A. Sommerfeld, Atombau und Spektrallinien (Vieweg 
Sohn, 1929), Vol. Il, p. 24-32. 

* Handbuch der Physik (Berlin, Springer, 1933), XXIV/1, 
p. 
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Here x1, x2 are the function densities, as the 
normalization condition is 


f (6.2) 
0 


x is an integer, positive or negative, but not 
equal to zero. E is the total energy and Eo=yuc° 
the rest energy. 

We introduce the following notation: 


yi=(k-—aZ)?; y2=(x+aZ)! 


) 


and replace x1, x2, by ¢1, ¢2 in (6.1) defined by: 


e=E Ey; b=mch 

(6.4) 
The result is: 

—bgi = (7, r+eby1/72) ¢2, (6.5a) 

(y/r—eby2/71) ¢1- (6.5b) 

¢2 by defined through: 

¥2= ¢1— ¢2, (6.6) 


Now, replacing ¢:, 


we obtain 


r y dr Y 


a=beaZ. (6.8) 


From (6.7) we easily obtain two equations of 
the second order in which the function y; and 
¥2 are separated. These are: 


y ad 

r y dr y dr 7’ (6.9a) 
y a@ dj{y ad 


(6.9b) 


But these are exactly the equations (5.2a) and 
(5.2b) for m=0 and m=1. Therefore: 


¥i=AF*,,; ¥2=BF',,, (6.10) 


where 


ex? a* 


The last equation, together with (6.8) and 
(5.3), gives the known formula: 


e= E/E,=(1+a°Z? (6.12) 


In (6.10) A and B are constants. The functions 
F’,, and F',, are detined by (5.5) and gained 
by descending from F',, to the lowest level 
(F°,.,) and to the level above (F';,,). 

But the ratio A/B is not arbitrary. We can 
find it by going back to the equations (6.7) of 
the first order. In terms of the H’s belonging 
to the generalized Kepler problem (Section 5) 
they can be written 


AH'-F*, ,=Bb(ex y—1)F'1,,,  (6.13a) 
F',,, =Ab(ex/y+1)F%,,.  (6.13b) 


Let us multiply these equations by the normal- 
ization factor N which changes the //’s into the 
K’s. Then we obtain from (6.13): 


AX '-F*,,, =AF',, 


= NBb(ex y—1)F'),,, (6.14a) 
BX F',,,=BF*,,, 
= NAD(ex/y+1)F%,,. (6.14b) 
N-? = (b?/y?) —1); (6.15) 


=[(ex—-y)/(ext+y) 
Thus the solution is 


The constant C is to be determined by the 
normalization condition (6.12). 

Is 1=0 permissible? For /=0, 2 vanishes, 
because the upper index is greater than the 
lower. From (6.12), and because «>0, we have 
= 7/ |x|. Thus ¥,=0 if and ¥,+0 if 
«<0. This result usually gained in a complicated 
way, follows immediately from our formulae. 

Going back to the functions x1, x2 from which 
we started, we have, because of (6.16), (6.6) 
and (6.4) 


(6.17b) 


— 
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7. THE KEPLER PROBLEM IN A HYPERSPHERE 


Schrédinger considered in his paper (quoted 
above) a very interesting case: the Kepler 
problem in a spherical space. This problem is of 
theoretical interest because its solution gives, 
unlike that in a Euclidean space, only a discrete 
spectrum. The equation which Schrédinger 
derives is: 


d da 
—(sin* +2y sin x cos xo 
‘x 
—m(m+1)o+dA sin? (7.1) 
h=(2u/h?)ER? (7.2) 


where 
v= (p/h?) RZe?*; 


and R is the radius of the hypersphere. 
A word about how (7.1) was obtained: The 
quadratic form of the space is assumed to be 


dr*+ R? sin? (r/R) (dd?+sin? dd 
= R*dx*+ R? sin? x(dd?+sin*? ddy*), (7.3) 
x=r/R. (7.4) 
After having written the Schrédinger equation 
Ay + (2m/h®)(E— V)y=0 (7.5) 
for this space, we introduce into (7.5) 
V = —(Ze*?/R) cot x, (7.6) 


because this V is a harmonic function in a 
spherical space, that is, it satisfies the equation: 


if 


(sint x— J= (7.7) 


The V in (7.6) corresponds for a small x, to the 
Coulomb energy —Ze?/r. To solve (7.5) we must 
express y as the product of the ordinary spherical 
harmonics and o(x) defined by (7.1), this equation 
replacing the old radial equation (3.1) for the 
Kepler problem in a Euclidean space. 

The transformation 


s=(sin x)o (7.8) 
brings (7.1) into the desired form: 


+1) 


m(m 
s’’+(2» cot x)s -————s+ (A+1)s=0. (7.9) 


sin? x 


Let us assume v=0. Evidently the term 
(2v cot x)s owes its existence to V, represented 


INFELD 


by (7.6). Therefore v=0 means a free particle 
in a spherical space. But then the equation 
(7.9) (for v=0) is identical with (1.4). Thus the 
physical interpretation of (1.4) is given. It is 
the radial part of the wave equation of a free 
particle in a spherical space. 

We can now turn towards solving, by our 
method, the more complicated equation (7.9). 
The factorization gives: 


v d 
{(m-+1) cot x— +—| 


m+1 dx 
d 
m+1 dx 
2 
(7.10a) 
(m+1)? 
v v 
cot m cot x -——-+— 
m m dx 
= (7.10b) 


A regular solution vanishing at the end of the 
interval (0, 7) exists, if 


A, =L(l+1) =1(1+2) —v?/(l+1)’, (7.11) 
L(l+-2)h? ety 


—Z? (7.12) 
2uR? 2h2(1-+1)2 


For small /’s and increasing R, (7.12) tends 
towards Bohr’s formula and for very large l’s 
towards the very dense spectrum of a free 
particle. The uniqueness can be proved similarly 
as in Section 1, and therefore the continuous 
spectra must be absent.® 

The top of the ladder s,°, s;'-++s;' is: 


s'=asin'*'x exp —(vx/1+1), (7.13) 


satisfying the differential equation 


d 
{ cot s'=0 (7.14) 
l+1 dx 


*For a given \, positive or negative, we can always 
find such a itive /’ that =L(/’+1). But our solutions 
exist only tL L(l’+1)—L(m+1) goes through zero for 
some integer m. Therefore /’ must be an integer if our 
solutions exist. 
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and eigenfunctions of all other levels are given 
by: 


v ad 
(m cot (7.15) 
dx 


8. OTHER EXAMPLES 


In our treatment of the most fundamental 
equations in quantum mechanics, one was left 
out, that of a harmonic oscillator. It seems not 
to fit into the developed pattern, because of its 
equation 


—x*yt+Ay=0, (8.1) 


the parameter m does not appear explicitly. 
I shall show briefly how to treat this case though 
the recurrence formulae to which it leads are 
well known. Instead of (8.1) we write: 


and the factorization gives: 
(x—d, dx)(x+d (8.2a) 
—(8.2b) 
therefore 
V=2+1; 
yi" =[2(l—m+1) 3(x—d 


mt exp (—x?/2); 
(8.3) 


or putting »=/—m we have the normalized 
eigenfunctions for \=2n+1: 
=m ‘exp (—x* 2); 
+1) dx), (8.4) 
(2n)-?(x+d, dx)y,. 
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From (8.4) we have 
(2(n+1)) Wari t(2n) (8.5) 


the application of which allows us quickly to 
calculate the intensities. Indeed we have 


f dda (8.6) 


and therefore 


(8.7) 


A few words now about Bessel’s equation: 
It can easily be factorized, but, of course, it 
cannot be treated in the described way. Indeed 


Z" —[(m?—}) x*]s+As= (8.8) 
can be written: 
| (8.9a) 
x 
m—} m— d 
=z. (8.9b) 
x 


The described method breaks down, however, 
because L(m) = L(m+1)=0.!° 

The fact that all the important quantum- 
mechanical eigenvalue problems can be treated 
by this method seems to suggest that it is 
something more than merely a mathematical 
trick. Unfortunately, I was unable to find a 
deeper reason for this. ° 

My thanks are due to Dr. Griffith, Professors 
Synge and Stevenson for reading the manuscript 
and for helpful criticism. 


°It may be noted in passing that we can deduce 
quickly the known recurrence formulae from (8.9) : 
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When a sinusoidal voltage is applied to a condenser containing an absorptive dielectric 
the absorption current can be split up into two components, one in phase and another in 
quadrature with the voltage. In the present paper, direct relations, suitable for numerical 
calculation, are established between these components, and a general relationship between 
the dielectric constant and the dielectric loss is given. 


F one admits the validity of the Hopkinson 

principle of superposition, the current J(t) 
in an absorptive condenser of geometric capaci- 
tance Cy and leakage resistance R=1/G» is 
given by! 


t)dr, (1) 


dU t dU(r) 
I(t) f 
dt 


where U(t) is the applied voltage and ¢(¢) the 
dielectric relaxation function (¢(*)=0). 
When the voltage is sinusoidal, 


U(t) = (2) 


it is easily shown, that the expressions for steady 
state current and for dielectric loss can be 
written formally’ in the same way as in the 
case of a non-absorptive condenser by intro- 
ducing a conductance function G(w) and a 
capacitance function C(w). 

The steady state current is given by 


I(w) = [G(w) + jwC(w) JU(w) ; (3) 


the dielectric loss factor is given by the tangent 
of the loss angle, 6 


tg 6=G(w)/wC(w), (4) 


where 
G(w) f Ge 
and 


C(w) = Cot f Ch 
0 


1J. B. Whitehead, Dielectric Theory and Insulation 
(McGraw-Hill, 1928); M. F. Manning and M. E. Bell, Rev. 
Mod. Phys. 12, 215 (1940). 


Let us put 
A (w) =G(w) —Go, (6a) 


B(w) = C(w) — Co. (6b) 


The functions A(w) and B(w) characterize the 
absorption current. A U is the part of this current 
in phase, and BU the part in quadrature with 
the applied voltage. 

Within the whole range of frequencies, up to 
about 10" cycles per second, concerned here A (w) 
and B(w) are often obtained by direct measure- 
ments. But it is impossible to measure directly 
¢(c) within the corresponding range of values of 
its arguments because this involves d.c. measure- 
ments for extremely short times.* The relations 
(5) are, therefore, of no use for experimental 
work at high frequencies. An examination of 
these formulae shows, however, that they can be 
inverted and that ¢(c) can even be eliminated 
from them. The purpose of this paper is to 
establish a direct relation between the functions 
A and B. 

The validity of this relation is not restricted 
to the theory of anomalous dielectrics. It also 
represents the expression of a new theorem in 
the theory of linear networks. And we believe it 
to be of interest even from a more general point 
of view, as it will prove to be an application of a 
very general integral transformation mentioned 
until now only in purely mathematical papers. 

The theory of linear networks starts formally 
from the same equation as does the theory of 
dielectric after-effects, i.e. from Eq. (1), if written 
in a slightly different form by introducing the 


2? Most measurements extend only to 10% or 107 sec. 
Only in the case of highly anomalous dielectrics does one 
succeed in extending the range of measurements up to 
10-* sec. (B. Gross, Zeits. f. Physik 108, 598 (1938)). 
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indicial admittance h(c) and omitting the term 
C,(dU/dt). But by separating the d.c. component, 
h(c) can be written always in the form f(c) 
+const where f(«)=0; and the theory remains 
more general retaining a term Co(dU/dt), apply- 
ing then to impulsive networks too. The main 
difference between the two theories consists in 
the shape of and h(a), respectively; does 
not change its sign and is continuously decreasing, 
whereas h(c¢) may, or may not be oscillating, but 
for dissipative networks every oscillating term 
contains a damping factor e~”’. In both cases, 
therefore, the integrals over | ¢(c)| and |h(c) 
—const|, respectively, exist from 0 to ©. Hence 


ff (7) 
0 


Q is the total residual charge of an anomalous 
condenser after application of a unit voltage 
during infinite time. A corresponding equation 
is valid for |h(c)—const}. 

All calculations based on Eq. (1), so far as 
they do not use any particular analytical ex- 
pression for g(c), but only the general property 
(7), are therefore valid both for anomalous 
dielectrics and linear dissipative networks. Only 
the names have to be changed. For the theory of 
networks we search for an expression which gives 
a relation between the real and the imaginary 
components of the steady state admittance. 

g(c) is certainly continuous for 0<e=~. 
Equation (7) shows that the integral over | ¢(c) | 
exists from 0 to «. This behavior enables us to 
consider (5) as Fourier integrals which can be 
inverted by Fourier’s transform, thus giving 


2 
=— f ) sin wodw, (8a) 
fa ) dw (8b) 
(w) cos wodw. 


The knowledge of any one of the functions A (w) 
or B(w) is sufficient for the calculation of ¢(c). 
The relations (8) suggest the possibility of deter- 
mining the dielectric relaxation function for 
values ¢X<1 by means of a.c. measurements. 
This method would not suffer from the experi- 
mental difficulties inherent in d.c. methods. It 


may furnish this function over a much wider 
range of o than known hitherto, but it must be 
said that the problem of evaluating the integrals 
(8) for numerically given functions A(w) or B(w) 
is difficult. 

The relations (8) need a further comment. 
The function ¢(¢) is equal to zero for all values 
¢<0. This is not the behavior of the functions 
defined by the integrals (8), which are, respec- 
tively, an odd and an even function. These 
integrals represent ¢(o) only for positive values 
of o, but this does not limit the validity of our 
calculations, because the particular form of Eq. 
(1) excludes the occurrence of negative argu- 
ments in ¢(c). 

In order to eliminate ¢(c), we substitute it in 
(5a) by its value (8b) and in (5b) by the ex- 
pression (8a). There results 


A(w) 2 
2 f sin | B(a) coscada, (9a) 
0 0 


y A(a 
B(w) =— f cos woda f sin cada. (9b) 
Qa 


These relations would be of no use for numerical 
applications if it were not possible to simplify 
them. Fortunately it is possible to do so. One 
obtains finally (see Appendix) 


The integrals are principal values. 

These are the desired relations between anoma- 
lous capacitance and conductance. Equation 
(10a) gives G(w) —Gp as a function of Cfw)— Co; 
Eq. (10b) gives C(w)—Cy as a function of 
G(w) —Go. In this form, the integrals are suitable 
for numerical application. In consequence of the 
factor 1/(a?—w*), these are mainly the values of 
the integrand in the vicinity of the point a=o, 
which contribute to the value of the integrals for 
a given w. The graphical or numerical evaluation 
of the integrals is therefore not difficult. 

The expression for the dielectric loss is now 
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obtained easily. Taking into account that 


f (11) 
0 


and observing that for technical, and higher, 
frequencies 


there results 
2 w * C(a) 
tg —— ——da. (13) 
C(w) w?—a? 


The knowledge of the capacitance function C(w) 
is therefore sufficient for calculating the di- 
electric loss. 

The form of the relations (10) is remarkable. 
These equations are a pair of integral equations 
of the first kind, each one of them inverting the 
other one. It can be shown that they apply to 
very general periodic or nonperiodic functions, 
which can be represented by these integrals. 
For periodic fynctions, the integrals may be 
simplified and result finally in Hilbert’s inversion 
formula.* But the investigation of the general 
properties of the system (10) seems to be a 
matter of purely mathematical interest and will 
be given elsewhere.‘ 


APPENDIX 


The relations (10) can be obtained directly from (9). 
Here we prefer to verify them by substitution. 
If we substitute in Eq. (10a) the value of B(a) given by 


Eq. (5b), that is 
B(a) = cos aede (14) 


and write the result explicitly, we have 


Aw) 


| e(e) cos (16) 


*D. Hilbert, Grundzuege einer <n Theorie der 
Integralgleichungen (Lei ig, 1912), p 
B. Gross, Ann. Aca ras. Sci. 13, (1941). 


But 


The integral (16) is absolutely and uniformly convergent. 
Now consider the second of the repeated integrals in (15), 
There exists the corresponding double integral ZL 


L=f" 


because the functions under the integral are all continuous 
and because 


COs (17a) 


The order of integration in Eq. (15) can therefore be 
inverted: 


2,. 


0 w?—a? wb w? 


Designating by F the expression in parentheses, one sees at 
once that it can be written in the form 


cos 
a 
+c0s au] dat f° 


cos 


a 
=sin ow[ —2Si (19) 
where Si is the sine integral. In consequence of the absolute 
convergence of the integral /(°¢(c) sin wede and because 
|Si (06) | <const and <const, for « ,0=8=w 
we have® 


Stim (e)|2 sin Si (08) + one | 


de=0, (20) 


and therefore finally 
A(w) sin wade, (21) 


in accordance with Eq. (5a). In a similar manner, one 
verifies (10b). 


5R. Courant, Differential-und Integralrechnung II 
(Berlin, 1931), second edition, p. 251, 252. 

6 E. W. Hobson, Theory of Functions of a Real Variable 
(Cambridge, 1907), first edition, p. 597. 
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Theory of Degenerate Non-Ideal Gases 


L. I. Scuirr 
Randal Morgan Laboratory of Physics, University of Pennsylvania, Philadelphia, Pennsylvania* 
(Received January 27, 1941) 


The partition function for a non-ideal quantum gas is calculated by considering the de- 
generacy exactly and the inter-particle forces only so far as binary collisions are concerned. 
In the absence of forces, this reduces to the ideal Einstein-Bose or Fermi-Dirac gas, as the 
case may be, rather than to the ideal classical gas. If the wave functions describing the two- 
particle collision are known, the partition function can be reduced to a set of quadratures. 
The theory is applied to the phase transition in liquid helium in the following paper, and 
other possible applications are discussed here. 


1. INTRODUCTION 


HIS paper develops the theory of non-ideal quantum gases in an approximation in which the 
degeneracy is taken into account exactly and the inter-particle forces are considered to an 
accuracy corresponding to the classical second virial coefficient. There would seem to be ample 
opportunity for application of this theory, especially to those non-gaseous systems where the ideal 
quantum gas provides a reasonable first approximation. Although the existing treatments of gases at 
low temperatures! consider the quantum statistics only so far as it affects the collision of a pair of 
particles of the gas, it is probably not worth the effort to regard the statistics as operative on all 
the particles unless the forces can also be considered exactly. For the contribution of the forces to the 
second virial coefficient of helium gas is greater than that of the statistics, at all temperatures.' 
However, there are some non-gaseous systems: the electrons in a metal, the neutrons and protons in a 
heavy nucleus, and the atoms in liquid helium,’ of which at least a qualitative understanding can be 
obtained in terms of a quantum ideal gas. It would be of interest to see what effect would be produced 
by the introduction of binary collisions if the exact treatment of the degeneracy were not sacrificed. 
The procedure followed here consists in calculating the partition function for the gas, by using 
completely symmetrized wave functions and considering those regions of configuration space where 
not more than two particles are within their interaction distance of each other (although many 
pairs may be close together at once). Thus in the limit of no interaction, the ideal Einstein-Bose 
or Fermi-Dirac gas appears. As we shall see, this procedure is feasible only if the particles of the 
gas are, apart from their mutual interactions, in a field-free region. A complication arises when one 
or more particle states of the ideal gas of the same temperature and density each contain on the 
average a number of particles large compared to one. The complication is not a serious one; and 
since it probably appears only in the theory of liquid helium, it is discussed in that connection.® 


2. THE PARTITION FUNCTIONS 


All of the thermodynamic properties of a system can be derived from the partition function or 
sum-over-states* 


Z=), exp—(8E,), (1) 


.where 8=1/kT and the index yu labels the quantum states of the whole system. Since we shall wish 


* This research was begun at the University of California, Berkeley; the writer wishes to thank Professor J. R. Oppen- 
heimer for discussion during its early stages. 
1G. E. Uhlenbeck and E. Beth, Physica 3, 729 sim f 4, 915 (1937). J. de Boer and A. Michels, Physica 6, 409 (1939). 
he My pe Phys. Rev. 50, 963 (1936) ; 55, 1095 (1939). 
ondon, Phys. Rev. 54, 947 (1938). 
%%.. Schiff, Phys. Rev. 59, 758 (1941); following paper. 
*R. C. Tolman, The Principles of Statistical Mechanics (Oxford, 1938), pp. 533, 567. 
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to separate out those regions of configuration space where various numbers of particles are close 
together, we write this® 


where the ¥, are any complete orthogonal set of properly symmetrized wave functions for the N 
identical particles, normalized in a (rectangular) box of volume 2. We assume that the Hamiltonian 


can be written 


Vin (3) 


i=1 >i 


where 7; is the kinetic energy operator of the ith particle, and the potential energy V,; depends 
only on the distance between the points x; and x;. This of course implies that there is no field within 
the box. From (1) and (2) it is clear that the operator e~*” is defined by the power series 1—8H 
+3@H?—---. It is convenient to use wave functions that obey periodic boundary conditions at 
the edges of the box; for spinless particles that obey Einstein-Bose statistics these are 


P exp[i k;-x;], (4) 


where A is a suitable normalizing factor and }p represents a sum over all different permutations P 
of the coordinates x;. The partition function then becomes 


f Pexp[-i k;- x; -exp [i k;-x; -+dXy, (S) 


where each state yu of the whole system is specified by a set of non-negative integers m, that designate 
the numbers of vectors k; that are equal to k, subject to the restriction }>; m.=N. 

We now suppose that the potential V;; is of short range, so that two particles interact appreciably 
only if one lies within a small volume w centered about the other.® The first term in Z, in which 
no particles are within w of each other, is obtained by replacing H in (5) by > 7;; this gives 


Zo=d, expl Liz mex], (6) 


where €,=h?k?/2m and m is the mass of the particle. The next term, in which only one pair of 
particles are within w of each other, is obtained by replacing e*” in (5) by Dis; fexpl—A8(L: 


T.+ Vii) T.]}: 
Z:=Zo 3ni(Mm— dim) f > P» exp[ ] 
P2 


fexp[ —8(71+T2+ Viz) ]—exp[ —8(7:+T2) J} Jéxidx». (7) 


Here, >> p2 means a sum over different permutations of the two coordinates; the integrals may be 
taken over all space, since the integrand practically vanishes outside the region w for the relative 
coordinates. If two relative coordinates are within w we obtain a two-pair contribution to Z by 


5 Cf. J. C. Slater, Phys. Rev. 38, 237 (1931). 
* This follows the treatment of H. D. Ursell, Proc. Cam. Phil. Soc. 23, 685 (1927); cf. R. H. Fowler, Statistical Mechanics, 
second edition (Cambridge, 1936), p. 241. We neglect the contributions from ternary and higher order collisions. 
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replacing e~®” in (5) by. 
* i(expl Vist Ver) ]—expl—8 — (exp T.+Vi,)] 
~exp[ - BY T+ Viv) ]—exp[— 8 
= A(T, +T;+ Vi) J]- 8(T,+T,;) })(exp[ aT +T. 
expl-6 T,])}. 


An induction is readily established, and one sees that the r-pair term is 


where 


Per 
X (expl —BHi2]—exp[ — BH 12° }) - - -(expl 2-1] 
—exp[ —BH"2,, 2,1 ]) -X2,) --dX2,, (9) 


and we have abbreviated 


Vig= Hijo + Vij. (10) 
Thus to the accuracy of the classical second virial coefficient, the partition function is given by 
iN 
Z=> 


r=0 
The sum over yz in (8) is readily carried out if we replace the canonical ensemble implied by (1) 
by a grand canonical ensemble,’ for which the total numbers of particles in the member systems 
are no longer the same. A new partition function may be defined: 


W= =>: W,, (11) 
N’=0 r=0 
where 
W.= (12) 
N’=0 


The thermodynamic quantities can be obtained in terms of W as follows.* The probability that a 
member of the ensemble has N’ particles and is in the energy state E, is exp(8¥—aN’—BE,). Since 
this probability is normalized to one, we have that BY= —InW. The average energy and number 


of particles are given by 
N=—(dInW/da)a,», (13) 


and the latter relation serves to fix the value of a. Suppose now we consider variations of the pa- 
rameters of the ensemble (Q, a, 8) subject to the restriction 6N =0. Putting 


— S/k=pY—aN —BE, (14) 


we see on making use of the normalization condition 8(é)+p62)+(y—E)s8—Néa=0, that 


Té6S=6E+p6Q. Here p=(—@E/dQ) is the pressure, and so S is the entropy. The Helmholtz free 
energy is then 
A=E-TS=y-aN/8, (15) 


7R. C. Tolman, reference 4, p. 510; the method of Darwin and Fowler can be used equally well. 
® This follows Tolman’ s treatment for Z (reference 4). 
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and the pressure can be put into a form more convenient for calculation: _ 
p= (0A /dQ)s,x =(1/8)(9 INW/dQ),, 4. (16) 


The specific heat at constant volume is a 
x. (17) 
With the help of the relation 


we obtain from (8) and (12): 


where the vectors 1,-++l:, are alike in groups of s; (}0; s;=27), and the »; are defined by 
(exp[a+ Be, ]— 1)-*. (19) 


The quantities »; are just the average numbers of particles in the various free-particle states in the 
absence of forces. Now although we cannot assign individual particles to definite states when they 
can interact with each other, we shall see that the parameter a is not changed appreciably by the 
introduction of binary collisions, and so we can put 


(20) 


Equation (20) will be used only to establish the relative orders of magnitude of the various terms 
that we encounter. 

The evaluation of the sum in (18) is complicated by two factors: the term II; s;!, and the large 
number of exchange terms that appear in G, as given by (9). Now G,(li- --1:,) is a relatively slowly 
varying function of its arguments, and for low temperatures the »; are appreciable only for small 1. 
Thus we can estimate orders of magnitude by replacing G, by some average value G, in (18). By 
making use of (20), the sum becomes simply G,N* if the term II; s;! is left out; the error that arises 
from omission of an s! when s of the 1; are the same is of order >; s!v,*/N*~s!(b/N)*', where i 
is an average value for the »;. Since we have assumed that none of the », is large compared to one, 
this error will be negligible for 1<s<«WN. It is conceivable that an appreciable error could be made 
by omitting the s! terms when s is of order N, and we return to this point below. 

Assuming that the omission is justified, the sum in (18) is seen to be composed of a large number 
of exchange terms of various kinds, one of which is yo", where 


y= v.viGi(k, 1) ; (21) 


this arises from that part of G, in which at most x; and Xe, x3 and x4, etc., are exchanged on the 
left side of the integrand of (9). The next most simple type of term is one in which the pair X,, X2 is 
exchanged with the pair x3, Xs, and x; and X¢, x; and Xs, etc., are also exchanged. This gives a con- 
tribution where 


exp[ —i(m-x,+n- x2) }} (exp [ —8Hi2 ]—exp[ — }) 
2 


and 4r(r—1) is the number of different ways of choosing two out of r pairs. The prime on the sum- 
mation of (22) indicates that k+1=m-+n; this restriction is valid because of the fact that an operator 
like exp—[8H2] does not affect the sum of the momenta of the two particles. Next there will be 
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a term in which » sets of two pairs will be exchanged with each other after the manner of y;: 
r![2"n!(r—2n)!]- (23) 


There will also be more highly exchanged terms in which several pairs will be mixed up among 
each other; these will have additional restrictions on them arising from the conservation of the 
total momentum of each pair of interacting particles. 

The relative orders of magnitude of the various terms is readily determined. We define a quantity 
x= Nw/Q; it is then a necessary condition on the validity of the use of the second virial coefficient 
that x1, although of course Nx>1.° We then have in order of magnitude: yo~ Nx, y1~ Nx", the 
latter being due to the restriction on the summation of (22). Similarly, all the more highly mixed 
terms will involve higher powers of x, although not of N, since the momentum conservation imposes 
more incisive restrictions on the summations in these terms. We need thus consider explicitly only 
the contributions from yo and y:; if the y: part can be neglected, this is a fortiori true for the higher 
terms. Summing (23) over m from 0 to $7, and substituting into (18) and (11), we obtain on reversing 
the order of summation: 


n=0 r=2n 
Since yi<vyo, we can neglect it in the physically interesting quantity: 


InW= In{1—exp[ — (a+ Bex) ]} +3 v.viG,(k, 1). (25) 


Now the grand canonical ensemble of (11) and (12) is a valid mathematical artifice only when 
the contribution to W from terms for which N’ is much larger than N is negligible. This is true of 
(25), since the main contribution to (24) comes from terms in the summand for which 27 ~ y)>~ Nx 
<«N. For these terms the omitted factor II; s;! is of no importance. But even for terms for which 
r~3N, the most this factor can contribute to the left side of (24) (this occurs when s~N) 
is ~(N!/(3.N) ~ (Px/e)**. This is small compared to as long as x1 
and # is not large compared to one. Our omission of the II; s;! factor is thus justified. 

A quantity exactly analogous to (25) can be obtained for a Fermi-Dirac gas of identical particles 
possessing spin 3h: 

InW’=2 In{1+exp[ —(a+ Be.) ]} +3 2 (k, 1; 11, (26) 


luyue 
Here 
(27) 


and 


F\(k, 1; 41, =exp[B(ex +e.) JQ-? f exp[ —i(k-x;+1-x2) ] 
—u*(2)u2*(1) expl —i(K-x2+1-x;) J} (expl — —exp[ 
X u;(1)u2(2) Jdx:dxo, (28) 


where the u’s denote normalized spin functions. The relations (13) to (17) remain valid. 

It should be noted that the formal simplicity of (25) and (26) derives essentially from the as- 
sumption that the particles are free except for their mutual interaction. Without this assumption 
there would be no momentum conservation and no justification for discarding the higher exchange 
terms. The correctness of (20) can now be established by putting (19) and (25) into (13). We then 
see that the left side of (20) differs from N only by terms of order Nix, which is small compared 
to N. Thus to the accuracy of this calculation, a is the same as for an ideal gas having the same 
temperature and density. 
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These results are readily extended to mixtures of different kinds of particles, or to particles of 
higher spin. Since such formulas would only be of interest for special problems, we do not give 


them here. 


3. EVALUATION OF G; 


The method outlined in this section for the explicit calculation of G; can easily be extended to 
F, and to the corresponding integrals for mixtures. It consists in expanding the plane waves that 
appear in the integrand in terms of the eigenfunctions of the two-particle Hamiltonian. Once these 
eigenfunctions are known, the exponential operator becomes quite tractable and G, can be reduced 
to a set of quadratures. Since with one exception the eigenfunctions can only be obtained by an 
analysis into spherical pat.ial waves, we discuss only this case. The exceptional case, in which the 


potential is of the Coulomb form, requires a separate discussion in any event. 
The integral G,(s, t) can be more simply expressed in terms of the sums and differences of the 
coordinates and momenta: x= }(x:+x:2), r=(x:—x:2), K=s+t, k=}(s—t). Performing the integral 


over x, we obtain for k~0 


Gi= (exp[2Bex ]/2) f (exp[ —7k- r]+exp[ik- — V)] —exp[8h?4/m)]] exp[ik- r jdr. (29) 


Substituting the expansion® 
exp[ikr cosé]= >>; (2/+1)i'j(kr) Pi(cosé) 
into (29), we obtain: 
Gi= (2 (2/+-1)? f ji(kr) Pi(cosé) [exp[6(h?A/m — V)]—exp[sh?A/m) 
Xjilkr)Pi(cosé)dr. (30) 
Suppose now that we know the regular eigenfunctions of 
1d 


r? dr dr h? r? 


+4) =0 


that have the asymptotic form for large r 
sin(kr+ 6,")/r. (31) 


We assume that there are no bound states for this potential; however, the theory can readily be 
extended to the case in which they exist. Then the u;'(r) form a complete orthogonal set of functions 


normalized so that 
f ux! *(r) ux = 6(k —k’), 
0 


and we can expand 
0 


where the integral is interpreted as the principal value, since ax,’ has a simple pole at k=k’. This 
is seen most easily by noting that 


f *(r) = RAR’, (33) 
0 


The “spherical Bessel function” j:(z) = (/2z)*J144(2). 
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where 
0 


is always finite for physically interesting potential functions. 

The quantity 8;' can be determined by fitting the asymptotic form of j;(kr), using (31) and (32); 
however, we now show that it does not enter into the final result. Since the integrand of (30) vanishes 
outside w we can extend the r-integration over all space. We substitute (32) for the second j,(kr) 
in the first exponential of (30), and for the first 7,(kr) in the second exponential of (30). The indicated 
operations in the integrand are then readily performed, and it is seen that the 6,' terms exactly 
cancel. The result is: 


Gx(s, t) = (Bx (2141) f | 


(exp — [Bh*k’?/m ]—exp —[Bh?k?/m k=}(s—t)¥0, (34) 


where again the integrand has a simple pole at k=’, and the principal value of the integral is 
implied. For k=0, the only nonvanishing contribution comes from /=0; it is easily seen to be 


Gi(s, s) = (4r/0) f | | (exp — [Bh?k’?/m (35) 
0 


Since u,/°(r) is proportional to k’ for small k’, this integral converges at the lower limit. 

For applications of this theory to situations where none of the states contains a large number of 
particles (no »,.>1), we can substitute (34) into (25) and not be concerned by the factor 2 difference 
between (35) and the limiting form of (34) as st. For this only makes an error of order }°; v.2/N’, 
which is negligible. This factor 2, which comes from the failure to exchange particles that are in 
the same state, only makes an appreciable difference when very many particles are in a single state. 
It is of course important to know that (35) is finite. 

It is worth noting that while (25) and (30) cannot be expressed in terms of the phase shifts 6,' 
alone, the limit of (25) for low densities (in which case (19) becomes the Maxwell-Boltzmann dis- 
tribution law) can be so expressed. The result is precisely that already obtained by Uhlenbeck 
and Beth.' In this limit, one regards the statistics as operative only between the pairs of colliding 
particles. 


4. APPLICATIONS 


It is clear from the foregoing analysis that there are three principal limitations on the applicability 
of this theory. (1) The wave functions for the elementary collision process must be known, preferably 
in analytic form; this restricts the type of two-particle potential that can be employed. More im- 
portant, it precludes the possibility of extending the formalism to the third virial coefficient, since 
accurate three-particle wave functions would be virtually impossible to obtain. (2) The interaction 
potential between a pair of particles must fall off sufficiently rapidly with distance so that the 
various integrals converge. This is not true classically for the Coulomb field, and it seems likely 
that an unscreened Coulomb field cannot be treated by this method; this would of course require 
a special investigation. (3) The particles must be free except for their mutual interaction. This 
restriction would be of little importance except in the case of the electrons in a metal, where one 
would have to use plane waves rather than Bloch-type wave functions. 

In spite of these limitations, the correction to the ideal gas models for such condensed systems 
as liquid helium and heavy atomic nuclei might be expected to yield interesting results, and the 
application to liquid helium is given in the following paper.* 
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On the Phase Transition in Liquid Helium 


L. I. Scuirr 
Randal Morgan Laboratory of Physics, University of Pennsylvania, Philadelphia, Pennsylvania 
(Received January 27, 1941) 


The method developed in the preceding paper for treating non-ideal quantum gases is 
applied to a gas of helium atoms having the density and temperature of liquid helium. The 
object is to investigate the influence of binary collisions on London's ideal Einstein-Bose gas 
model of liquid helium. It is shown that the condensation temperature and the continuity of 
the specific heat at this temperature are not affected. These results are independent of the 
force law assumed to hold between a pair of atoms. Explicit calculations for a rigid sphere 
interaction of small radius show that the specific heat is lowered near the condensation temper- 
ature. It seems possible that the phase transition observed in liquid helium is characteristic 
of a condensed system, and cannot be obtained from a gas approximation method such as 


that employed here. 


1. INTRODUCTION 


ONDON has suggested that the condensation phenomenon of the ideal Einstein-Bose gas, dis- 
covered by Einstein, is connected with the “‘A-point”’ phase transition observed in liquid helium.' 
This is supported by the remarkable agreement between the temperature of condensation, 3.14°K, 
and the temperature of the phase transition, 2.19°K. On the other hand, the specific heat of the 
Einstein-Bose gas exhibits a discontinuity of slope, but not of magnitude, at the condensation tem- 
perature, while the specific heat of liquid helium probably possesses a discontinuity in magnitude at 
the transition temperature.? In an effort to bring the theoretical and experimental specific heat 
curves into agreement, London postulated that the interatomic forces can be approximately taken 
into account by altering the distribution of the lowest energy levels of the (independent) atoms.’ 
By empirically fitting three parameters and adding a Debye term, he was able to reproduce the 
experimental specific heat curve quite well. This procedure appears to be open to two basic objec- 
tions. (1) The lowest energy states of a particle moving in a periodic field have a density proportional 
to the square root of the energy, and not to a higher power, as assumed by London. (2) The condensa- 
tion temperature is no longer regarded as a quantity to be derived from the theory, but merely as an 
adjustable parameter. It would thus seem that the greater part of the physical content of the earlier 
theory has been discarded. 

A method for considering interatomic forces consistently to an accuracy corresponding to the 
classical second virial coefficient while treating the degeneracy exactly is provided by the preceding 
paper.‘ The latter point is essential, since the appearance of the condensation phenomenon depends 
on the exact treatment of the degeneracy (use of completely symmetrized wave functions in calcu- 
lating the partition function). While it is not to be hoped that this inclusion of binary collisions will 
account for the properties of a condensed phase, one might expect that such a treatment will at least 
indicate the nature of the deviations from the ideal gas model caused by interatomic forces. The 
method was fully developed in the preceding paper only for the cases in which no states of the corre- 
sponding ideal gas contain a number of particles large compared to one. However, below the con- 
densation temperature of the ideal gas, the lowest energy state contains on the average a finite 
fraction of all the particles. This circumstance not only necessitates an extension of the method of 
the preceding paper, but also requires that the use of the grand canonical ensemble be abandoned. 


1F, London, Phys. Rev. 54, 947 (1938). 

2? W. H. Keesom and A. P. Keesom, Physica 2, 557 (1935). 

*F. London, J. Phys. Chem. 43, 49 (1939). , 

*L. I. Schiff, Phys. Rev. 59, 751 (1941); preceding paper. Primed equation numbers refer to this; the same notation is 
used in both papers. 
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The reason for this is that the grand canonical ensemble predicts fractional density fluctuations of 
the order of one in this case, and hence does not provide a physically satisfactory description of the 
system.’ The saddle-point approximation in the method of Darwin and Fowler also breaks down 
here.° We therefore adopt a procedure which is closely related to that used in finding the most 
probable state of a microcanonical ensemble.’ 


2. THe Speciric HEAT 


We shall deal with the ordinary partition function 


iN 


r=0 


where Z, is given by Eqs. (8’) and (9’) ;4 this may also be written 
Z(nx). (1) 


Now In Z(n,) can be evaluated to the accuracy of terms of order Nx by the method used in Section 2 
of the preceding paper. If mo at most is of order N, all other n, being small compared to N although 
not necessarily of order one, it can be shown that 


In Z(nx) = —B +3 n,niG,(k, 1) +7 k) + 3,°G,(0, 0). (2) 


k, k#0 


We now divide the independent particle states into groups g,>1 having a mean propagation vector 
k, for each of which 
= > 


king 


To the accuracy employed, In Z(n,)=In Z(n,’), and is given by (2) with n,. replaced by n,’. Instead 
of (1), we now have 


n, ’ 7 
« m,'\(g.—1)! 
where the index y’ represents a set of integers n,’ that are subject to the restriction 
>. =N. (4) 


The summand of (3) has a sharp maximum for a set of values 7,’ for which its variation and the 
variation of (4) both vanish. Combining these requirements by Lagrange’s method of undetermined 
multipliers and using Stirling’s formula on the factorials, we obtain at the maximum 


i,’ =g, expla’ +Be,— 1) (5) 


The undetermined multiplier a’ must be chosen so that the relation obtained by substituting (5) 
into (4) is valid. The physical significance of (5) is that any quantity that can be calculated from Z 
can be estimated by considering only the contribution to Z from the neighborhood of the maximum 
of its summand. 

Equations (5) and (19’) agree in form except for the correction term in the exponent of the de- 
nominator of (5), which is of order x1. Thus we have fi,’ = g,», above the condensation temperature, 
which is not changed to this accuracy. Below the condensation temperature, we must have fio’ of the 


5 The grand canonical ensemble gives ((N*)ay—(N)a?)/( Na? = 1/ (nya? /(N)x?, which is of order one if any of 
the are of order (N)yy. 
_ In the contour integral method of Darwin and Fowler it is no longer convenient to expand the logarithm of the 
integrand about the saddle-point, since the resulting series converges slowly, if at all. 

7Cf. R. C. Tolman, The Principles of Statistical Mechanics (Oxford, 1938), Chapter 10. 
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order of magnitude of N in order to satisfy (4). This means that 
a’ k), 


instead of a~0, as it is for the ideal gas. Then 7,’ =g,[Be:—7i0'G,(0, 0) |" (which is positive: see 
Eq. (35’)); this is somewhat smaller than the ideal gas value giv1=gi(8e:)~!, but this makes little 
difference since in any case it is very small compared to N. It can now be shown that it is actually 
vio that is of order N, and not 7, etc. This is done by treating the lowest group go separately and 
without using Stirling’s formula on it; since vio’ turns out to be independent of go for go=1, it is clear 


that only fio is very large. 
Thus in general, most of the Z in (1) comes from nm; ~v,;, where »; is given by (19’) with 


a=a'—>, 7,'G,(0, k). 
The average energy of the system may then be written 
f= E,Z(nx) dX. 
Z(nx) Lu Z(nx) 
From (6) and (2) we obtain 
E=> 1)/0B—v0 v.0G,(0, —302dG,(0, 0) /a8. (7) 
k k, 10 k#0 


= In vn. (6) 


The specific heat at constant volume is . 
C,=dE/dT. (8) 


3. NUMERICAL CALCULATIONS 


In this section we calculate E explicitly, making the not uncommon assumption that the inter- 
action between a pair of helium atoms at a distance r can be represented by the potential 


r<a; V(r)=0, r>a; (9) 


a is supposed to be small compared to the mean wave-length h/(mkT)} of the temperature motion of 
the atoms. This is no less correct for liquid helium than the assumption of binary collisions that 
underlies the use of the second virial coefficient. Moreover, it enables us to simplify the calculation 
by considering only S collisions. We show later that the main results of this calculation carry over to 


a more exact treatment. 
Referring to Section 3 of the preceding paper,‘ we have from (31’), taking /=0: 


sink(r—a)/r. (10) 
From (33’), with the limit V-+ « for r<a, when a is assumed to be small, we have 
= (2/2) sinka/k (2/1) 
With the substitutions y? = Bh?k"?/m, b=Bh*?k?/m>0, we obtain from (34’): 


G,(s, t) = (11) 
where 


400) = My. (12) 


The second term in the integrand is added for convenience, since it removes the singularity at y?=) 
and its principal value vanishes. The integral in (12) can then be expanded about 6=0 by Maclaurin’s 
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theorem, vielding a series that converges quite rapidly: 


f(b) =¥ (13) 
We also have 
Ga(s, 8) = — (14) 


We can replace }>; in (7) by fQdk/(27)* and take the lower limit of the integral to be zero, since 
the contribution from this neighborhood is very small. We can also take »=0 for T>T7 ) 
=(N/2.612V)8h?/2xmk,! and and given by (19’) with a=0 for In 
calculating C, from (8), we must remember that for 7>T>, d/dT =0/0T+(da/dT)d/da, where 


da/dT = —(@N/dT)/(0N/da), (15) 
and 


N=E n= f Odk (16) 


It is readily seen that 0N/da@ becomes infinite as a—0 (T—T> from above). This means that in calcu- 
lating C, for T just above Ty we need retain only those terms in 8E/da which diverge in the limit 
a—0, and cancel the divergent terms before taking the limit. The rest of the calculation is straight- 
forward, although tedious, and gives: 


C,=R[1.932!+ ]; (17) 
(18) 
(dC,+/dT)r=ro= — (R/T») [0.77 +20.6(Na*/Q)#]; (19) 


here, z= 7/7 )=1. The superscripts + and — denote values just above and just below 7», respec- 
tively ; the first terms in (17) and (19) are the ideal gas values of London.! 

The appearance of the term z~! in (17), and the validity of (18), are the two striking features of 
this rigid sphere calculation that we should like to consider more generally. The z~} term comes from 
the last factor in (7), which contributes (aside from another part) an amount (v¢?/2k7T*)(8°G,(0, 0) /d8*). 
From (35’), we see that in general this is 


f | Vox®|? exp—[Bh2k’2/m (20) 
0 


For small T (large 8), practically all of the integral in (20) comes from very small k’, for which V,-° 
is proportional to k’; the integral is then proportional to 8-!. Thus for small 7, the contribution (20) 
to C, is positive and varies as 7~*. A similar argument shows that no other term in C, becomes 
infinite as 7-0. This infinity in C, (and hence in the entropy) at T=0 must be regarded as an 
indication that the second virial approximation breaks down here. Indeed, from (2) and (35’) one 
sees that for small T the ideal gas part and the correction term in In Z are of orders Nz! and Nxz-}, 
respectively. Since these represent the first two terms of a series in virial coefficients, one would 
expect the expansion to be useful only when 2? >x. The calculation is thus valid only for T>7T.~x'T9; 
since we always assume that xX1, we have that 7.&T». 

To examine the validity of (18) in general, we need consider only those parts of (7) whose contribu- 
tions to (8) are different just above and below 7». The first term of (7) is the ideal gas energy, and its 
specific heat has been shown by London! to be continuous at 7». The last term of (7) gives no contri- 
bution to C, on either side of To, while the third term gives a contribution only below To, of magnitude 


f ]—1), (21) 
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where ¢(s—t) = —0G;,(s, t) 08. The 0/08 of the second term of (7) contributes the same amount to 
C, on both sides of 7», while the 0/da of this term gives a contribution only above 7». We thus wish 


to compare (21) and 


lim f f (22) 


in the general case. With the help of (15), we can rewrite (22): 


im fa expLat+ fe, ]—1)? ] 


Since both the q and 1 integrals in (23) diverge for small values of their arguments in the limit a0, 
the only nonvanishing contribution to (23) will come from very small 1. Expanding ¢(k—1) about 
1=0 and canceling the divergent integrals, we obtain in the limit: 


f [Boex ]—1). (24) 


For T=T> we have that N+ »o=constant (the true total number of particles), where N is given by 
(16) with a=0. Differentiating this with respect to T, we see that 0N/dT+év)/dT=0. Thus (24), 
and hence (22), is equal to (21), and the specific heat is continuous at 7». 


4. Discussion 


The information contained in Eqs. (17), (18), and (19) is plotted in Fig. 1 for three choices of the 
rigid sphere radius a. The largest value of a used (1.3A) is still considerably smaller than the best 
known radius of the repulsive interaction.* The 
more nearly correct value of a@ would make C, 
gy negative near 7; this is not surprising in view of 
n . a0 the fact that we assumed a to be small in calcu- 
lating (11) and (14). It would be perfectly feasible 
to calculate C, without this assumption (since the 
wave functions for a rigid sphere potential are 
elementary functions for all partial waves) ; how- 
ever, this does not appear to be worth the effort in 
the light of the general results proved in the last 
section. These may be summarized: 

A calculation made to an accuracy correspond- 
ing to the classical second virial coefficient, in which 
quantum degeneracy and diffraction effects are 
considered exactly (use of completely symmetrized 
wave functions), shows that the condensation tem- 
=" perature is unchanged and that the specific heat at 

” constant volume iscontinuousat this temperature, 

Fic. 1. The specific heat at constant volume, regardless of the form of the (spherically symmet- 

considering only binary collisions. Although 
ric) interaction between a pair of helium atoms. 


these curves eet a rigid sphere interaction 
of radius a in A, the continuity of C, at T=7> -- 
8 H. Margenau, Phys. Rev. 56, 1000 (1939). 


remains for any interaction. 
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The marked disagreement between these results and the observed specific heat of liquid helium* 
suggests the possibility that the experiments cannot be understood in terms of a finite number of 
virial coefficients, that is, in terms of a gas model in which the helium atoms are supposed to be 
free most of the time, with corrections being made for binary, ternary, etc., collisions. It would 
seem rather that the condensed character of the system must be taken into account from the begin- 
ning, by using an approximation based on a liquid model. In any event, this would appear to be the 
more constructive approach, since the greatest difficulties would be encountered in extending the 
gas approximation to virial coefficients higher than the second. It seems likely that the decisive role 
played by the Einstein-Bose statistics, and the remarkable agreement between the zero-point energy 
of the atoms and the thermal energy at the transition temperature, will also appear in a successful 


liquid model. 


MAY 1, 1941 


PHYSICAL REVIEW 


VOLUME 59 


The Energy Distribution and Composition of the Primary Cosmic-Ray Particles 


NORMAN HILBERRY 
University of Chicago, Chicago, Illinois 
(Received April 3, 1941) 


HE cascade computations for the experi- 

ment on the variation of the number of 
extensive cosmic-ray showers with altitude as 
previously reported' have been repeated for a 
wider range of variables and of the exponent in 
the expression for the energy distribution of the 
primary particles. The best value for this expo- 
nent as indicated by these calculations is 2.70 
or slightly larger, although values from 2.60 to 
2.80 cannot be excluded. If the value 2.75 is 
chosen, the energy distribution is given as 


particles /(cm* Xsec.). 


With the above choice for the energy distribution 
and with the assumption that the particles are 
isotropically distributed at the top of the atmos- 
phere, integration gives 0.09 particle per cm* 
per sec. from all directions on top of the atmos- 
phere as the number with energies in the range 
0.67 X10" ev to 1.7 X10" ev. Integration for the 
energy carried by the particles described by this 
expression shows that for those with energies 
above 1.70X10'° ev the energy carried into the 
atmosphere per cm? per sec. from all directions 
is 8.9X10° ev. Bowen, Millikan, and Neher 
give** for these quantities the experimental 
1'N. Hilberry, Phys. Rev. 59, 112A (1941). 


? Bowen, Millikan and Neher, Phys. Rev. 52, 855 (1938). 
*R. A. Millikan, Cosmic Rays (Macmillan, 1939). 


values of 0.09 particle /(cm*Xsec.) and 9X 10° 
ev /(cm* Xsec.), respectively. Of this experimen- 
tally observed energy they tentatively ascribe 
one-half to a non-ionizing primary radiation. 
The present calculation suggests that this as- 
sumption may not be necessary. 

No claim is made for the precise validity of the 
expression given here. However, when our new 
data are used, the determination of the coefficient 
as a function of the exponent is considerably 
more accurate than the choice of a particular 
value of the exponent. The value 2.75 was chosen 
since the expression then not only describes 
satisfactorily the high energy shower phenomena, 
but also yields a quantitative agreement with 
observations at lower energies. The extensive 
showers observable with the present apparatus 
arise almost entirely from primaries with energies 
above 5X10" ev. The fact that an energy dis- 
tribution determined for primary particles in 
this energy range is capable of giving quantita- 
tive agreement with observations depending 
upon primaries with energies of 10'° ev indicates 
that there is but one kind of primary cosmic-ray 
particle in any significant number. The experi- 
ments of Schein, Jesse and Wollan just reported‘ 
agree with this conclusion and show that this 
particle is most probably the proton. 


‘Schein, Jesse and Wollan, Phys. Rev. 59, 615 (1941). 
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The choice of the proton as the primary would 
appear to account for the observed number of 
extensive showers at low altitudes, which is 
definitely too high for cascades caused by pri- 
mary electrons. At the large energies required for 
the observed showers (>5 X10" ev) it is not un- 
reasonable to suppose that the radiation cross 
section of the proton would be large. Under some 
such assumption, it may be supposed that most 
of the energy of the proton is spent in exciting 
mesotrons and photons high in the atmosphere. 
The energy which the proton transmits into 
a typically developing cascade shower should 
then be nearly the same as that transmitted 
by a primary electron of the same initial energy. 
This equivalence should be even more pro- 
nounced since any mesotrons produced by the 
primary proton in this energy range should also 
radiate until their energies fall below some 
critical value.5 Corresponding to this approxi- 
mate equivalence in the energies transmitted into 
the cascade, is an equivalence in the resulting 
spatial distributions of particle density. In the 
regions in which the showers are fully developed, 
the differences between the two distributions 
would probably be less than the present experi- 
mental inaccuracies. Now, however, due to the 
production of mesotrons by the proton, there 
should be associated with the shower a penetrat- 
ing component consisting largely of mesotrons 
with energies below their critical energy. This 
would account for the penetrating component in 
a previously reported experiment on extensive 
showers.* These mesotrons, together with any 
which may be produced by photons,’ furnish a 
mesotron distribution in the shower which will 
persist to much greater depths in the atmosphere 
than the cascade electron distribution. This 
would account for the increase in the experi- 


Christie and S. Kusaka, Phys. Rev. 59, 414 
(1941). 
, ®N. Hilberry and V. H. Regener, Phys. Rev. 59, 471A 
1941). 

7 Schein, Jesse and Wollan, Phys. Rev. 57, 847 (1940). 
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mental counting rate for extensive showers ob- 
served at low altitudes over that predicted by 
the electron cascade calculation. 

The primary energy distribution, together 
with the modifications produced by proton and 
mesotron radiation processes, would appear 
capable of accounting qualitatively for the ob- 
served mesotron energy distribution at high 
energies. Analogy with photon production by 
electrons, and with pair production by photons 
suggests that the mesotron energy distribution 
should have the same power law exponent as the 
primary protons. The measurements of V. C. 
Wilson® represent an integrated mesotron energy 
distribution with an exponent of 1.77 for the 
energy range corresponding to depths down to 
250 m water equivalent compared to the value 
of 1.75 for the present case. For greater depths, 
Wilson gives the exponent 2.52. Qualitatively, as 
Wilson points out to me, this change in the 
exponent is what would be expected if the 
mesotrons lose energy by radiation when pro- 
duced with energies above their critical value. 
That is, the number of mesotrons in any energy 
range above this value would be reduced, the 
reduction depending upon the magnitude of 
their radiation cross section at that energy. As 
the critical energy is approached closely enough 
for the contribution to the total number of 
particles from the energy range above this 
value to become important, the exponent of the 
integrated mesotron energy distribution must 
increase. Actually, the change of exponent in 
Wilson’s observations is not abrupt but increases 
rapidly from 1.77 to 2.52 in the neighborhood 
of 250 m water equivalent. 

These results will be published in detail in the 
near future. 

In conclusion, the author wishes to express his 
indebtedness to Professor A. H. Compton for his 
interest and advice. 


8 V. C. Wilson, Phys. Rev. 53, 337 (1938). 
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Letters to the Editor 


ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. Closing dates for this department are, for 
the first issue of the month, the eighteenth of the preceding 
month, for the second issue, the third of the month. Because 
of the late closing dates for the section no proof can be shown 
to authors. The Board of Editors does not hold itself responsible 
for the opinions expressed by the correspondents. 
Communications should not in general exceed 600 words 
in length. 


A Directional Geiger Counter 


H. G. STEVER 
California Institute of Technology, Pasadena, California 
April 8, 1941 


URING the last half-year in a study of the discharge 
mechanism of fast or self-quenching Geiger counters 
a new directional Geiger counter was discovered. Out of 
this study, among other things, come two facts which 
explain the operation of the directional Geiger counter. 
First, the discharge of a counter spreads along the unob- 
structed wire such that the charge collected on the wire 
in a given pulse is proportional to the length of the counter. 
Secondly, a small glass bead on the wire of diameter 
(approximately) a few times the diameter of the wire is 
sufficient to stop the spread of this discharge along the wire. 
The directional Geiger counter consists of an ordinary 
Geiger counter with beads of glass on the wire dividing the 
counter into sections, preferably equal. It is to be under- 
stood that the beads do not insulate the sections of wire 
but are merely put around the single central wire. Also 
it is to be noted that this is an organic-vapor-filled counter 
with treatment of the copper wall of the cylinder to give 
very low photoelectric emission efficiency and high work 
function. Oxidation by heating the counter with NO» gas 
will give such a surface. 


TABLE I. Triple counts per five minutes. 


WITHOUT 7-RAY Wi1TH 20,000 

BACKGROUND 7-RAY SINGLES 
1140 1145 
1131 1144 
1139 1117 
1145 1150 
1143 1202 

Av 1140415 1152315 


If an ionizing particle passes through one section of the 
counter this section by normal counter action furnishes a 
charge g to the wire, this charge depending on the length 
of the section. The other sections do not enter into the 
discharge. If an ionizing particle passes through two 
sections, then the charge furnished to the wire is 29 (for 
equal length sections) and for three sections 3g and so on, 
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Hence for a given wire capacity, C, the voltage pulse is 
nq/c where n is the number of sections through which the 
particle passed. From this the directional effect is im- 
mediately obvious, for if, in the recording circuit, the first 
amplifier tube is biased so that it passes on only pulses of 
certain minimum voltage, mog/c, where mo is the total 
number of sections, then only ionizing particles which 
have passed along the tube through all sections will be 
recorded, This then is a true Geiger counter telescope. 

To show this directional effect an electron source was 
placed at the end of a two-beaded counter, so that singles, 
doubles and triples were expected. The voltage pulse from 
the wire was put on the vertical plates of an oscilloscope 
and a linear sweep on the horizontal. Figure 1 is a one- 
minute exposure of the oscilloscope screen showing the 
three heights of pulses. The individual pulses are not 
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Fic. 1. Time exposure of oscilloscope screen showing singles, 
doubles and triples. 


distinguishable since it is a time exposure. It is important 
to note the sharp dividing line between single, double and 
triple counts. 

Table I gives results of a simple experiment to show 
that this counter may be used to count a directional 
ionizing ray where there is a large general background 
radiation with no hindrance from the background. With 
the electron source at the end of the tube the number of 
triples was recorded for a five-minute interval. Then on 
the perpendicular bisector of the axis of the tube a y-ray 
source was placed such that the increase in single counts 
was 20,000 counts per five-minute interval and again the 
triples were recorded. There were 22+4 triples from the 
y-ray source with electron source removed. 
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On the Theory of the Temperature 
Diffuse Scattering 


W. Hl. ZACHARIASEN 
Ryerson Physical Laboratory, University oj Chicago, Chicago, Illinois 
March 20, 1941 


ig a recent article! G. E. M. Jauncey has followed up 
in greater detail earlier suggestions by G. D. Preston 
and by W. H. Bragg that the maxima of the temperature 
diffuse scattering can be “explained” on the assumption 
that the thermal agitation breaks the crystal lattice up 
into small groups of atoms. The present writer finds it 
difficult to understand how it is justifiable to assume the 
existence of such independently scattering small groups of 
atoms since this assumption is in direct contradiction to 
the observed sharpness of the Laue-Bragg reflections. 
Clearly one cannot use one crystal model for the sharp 
maxima and another for the diffuse maxima. One is led 
to believe, therefore, that both crystal models coexist, 
i.ec., that the crystal is composed of two types of con- 
stituent blocks, one type of blocks containing millions of 
atoms, the other type of blocks containing only one 
hundred atoms or less. This conclusion is improbable to 
say the least and cannot merit serious consideration until 
it is supported by convincing evidence. 

Although one cannot accept the physical picture upon 
which the Preston-Bragg considerations are based it is 
nevertheless true—as Jauncey demonstrates—that their 
assumption gives the approximate positions of the diffuse 
maxima. The reason for this apparent success is readily 
understood. The explanation is indeed obvious when we 
compare the equations used by Jauncey with the ones 
deduced by the writer.? Equation (8) of Jauncey’s article 
can be written in the following form 

T=const X exp(is-r) |, [J8] 
where s = 27(k—Ko), r= D/;a; and where (actually /3) 
assumes values 0, 1, --» N—1. On the other hand Eq. (8) 
of reference 2 may be written in a similar manner as 


J»=const X | Lic: exp(is-r)) |’, [Z8] 


where the coefficient c; is so defined that cicv = Pir (com- 
pare reference 2). As Py rapidly decreases as /—l’ in- 
creases (see Eq. (11) of reference 2) it is immediately 
apparent that Eq. (J8) is a somewhat crude approximation 
to Eq. (Z8). 

Jauncey in his article uses the term modified reflection 
of x-rays by crystals. This writer sees no apparent reason 
for rejecting the well-established term diffuse scattering. 
Diffuse scattering as used in the past by Jauncey and 
others is defined as the difference between the total 
scattering and the Laue-Bragg scattering, i.e., as the 
difference between the mean square and the squared mean 
amplitude. In normal crystals where the only disorders 
correspond to the electronic motions within the atoms and 
to the heat motion of the atoms the diffuse scattering 
accordingly consists of two terms, namely the Compton 
diffuse scattering and the temperature diffuse scattering. 
This definition leads directly to the expression (Z8) for 
the temperature diffuse scattering. As shown by this 
formula the intensity is not zero for any scattering di- 
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rection (except, of course, in the exact forward direction) 
and the adjective diffuse is thus truly descriptive of the 
phenomenon. 

The formula for the shift of the diffuse maxima from 
the Bragg positions given as Eq. (23) of reference 2 should 
not be considered exact since it was deduced from the 
approximate intensity formula (Eq. (18) of reference 2) 
rather than from the accurate formula (Eq. (17) of 
reference 2). Indeed, more rigorous treatments (which 
will soon be reported) show that for NaCl the shift is 
about 40 percent smaller than demanded by Eg. (23) 
and also that the half-width of the maximum is about 
30 percent smaller than indicated by Eq. (25). 


!G. E. M. Jauncey, Phys. Rev. 59, 456 (1941). 
2W. H. Zachariasen, Phys. Rev. 57, 597 (1940). 


On the Absence of the M8 Satellite 
Intensity Anomaly 


FREDERICK R. Hirsu, Jr. 
Pasadena, California 
April 10, 1941 


ARLY in 1936 the writer predicted intensity anoma- 
lies, due to the Auger effect, for the satellites of the 
x-ray diagram lines Ma; and MB.' Very recently, the 
present author was able to demonstrate by some M series 
plates? that there actually was an intensity maximum, 
due to the Auger effect, for the Ma satellites in the neigh- 
borhood of Z =82.' More recently, the writer again sought 
out these same plates? to investigate the intensity of the 
MB satellites in the region of Z=84, where the intensity 
anomaly should exist. (The Coster-Kronig diagram, Fig. 1, 
shows that the ionization energy curve for the Nyy, shell 
is intersected by the curve for the energy released by the 
radiationless transition M);;—~Myy at about Z=84.) Thus 
the M8 satellites originate from the single electron transi- 
tion between doubly ionized states: Mjy Niy,v>Nvi 
Moreover, there is every reason to believe that this is the 
correct transition, since in the case of the Ma satellites 
the correct transition is known to be My Niy,v~>NvyiNiv.v, 
and the Semi-Moseley graphs for these two groups of 
satellites have been shown by the writer to be nearly 
identical.? In addition, the Ma; and Mg line profiles have 
been shown by the writer to be very similar.' 


Atomic Number 


Fic. 1. The energy released by the radiationless transitions 
and and the ionization energies for the 


N shells (for Z +1) plotted against atomic number. 
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It was, therefore, an occasion of great surprise when 
these plates* definitely showed that the predicted anomaly was 
not present. Inspection of Fig. 1 shows the presence of 
another radiationless transition energy curve (My,;~Mry) 
intersecting the ionization energy for the N, shell at the 
same atomic number where the energy curve for the 
radiationless transition M,;;—>Myy intersects the ioniza- 
tion-energy curve for the Nyy,v shell. Inspection of line 
widths, as recently measured by Dr. J. N. Cooper‘ reveals 
the fact that the radiationless transition M;;~Myy is 
clearly more probable than the radiationless transition 
Mur>Mvw, and in robbing the less probable transition 
(Mur>Mr), effectively prevents the occurrence of the 
MB satellite intensity maximum below Z=84. The facts 
will be presented more adequately in a forthcoming paper 
in this journal. 


1F. R. Hirsh, Jr., Phys. Rev. 50, 191 (1936). 
2? F. R. Hirsh, Jr., Phys. Rev. 38, 914 (1931). 
3 F. R. Hirsh, Jr., Phys. Rev. 57, 662 (1940). 
4 J. N. Cooper, Thesis, Cornell University, 1940. 


The Influence of End Effect on the Theoretical 
Shape of Intensity Curves in Radio- 
frequency Spectra 
A. F. STEVENSON 
Department of Applied Mathematics, University of Toronto, 
Canada 
April 7, 1941 


I‘ a recent paper! on the theory of Rabi’s radiofre- 
quency spectrum method of determining magnetic 
moments, the “end effect’’ considered by Millman? was 
worked out quantitatively. It was pointed out, however, 
that the results required some modification on account 
of the inhomogeneity of velocities in the beam. The effect 
of this inhomogeneity, if a Maxwellian distribution is 
assumed, has since been worked out by Torrey,’ whose 
results enable the end effect analysis to be completed. 
Torrey has also shown how the calculations for the simple 
case J=0, J=} may be extended to the general case. 
Torrey’s results (with obvious extensions and with a 
change of notation) combined with those of reference 1 
lead to the following relation‘ for the ratio of beam inten- 
sity (J) to intensity in the absence of transitions (Jo), 
in the neighborhood of resonance for the transitions 
n<m: 
(1) 
where 


y= | |L/2hv, B=(|E,—En| —hw) / | Vam|. 


Here @ denotes the ratio of the number of molecules in 
states m, m to the total number of molecules in the beam 
(measured in terms of intensities), E,(H), En(H) the 
energies of states m, m in the uniform field of strength H, 
Vam Cos(wt+5) the matrix element of that part of the 
Hamiltonian due to the oscillating field of frequency 
w/2r, 1 the length of the oscillating field, v the average 
velocity’ in the beam, and y a small constant, determined 
by the magnitude (distance between wires)/length of 
oscillating circuit (y= 2Cb/I in the notation of reference 1). 
It is assumed that there are no other states in resonance 
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with m or m at the frequency w/2. The function f is 
defined by 


Vat, g(x) = cos(x/y)dy. 


Torrey discusses his results in terms of the function g, but 
it seems rather simpler on the whole to use the function /. 
A table of values of f(x) may easily be constructed with 
the aid of Torrey’s tables. 

The factor (1+7y6y) in (1) is due to the end effect. It 
is assumed that the oscillating circuit is of the type given 
in Millman's paper,’ and that the direction of the uniform 
field is that referred to by Millman as the “normal” one. 
The + or — sign in (1) is then to be taken according as 
E,,(H) is greater or less than E,,(H), where m is assumed 
to be the state which has the smaller component of angular 
momentum in the direction of the uniform field. In many 
cases this is equivalent to the statement that the + or — 
sign is to be taken according as the magnetic moment is 
positive or negative. We shall assume this to be so in the 
following, and that the direction of the uniform field is 
the normal one. If the uniform field is reversed, the signs 
are reversed. The result (1) does not hold for very small 
or very large values of y. 

Working to the first order in y, we find from (1) that 
resonance (minimum J) now occurs, instead of at 8=0, at 


B=Bres = (2) 


Since f’(¥) <0 for all y, the shift is towards the high field 
side if the magnetic moment is positive,® in agreement with 
experiment.2 For y=0.6r=1.89, which according to 
Torrey is the optimum value, (2) gives Bres = +0.94y. 

As regards the asymmetry of the curve defined by (1), 
the difference in the width, measured from resonance, on 
the high field and the low field side at half-minimum is, in 
terms of 8, 


A= /f') — (3) 


where @ is defined by f(¢)=43f(y). The expression in 
square brackets in (3) is found to be positive for all y, so 
that the curve is steeper on the low field than the high 
field side if the moment is positive, also in agreement with 
experiment.2 This clears up the discrepancy noted in 
reference 1. For y=1.89, we find A&+0.6y. This is close 
to the minimum value of A. 

Equations (2) and (3) enable the influence of end effect 
on the position of the resonance point and the asymmetry 
to be predicted for any particular case, since y will be at 
least approximately known and the constant y can be 
found once and for all for a given apparatus by comparison 
with experiment. 

1A. F. Stevenson, Phys. Rev. 58, 1061 (1940). 

2S. Millman, Phys. Rev. 55, 628 (1939). 
?I am much indebted to Dr. Torrey for kindly letting me see his 
paper in advance of publication. 

4In deriving (1) from Eq. (22) of reference 1, the factor (1+ 8y) 
has not been averaged over the velocity distribution. This is justifiable, 
since it is the oscillatory factor in (22) which is the important one to 
average, and in any case any error committed could be approximately 
compensated by multiplying y by a factor. 

5 More precisely, » =(2k7/m)% is the most probable velocity in the 
Maxwell distribution. 

6 Since, with our assumptions, @ is of the form K(H —Ho) in the 


neighborhood of resonance, where K is a positive constant, “High 
field"’ corresponds to ‘‘low frequency”’ if w is varied instead of H. 
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Young’s Modulus of Elasticity and Its Change 
with Magnetization in Iron-Cobalt Alloys 


Mikio YAMAMOTO 


Research Institute for Iron, Steel, and Other Metals, 
Téhoku Imperial University, Sendai, Japan 
March 24, 1941 


HE only available data on Young’s modulus of 

elasticity and its change with magnetization, i.e., 
the so-called AE effect, measured in the whole composition 
range of iron-cobalt alloys, are those obtained by Honda 
and Tanaka! fifteen years ago. The method of measurement 
adopted by them is the static method of bending. The 
author has investigated the same subject by the use of 
the dynamic method of magnetostrictive oscillation.? The 
specimens used were those on which Professor H. Masu- 
moto? previously measured the intensity of magnetization. 
They were prepared in the following way: weighed Armco 
iron and electrolytic cobalt were melted and cast, and the 
ingots were turned to cylindrical rods of diameter about 
0.4 cm and of length about 10 cm. These specimen rods 


had previously been annealed at 1000°C for one hour in a 


hydrogen atmosphere, and the author again heat treated 
them at the same temperature in vacuum before the 
measurements were carried out. The measurement of the 
AE effect was first conducted in the magnetic field up to 
900 oersteds, and then, by extrapolating the measured 
data to the null magnetic field, the absolute values of 
Young’s modulus of elasticity were determined. 

The relation between the measured maximum values of 
the change of Young’s modulus with magnetization, 
relative to that in the nonmagnetized state, (AE/E) max, 
and the composition (in weight percent) are shown in 
Fig. 1. (AE/E) max is +0.234 percent for Armco iron, and 


We 
3 

YAMAMOTO 
18 
por Fic. 1. Relation be- 
tween the maximum 
We values of the relative 

ulus of elasticity wit 

vo} magnetization and the 
| cobalt alloys. 

33 72 HONDA-TANAKA’ | 

10 20 30 @ 


increases more and more rapidly with cobalt content; 
finally it makes a sharp maximum near 50 percent cobalt, 
which amounts to over +22 percent.‘ Then the relative 
change of the modulus decreases quickly, but shows two 
more small maxima before reaching the value of +0.168 
percent for electrolytic cobalt. The measured data (for 
about 400 oersteds) of Honda and Tanaka, as shown by 
dotted lines in the same figure, agree roughly, though 
qualitatively, with our result, but there are two remarkable 
differences between them: (1) A high and sharp peak 
exists near 50 percent cobalt in our curve, while: this 


768 LETTERS TO THE EDITOR 


maximum appears as a rather flat and broad one in 
Honda and Tanaka’s measurement, and (2) Honda and 
Tanaka found the AE effect of pure cobalt to be negative 
in the whole of the magnetic field range up to 400 oersteds, 
but, contrary to this, our measurement indicates that it is 
always positive. These discrepancies are considered to be 
primarily due to the difference in the methods of measure- 
ment adopted. As is well known, in the static method for 
measuring the modulus of elasticity, as used by Honda 
and Tanaka, it is, in general, necessary to load beforehand 
a considerable weight as a dead weight on a specimen rod. 
Accordingly, what this method measures is the modulus of 
elasticity or its change in a considerably stressed state. 
Next, Young’s modulus of elasticity at nonmagnetized 
state is plotted against composition in Fig. 2. The curve 


Fic. 2. Relation between 
Young's modulus of elas- 
ticity and the composition 
in iron-cobalt alloys. 
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shows that Young’s modulus for iron increases with the 
addition of cobalt up to the composition of about 30 
percent cobalt, then decreases remarkably and is a mini- 
mum at about 85 percent cobalt, afterwards recovering to 
the value for pure cobalt. (Since, as mentioned above, 
the specimens employed in this measurement are made 
from cast rods, several of them have too low densities, 
due to the internal inhomogeneities. Consequently, calcu- 
lated values of Young’s modulus of such specimen rods 
have also been found to be too small, and points for them 
are omitted in the figure. The value for pure cobalt, as 
indicated by the double circle, is that obtained with the 
other specimen rod of electrolytic cobalt, prepared by 
hammering.) This result shows a good agreement with 
that of Honda and Tanaka and also with that of Nishi- 
yama® who measured in the low cobalt region with the 
static method of bending. It should be noticed, here, 
that our result and Nishiyama’s result do not show a 
minimum near 10 percent cobalt which was recognized 
by Honda and Tanaka. Further, the position of a minimum 
in the high cobalt region could not be determined as 
accurately by our measurement as by Honda and Tanaka’s 
measurement, because of the lack of measured points, but 
it would be needless to say that this minimum has a 
correlation with the phase boundary between a and y 
solid solutions which exists near 80 percent cobalt.* 

1K. Honda and T. Tanaka, Sci. Rep. Tohoku Imp. Univ. 15, 1 
a) Rep. Tohoku Imp. Univ. 27, 115 (1938). 

3H. Masumoto, Sci. Rep. Téhoku Imp. Univ. 18, 195 (1929). 

4 This value of the AE effect exceeds far beyond some 10 percent for 


pure nickel, which is, so far, the greatest value measured at room 
temperature. 

§Z. Nishiyama, Sci. Rep. Téhoku Imp. Univ. 18, 359 (1929). 

‘= Hansen, Der Aufbau der Zweistofflegierungen (Berlin, 1936), 
p. 485, 
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Cosmic Rays in the Bering Sea 


Paut F. Gast AND D. H. LouGHRIDGE 
University of Washington, Seatile, Washington 
April 10, 1941 


HE authors have previously reported! an investi- 

gation of the latitude effect between the magnetic 
latitudes of 53° 30’ N and 61° 36’ N with a precision type 
ionization chamber. The data there obtained showed the 
intensity, after reduction to a standard barometer and 
temperature, to be approximately independent of latitude. 
There were, however, some disturbing features ascribed to 
the “horizon effect,’’ the principal one being a drop of 
2 percent at the high latitude end of the curve. In order 
to check whether this was a purely local effect, and also 
in order to extend reliable measurements to as high a 
magnetic latitude as possible, measurements of cosmic-ray 
intensity were carried out along the northwest coast of 
North America from Seattle through the Gulf of Alaska 
and the Bering Sea and along the Alaskan Arctic coast as 
far as Point Barrow. 

The apparatus used was a Millikan-Neher type of 
electroscope® shielded with 11 cm of lead. Unfortunately 
there were two factors which materially reduced the 
precision of the present measurements over those of the 
preceding investigation.! The first factor was the lower 
intrinsic sensitivity of this instrument as compared with 
the Carnegie Model C* meter used previously and the 
second factor was the inability to take repeated readings 
over the same course, as had been done previously. The 
meter was mounted aboard the Coast Guard Cutter 
Itasca during the period from July 26 to October 10, 1940 
while it was on patrol duty in the Bering Sea. 

The hourly values of the cosmic-ray intensity were 
averaged over six hourly periods and reduced to a standard 
barometer of 760 mm of mercury by using a barometric 
coefficient of 0.17 percent per mm of Hg as determined by 
Millikan and Neher.? The daily averages of these intensities 
were then taken and plotted against the daily mean 
geomagnetic latitude of the ship. These results are shown 
in Fig. 1 where the vertical scale is arbitrary but propor- 


INTENSITY 


50° _|sse_ 60° 65° 70 


GEOMAGNETIC LATITUDE 


Fic. 1. Daily mean cosmic-ray intensities as a function of 
geomagnetic latitude. 
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tional to the intensity and the point 9.0 represents 1.90 
ions/cc/sec./atmos. of air. From the figure it can be seen 
that the values of intensity at any given magnetic latitude, 
which in general were taken at very different geographic 
locations, agree as well as the intensities taken on different 
days at a definite geographic location such as Unalaska 
or Nome. This would indicate that there is no significant 
difference in intensity over the whole region investigated. 
The results were next averaged by degrees of magnetic 
latitude and all reduced to a standard atmospheric temper- 
ature of 50°F by the use of the temperature coefficient of 
—0.05 percent per degree previously found by the authors! 
and other investigators‘ for rays of the type measured by 
such a lead-shielded ionization chamber. The temperature 
correction made no significant difference in the intensities 
with the exception of the intensity measured at Point 
Barrow where the temperature was 36°F. 

The mean intensity vs. magnetic latitude curve after 
reduction to standard temperature is shown in Fig. 2 


INTENSITY 
4 


160° 
GEOMAGNETIC LATITUDE 


Fic. 2. Mean cosmic-ray intensities for each degree of latitude plotted 
against latitude. 


where the vertical scale shows percent deviation from the 
mean. Since, under the conditions of measurement, the 
intensities are accurate to within 1 or 2 percent,? the 
intensity can be regarded as substantially independent of 
magnetic latitude over the region investigated. 

The results of the present investigation are two. First, 
it is shown that the drop at the high latitude end of the 
intensity vs. latitude curve previously obtained by the 
authors is not the beginning of a downward trend but a 
purely local phenomenon due, evidently, to the “horizon 
effect.” Secondly, the investigation of the latitude effect 
has been carried to higher magnetic latitudes (68° 22’ N) 
than have previously been reached and the intensity found 
to be independent of latitude up to this point. The accuracy 
of the measurements was not such as to afford any reliable 
information on the “atmospheric latitude effect’’ of 
Compton and Turner® nor of the allied “air mass effect’’* 
in these regions. 

The authors wish to express their sincere appreciation 
to Commander F. A. Zeusler, Commandant of the Juneau 
District, U. S. Coast Guard, and to Commander L. W. 
Perkins, commanding officer of the cutter Jtasca, for 
making possible a trip into the Bering Sea and Arctic; 
to President L. P. Sieg of the University of Washington 
for making available the necessary funds; to Dr. J. A. 
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Fleming of the Department of Terrestrial Magnetism for 
supplying the apparatus, and to Dr. H. V. Neher for 
carefully rechecking the quartz electroscope. 


1 Paul F. Gast and D. H. Loughridge, Phys. Rev. 59, 127 (1941). 

2 R. A. Millikan and H. V. Neher, Phys. Rev. 50, 15 (1936). 

3A. H. Compton, E. O. Woldan and R. D. Bennett, Rev. Sci. Inst. 
5, 415 (1934). 

4V. F. Hess, Phys. Rev. 57, 781 (1940). 

5 A. H. Compton and R. N. Turner, Phys. Rev. 52, 799 (1937). 

6D. H. Loughridge and Paul F. Gast, Phys. Rev. 58, 583 (1940). 


The Magnetic Anisotropy of Gadolinium 
Sulphate Octahydrate 


K. S. KRISHNAN AND S. BANERJEE 
Indian Association for the Cultivation of Science, Calcutta, India 
March 11, 1941 


N a previous publication! we reported preliminary 

measurements on the magnetic anisotropy of gado- 
linium sulphate octahydrate, Gd2(SO,)3-8 HO, made with 
a moderately pure specimen. The difference between the 
two extreme principal susceptibilities of the crystal was 
found to be about 1.4 percent of its mean susceptibility. 
This value for the magnetic anisotropy, though small, is 
much larger than should be expected from the narrow 
separation of the 8S levels of Gd*** that would occur 
under the crystalline electric fields in the neighborhood of 
the Gd*** ions in the crystal. 

Recently we have made measurements with a specimen 
of known purity, kindly presented to us by Professor 
Trombe, of the Paris University. The specimen was quite 
free from both samarium and terbium, and its europium 
content was less than 0.1 percent. The main impurity 
present was the diamagnetic yttrium salt, about 1.4 
percent, which will not affect the magnetic anisotropy 
appreciably. The results of the magnetic measurements are 
given below. 

The crystal of Gd2(SO,)3-8 H2O is monoclinic, and has 
the axial elements a : 6 : c=3.009 : 1 : 2.007, 8B=118°.0. 
Denoting the maximum and the minimum susceptibilities 
in the (010) plane, for the above formula weight (747 
grams, containing 2 gram ions of Gd***), by xi and x2, 
respectively, and the susceptibility along the 5 axis by xs, 
it is found that at room temperature (303°K): xi—x2 
= 36X 10-8; xs— x1 = 16 10-8; x = (x1 +x2+x:3)/3 =52,000 
X 108 e.m.u. The x; axis makes an angle of about 17° with 
c and 45° with a. The difference between the extreme 
susceptibilities, namely x:s—x2=52X10-§, is only 0.1 
percent of the mean susceptibility. This value of the anisot- 
ropy will correspond to a separation of the order of « 
between the adjacent levels in the Stark pattern of the §S 
state of Gd*+**, where 

or «~kX0.3°, or 0.2 
which is of the same order as should be expected from the 
demagnetization and the specific heat measurements on 
the salt at very low temperatures.? 


1 Krishnan, Mookherji and Bose, Phil. Trans. Roy. Soc. A238, 133 


(1939). 
2 See M. H. Hebb and E. M. Purcell, Phys. Rev. 51, 384A (1937). 


770 LETTERS TO THE EDITOR 


A Single Component for the Primary 
Cosmic Radiation 
W. F. G. Swann 


Bartol Research Foundation of The Franklin Institute, Whittier Place, 
Swarthmore, Pennsylvania 


April 17, 1941 


Wwe reference to the Letter to the Editor entitled 
“The nature of the primary cosmic radiation and 
the origin of the mesotron,” by M. Schein, W. P. Jesse 
and E. O. Wollan,' I wish to point out that identical 
conclusions have been cited by the present writer from 
other considerations. These conclusions, which are covered 
in three published communications,’ are as follows: 

(1) There is only one type of primary radiation, a charged 
particle radiation—probably protons—comprising _ par- 
ticles of heavy mass. 

(2) By processes at present unknown, the primary 
radiation gives birth, probably indirectly, in the upper 
atmosphere, to mesotrons. 

(3) Those mesotrons which are born approximately at 
rest will have such short lives that they will disintegrate 
before they have traveled more than 300 meters. They 
will, in fact, disintegrate in the stratosphere, and in so 
disintegrating, will give rise to electrons which, on account 
of the disintegration occurring from mesotrons at rest, will 
emerge on the average equally in all directions. 

(4) The mesotrons formed with higher energy will 
disintegrate at lower altitudes, because of their longer 
lives, and because they disintegrate at high energy, will give 
rise to electrons which possess on the average a forward 
component at these lower altitudes. 

(5) The assumptions (1) and (2) lead, through the 
logical consequences (3) and (4), to an explanation of the 
following facts: (a) As shown by the Bartol Foundation’s 
observations in the two National Geographic—U. S. 
Army Air Corps Stratosphere Flights** and in the Jean 
Piccard Flight,5 the curve of intensity versus zenith angle 
flattens out with increasing altitude to a condition in which, 
at a depth of 0.5 meter in the water equivalent atmosphere, 
there is only a 20 percent change from vertical intensity 
to horizontal intensity. (b) As shown by our stratophere 
observations in Explorer II* and independently by the 
observations of T. H. Johnson and J. G. Barry,® there is 
no appreciable azimuthal asymmetry at high altitudes. It 
was, of course, to provide for the experimental facts cited 
under (a) and (b), that the hypotheses (1) and (2), with 
the consequences (3) and (4), were formulated. 

(6) Any incoming electrons of primary origin and of one 
sign would necessitate azimuthal asymmetry and primary 
electrons of one sign or of both signs would cause strong 
dependence of intensity upon zenith angle. The absence 
of such effects prohibits assumption of the existence of 
any primary electrons and necessitates that all the cosmic-ray 
electrons be born from the mesotrons in accordance with (2)." 
In other words, there can be no primary electrons. 

That the primary cosmic rays are in part protons has 
been suggested by A. H. Compton and H. Bethe,’ and also 
by T. H. Johnson.* A single primary radiation of protonic 
nature giving rise to mesotrons and through them to the 
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cosmic-ray electrons was suggested by the present writer! 
prior to his more detailed development of these ideas in 
papers cited. T. H. Johnson" also suggested protons as 
the parents of the mesotrons but invoked a primary 
electronic radiation in addition. 


4 Schein, W. P. Jesse, and E. O. Wollan, Phys. Rev. 59, 615 
1941 
¢ 2W. F. G. Swann, Phys. Rev. 56, 209 (1939); Rev. Med. Phys. 11, 
242 or. nate in particular, pp. 251-254.) Phys. Rev. 58, 200 (1940). 
. Swann and G. L. Locher, Nat. Geog. Soc. Contr. Tech. 

Pap. Series No. 1, pp. 7-14 (1935). 

«W. F. G. Swann, G. L. Lou and W. E. Danforth, Nat. Geog. 
Soc. Contr..Tech. Pap., Stratosphere Series No. 2, pp. 16-25 (1936). 

a. Swann, J. Frank. Inst. 222, 1 (1936). (See, in particular, 

4 
PP; W. F. G. Swann, Phys. Rev. 56, 209 (1939). 

7And to a negligible extent by shower production, which phe- 
nomenon, however, becomes relatively more important in the lower 
atmosphere. 

8A. H. Compton and H. Bathe, Sa 134, 734 (1934); also A. H. 
Compton, Rev. Sci. Inst. 7, 71 (19 rr 

*T. H. Johnson, Phys. Rev. 45, 569 (1934) 

Frank. Inst. 226, 757 1938). (See, in particular, 
pp. an 

 T. H. Johnson, Rev. Mod. Phys. 11, 208 (1939). 


Vibrational Structure of the ‘=,~(O.*)—*z,~ 
Rydberg Series of 
Y. TANAKA AND T. TAKAMINE 


Department of Physics, University of California, Berkeley, California 
March 24, 1941 


HE far ultraviolet absorption spectrum of molecular 

oxygen was studied by Price and Collins! in 1935. 
Our recent study in Berkeley, while confirming most of 
their observations, leads us to propose certain changes in 
their classification. The experimental arrangement was 
the same as in the case of N2 reported recently,? except 
that the pressure of oxygen in the spectrograph was 0.001 
mm. 

Especially significant are the features we find in the 
region between \740A and 660A. Stated briefly, it appears 
that the bands designated as the Q progression, by Price 
and Collins together with some of their V and W bands,’ 
can be arranged into three different progressions, each 
forming a Rydberg series and converging respectively to 
the normal (v=0), v=1, and v=2 vibrational states of an 
excited O,* electronic state which Mulliken‘ tentatively 
called b42,-, thus constituting the (0,0) (1,0) and (2,0) 
bands of a *2,-+*S,~ Rydberg series of Oz. 

If we assign the three strong bands at \732.22A (int. 4) 
725.98A (4) and 720.07A (3) as the (0,0) (1,0) and (2,0) 
bands for n=3,5 thirteen members (n=3 to 15) were ob- 
served for each series. The frequencies of the bands in 
each of the three series follow a Rydberg law closely, 
giving as the series limit voo= 146,548 cm™ for (0,0) band, 
147,705 cm for (1,0), and 148,831 cm for (2,0). The 
successive differences between the above numbers, namely 
AG, = 1157 cm™, and AG; +4= 1126 cm™ agree closely with 
the values 1162 cm™, and 1127 cm™ obtained from the 
known constants for the 6 ‘2,~ state of O,*. 

The ionization potential corresponding to the O,* 
b‘Z,- normal state comes out as 18.08+0.01 volts, which 
is close to the value 18.1+0.1 volts obtained by combining 
the data® for the )‘=,-—+a‘II, bands of O2* with the 
ionization potential of the a ‘II, level given by Tate and 
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Smith and Hagstrum and Tate’ from electron impact 
experiments. 

The value obtained by Price and Collins for the b *2,~ 
state is 18.2+0.1 volts, which would seem to Powetetiee. 1 
to the mean of our values for the heights of the normal, 
first, and second vibrational states of b ‘Z,~. 

Since there has heretofore been no spectroscopic con- 
nection between the quartet and doublet states of O2*, 
our result would be of significance in supplying the accurate 
spectroscopic data for the energy of the quartet levels. 

Unlike the case of N2 that we reported in these columns, 
the (1,0) and (2,0) bands in O: are even stronger than the 
(0,0) band. Those bands that lie just beyond the limits of 
the series appear very diffuse, likely because of auto- 
ionization. 

A more detailed report will appear soon in the Scientific 
Papers of the Institute of Physical and Chemical Research 
(Tokyo). 

We wish to express our sincere thanks to Professor F. A. 
Jenkins for his continued interest. 

1W. C. Price and G. Collins, Phys. Rev. 48, 714 (1935). 

2? Y. Tanaka and T. Takamine, Phys. Rev. 59, 613 (1941). 


3 In these groups Price and Collins mention that their classification 
is uncertain. 

*R.S. Mulliken, Rev. Mod. Phys. 4, 56 (1932). 

5 Our assignment is tentative. 

*R. Frerichs, Zeits. f. Physik 35, 683 (1926); T. E. Nevin, Trans. 
Roy. Soc. 237, 471 (1938). 

3. Ee Tate o—e P. T. Smith, Phys. Rev. 39, 270 (1932); H. D. 
Hagstrum and J. T. Tate, Phys. Rev. 59, 354 (1941). 


Variations in the Relative Abundance of the 
Carbon Isotopes 


Byron F. MURPHEY AND ALFRED O. NIER 
Department of Physics, University of Minnesota, Minneapolis, Minnesota 
April 18, 1941 


EVERAL years ago it was shown! that variations are 
present in the C#/C™ abundance ratio of carbon 
obtained from various natural sources. The results have 
now been confirmed and extended in a new investigation 
with the 60° mass spectrometer? of 57 specimens of 
widely different geological age and geographic origin. 

The samples were all converted to CO, and the C#/C" 
ratio determined from the relative heights of the mass 44 
and 45 peaks after correction had been made for C#O0"*O" 
molecules.’ 

As the measurements were made over a period of several 
months a control sample was included each time analyses 
were made. The difference between the control and the 
sample in question was then noted. All computed ratios 
are based on these differences and the assumed constancy 
of the control. Frequent interchecks were made. This 
method eliminates possible errors due to variations in the 
mass spectrometer. 

In Table I are recorded the results obtained. Each 
sample was analyzed at two different times and the final 
number adopted was the average of the two separate deter- 
minations. In over ninety percent of the cases the two 
separate values found for the C"/C" ratio agreed within 0.4 
percent. In only three cases was the agreement worse than 
0.6 percent, the poorest being 0.9 percent. The differences 
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TaABLeE I. C!2/C'3 ratio for various carbon sources. 


RANGE OF 
VALUES OF Av. RaTIO 
No. OF C®2/C8 FoR SAMPLES 
SOURCE SAMPLES FOR GROUP IN GROUP 

Limestone 10 88.8-89.4 89.2 
Coal (humic origin) 10 91.3-92.0 91.8 
Wood 7 91.6-92.2 91.8 
Petroleum 6 92.0-92.8 92.5 
Bituminous shales 7 92.1-92.7 92.5 
Torbanite and kerosene shales 3 91.3-92.0 91.7 
Meteoric carbon 7 89.8-92.0 91.3 
Graphite 1 90.2 
Zeolitic calcite 1 89.9 
COs in air 1 91.5 
Minneapolis, Minn., 3/1/41 
Lycopodium spores 1 93.1 
‘*Balkashite"’ algae 1 92.8 
Marine shell 1 89.5 
Sea water 1 89.3 


are thus significant although all of the values may be 
too high or too low owing to constant discrimination in 
the spectrometer. This constant error in the absolute 
numbers is probably less than 1 percent, although it may 
be as great as 2 percent. 

A glance at the table shows that variations up to 5 
percent exist in the C'/C ratio. As was shown in the 
earlier work, C™ is concentrated in the limestones whereas 
plant forms have a preference for light carbon. There 
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appears to be no “‘age effect.’’ The limestones varied in 
age from Pre-Cambrian to modern, the coal from carbon- 
iferous to Pliocene, the wood from Late Pleistocene to 
present. Although there are individual variations within 
the groups no trends with age were noted. 

It is interesting to note that the values for sea water 
and the marine shell check closely with those for the 
limestones. As sea water is considered to be the regulator 
of atmospheric CO:, it would have been interesting to 
analyze CO, collected above the sea. The high ratio found 
for the sample collected in Minneapolis can, no doubt, be 
attributed to the CO», produced when coal is burned. 

We are indebted to Drs. H. Berman and W. C. Darrah 
of the Harvard University Minerological and Botanical 
Museums for many of the samples. The meteor samples 
were kindly given to us by Professor G. P. Baxter of the 
Harvard Chemistry Department. The bituminous shales, 
plant spores and algae were supplied by Dr. Taisia 
Stadnichenko of the U. S. Geological Survey. Mr. R. O. 
Belkengren of the Botany Department kindly permitted 
the use of his Van Slykes apparatus for converting samples 
to CO». Through his efforts the old wood samples were 
obtained. 

1A. O. Nier and E. A. Gulbransen, J. in. Chem. Soc. 61, 697 (1939) 


2A. O. Nier, Rev. Sci. Inst. 11, 212 (1940). 
3B. F. Murphey, Phys. Rev. 59, 520 (1941). 
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HE second meeting of the Metropolitan Section of the American Physical 
Society for the season 1940-1941 was held at 3:00 P.M. on Friday, 
March 28, 1941, in the Pupin Physics Laboratories of Columbia University. 


The following papers were presented: 


The Structure of Black Carbon. A. H. WHITE. 
Microphonic Contacts and the Properties of Pyrolitic Carbon. R. O. 


GRISDALE. 


Quantum Relationships in Vision. SELIG HECHT, SIMON SHLAER, AND 


H. PIRENNE. 


The following officers were elected for the season 1941-1942: 


Chairman, BERNHARD KURRELMEYER 

Vice Chairman, RICHARD T. Cox 

Secretary- Treasurer, W. S. GORTON 

Members of the Executive Committee, J]. M. B. KELLOGG 


Davip B. LANGMUIR 
W. S. GorRTON 
Secretary- Treasurer 
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